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The more exact dissociation theory of Part II, is applied to eighteen different gases and 
vapors, ranging in molecular weight from Hz to Xe. Interaction energies per molecule in 
clusters formed during condensation of the saturated vapors are found as functions of tem- 
perature. From the differences between surface energy and measured surface tension, which 
remain finite at 7,, estimates are made of the numbers of heavy clusters maintained in the 
saturated vapor over the liquid surface, also as functions of temperature. In the case of eight 
vapors, the theory predicts measurable differences between the two critical points, T, at which 
the meniscus blurs, and Ty at which the densities equalize; the largest differences being for 
water, 6.5°K; methyl alcohol, 5.1°K; and nitrous oxide, 4.1°K. True critical volumes at Ty 


are calculated. 


INTRODUCTION 


N Part III,! the simple theory of Part I has 
been applied to the vapors of COs, 
and CH;Cl. The theoretical formula for PV/NkT 
of the saturated vapors was assumed to agree 
approximately with the facts, even up to the 
critical temperature, and from this assumed 
agreement, curves were obtained for the surface 
energies against temperature. As was to be ex- 
pected, these curves agreed well with surface 
tension data only for temperatures sufficiently 
below T,. Actually, however, the range of tem- 
peratures showing agreement should have been 
higher and the discrepancies appearing between 
0.57, and about 0.97, required special interpre- 
tation. It was there suggested that the clusters 
involved in the liquid vapor equilibrium were 
actually small clusters, and that, therefore, the 


1 Previous papers in this series, hereafter ree to as 


Part I, Part II, and Part III, respectively: W. Band, J. 
Chem. ’Phys. 7, 324-326 (1939); J. Chem. Phys. 
7, 927-931 (1939) ; T. P. Tseng . Ch’eng 
and W. Band, J. Chem. Phys. ‘oad 


surface energy computed by the theory referred 
rather to these small clusters than to the true 
surface tension of the liquid phase. A very rough 
estimate was given of the size of the clusters 
required to account for the observed discrep- 
ancies from this point of view. In the light of 
later work to be reported in the present paper, 
a much better explanation has been developed, 
and it now seems clear that the improved theory 
of Part II is sufficiently accurate to give true 
surface energy terms right up to the critical point 
without any such special assumption. 

In the simple theory of the saturation condi- 
tions of a vapor, Part I, the internal energy W, 
of a cluster of type s was expressed in the form, 


1(10) 
W.=(s—1)W+x(s, T) (1) 


because the part of W, linear in (s—1) controlled 
the convergence of the series of terms which gave 
the concentrations of the several types of cluster, 
see 1(13). The sum of the series at saturation 
then depended only on the nonlinear part 
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x(s, T), written briefly as x in I(16); finally the 
linear part, W, was eliminated from the theo- 
retical expression for PV/NkT at the saturated 
state, 1(17). 

In the second approximation which includes 
cluster sizes, Part II, it remains necessary to sep- 
arate W, in the above manner, and for the same 
reasons; but the term W is no longer eliminated 
exactly from the expression for P(V—V1)/NkT 
at saturation, see II(23) and II(9). 

On the other hand, no such separation of W, is 
required for the equations of an imperfect gas, 
see 11(38) and I1(39), where W+ x2= Ws. 

The separation of W, is essentially a mathe- 
matical device required for the handling of a 
purely mathematical problem arising in the study 
of saturated conditions of a vapor; but it has an 
obvious physical significance. In Part III, (6) 
it was found convenient to write 


x(s, T)=aA(T)(s*—1), (2) 


and the quantity aA(T) was evidently the sur- 
face tension per molecule of the liquid as or- 


dinarily defined for large clusters. For small 


clusters there is no definition of surface tension, 
and formula (2) is to be regarded simply as a 
definition of x, in terms of the surface tension 
per molecule, aA(T) of large clusters. It may be 
regarded as an acceptable definition of the sur- 
face energy of small clusters, but, in fact, this 
idea will not be used here. 

Now it is physically reasonable that the inter- 
action energy W2 between pairs should not be 
separated in the above treatment of an imperfect 


TABLE I. Molecular interaction energy at critical point. W 2(a) 
from Eq. (5); (b) from Eq. (6); (c) from Eq. (7) and 
(16), (d) from Eq. (29); units are 10-4 erg. 


P, We 
Gas Te °K | ATMOS Ve (a) (d) 
He 33.2 12.80 64.92} 0.88) 0.94) — _ 
He 5.2 2.26 58.13} 0.01) 0.02) — _ 
NH; 406.0 | 115.0 72.94) 77.1 | 52.0] 78. 58. 
653.5 | 217.7 51.6 | 87. 84. | 119. | 106. 
Ne 44.4) 26.9 42.2 3.3 3.4 2.4) 2.9 
Ne 126. 33.5 91.0 | 14.3 | 14.3 8.7) 14. 
154. 49.7 75.0 | 18.1 | 18.3 | 25 15. 
CH;OH 513.1 79.7 117.7 | 99. 77. | 101. | 166. 
155.6| 48.0 75.8 | 184) 185) — 
N:0 309.6 77.5 96.9 | 44.0] 443 | 58. 62.8 
304.2 73.0 94.8 | 43.0 | 43.4] 57. 10.6 
C:H;OH 516.2 62.7 167.0 | 86.6 | 83.7 | 107. | 123. 
aC 416.3 66.7 138.0 | 91.6 | 65.4 | 83.7| 35. 
CH;COOH | 594.7 58.0 171.2 | 142. | 102. | 131. | 196. 
C:H:0H 264.0] 50.6 219.7 | 56.0 | 41.6] 69. | 90. 
2(C2Hs)O 194.0} 36.0 282.3 | 40.8 | 30.8] 80. | 55 
Hg 1500.? | 203.? 50.? | 390.? | 300.2 | — | 230.? 
Xe 287.8 57.2 113.6 | 48.1 48.2) — _ 


gas, where there is no liquid phase; and equally 
reasonable that the surface term should play a 
distinct part, as it is found to do, in the behavior 
of the saturated vapor where a surface must 
exist. 

Now the critical point conditions can be ap- 
proached either from the equations of the satu- 
rated vapor or from those of the imperfect gas, 
and since the surface energy occurs in the one 
approach, and not in the other, a comparison of 
the two at once leads, as will be shown below, to 
conditions on x(s, T), or rather on x2=x(2, T-) 
at the critical temperature. 


MOLECULAR INTERACTION ENERGY AND 
SURFACE ENERGY AT 7, 


In Part II (38), (39) the equation of an im- 
perfect gas was obtained approximately in the 
form of Dieterici: 

P(V—B)=NhkT exp (—A/NkTV), 


where 
A 


and B is a function of T not given very precisely 


by the theory. If B and AT? are approximately 
constant near the critical point, the latter is 
given by the usual equations: 


V.=2B, P.=A/4e?B*, T.=A/4NkB. (4) 
The function B can be eliminated : 

= exp (W2/kT.), (5) 

exp (—Ws/kT.). (6) 


Values of W:2 obtained from (5) and (6) using the 
experimental values of P., V. and T, are shown 
in Table I. 

Since the equations of the saturated vapor are 
extended up to T, only formally, they lead to no 
perfectly clear and definite condition for that 
point. Three different ways of applying the 
formal relations will be given leading to results 
which do not differ too badly from each other. 

Elimination of (V— V1) from between II(21) 
and II(23) for the saturated vapor leads to 


(7) 
Comparison with (5) shows that 
exp (x2/kT.) = (e?/2)(1+¢0e), (8) 


and 
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which determines both x2/kT, and ¢o-: 
x2/kT.=2.075, ¢0-=0.429. (9) 


On the other hand I1(23) can itself give a critical 
point at which V= V;, namely 


Bod = (10) 
which is consistent with (6) only if 
= V. exp (x2/kT-). 
But Eq. II(16) gives V.=Bod/(1+¢3-) and hence 
exp (x2/kT.) (11) 
Again this determines both x2/kT, and ¢3-: 


X2/kT.=1.666, 3-=2.754. (12) 


Finally (9) and (11) lead to the relation 
(1+¢3) = (e?/v2)(1+¢.), 


which determines both sums: 
doc=4.1, 3-=16.2, x2/kT.=3.2. 


The general inference would be that at the 
critical temperature the term x2/kT. should be 
roughly the same for every substance, somewhere 
between, say, about 1.5 and 3.5. Any attempt to 
make a zero value satisfy the various equations 
results in violent disagreements. The circum- 
stance that the theory demands a positive value 
of x2 at T, indicates, through (2), that aA(T) 
and, therefore, the surface energy of large clusters 
does not vanish at 7.. 


(13) 


CURVES OF INTERACTION ENERGY vs. 
TEMPERATURE 


From Eq. (8) it is possible to find Wasa func- 
tion of T right up to 7; if ¢o can first be found 
throughout this range. This is done by means of 
the data on P(V—V,)/NkT and the Eq. II(23). 
From Part II, Eqs. (16) and (22) we have: 


V.= Np/(1 +3) 
Ve=(V./ V1) (Bob) /(1+¢s), 
so that 
and II(23) simplifies to 
P(V—V;)/NRT =(1+¢0)/(1+¢s) 


and 


(15) 
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This relation is at once amenable to computation, 
leading from the data on P(V—V.)/NkT and 
(V:/V-.) to the temperature dependence of $9 and 
$3. The critical point occurs in (15) formally as 
the point at which (V;/V.) is actually unity, and 
the whole expression zero. In view of the known 
discrepancy at 7., and the existence of T, above 
T., the critical point in (15) will be interpreted 
as T, rather than 7., so that the equation will 


TABLE II. Revised values of fe $s, (1+¢0)/(1+¢s3) in terms 
of w. 


(1+¢0) (1 +¢0) 
w (1+¢3) || w os (1+¢3) 
0.5 | 145.2 3311. 0.044 5 |0.1783 0.3913 0.847 
1.0} 16.21 141.2 0.121 6 |0.09231 | 0.1947 0.914 
1.5 5.391 27.17 0.227 7 |0.04924 | 0.1016 0.953 
2.0 2.482 9.247 | 0.340 8 |0.02671 | 0.05443 0.974 
2.5 1.368 4.215 | 0.454 9 |0.01463 | 0.02959 0.986 
3.0 0.8337 2.269 | 0.561 10 | 0.008058 | 0.01623 0.992 
3.5 0.5401 1.3510} 0.655 11 | 0.004456 | 0.008949 0.996 
4.0 0.3644 0.8588] 0.734 | — 
* w=a@A(T)/kT. 
be written as: 
(16) 


The direct summation of ¢o and ¢3 for smaller 
values of aA(T) has been carried out, and Table I 
of Part III improved. These results are shown in 
Table II. Tables were next constructed showing 
the quantities on the right side of (16) as func- 
tions of aA(T) and (V:/V,) for several choices 
of ¢3,. These tables then directly gave aA(T) at 
all temperatures at which data were available 
both on P(V—V;)/NkT and (V:/V,), for any 
choice of ¢3,. Actually Vy is unknown in most 
cases, but then never differs much from V,, which 
was, therefore, used instead of V, as a first 
approximation. 

The actual choice of $3, is somewhat arbitrary, 
and was made in the following manner. In most 
cases the data on P(V—V,)/NkT and (Vi/V.) 
show sudden drops to zero and unity, respec- 
tively, at T., and, in the light of Callendar’s work 
on H,0,? this discontinuity is probably spurious. 
Therefore, the two quantities were determined 
at T, by a very slight extrapolation which ignored 
this discontinuity in both curves, adjusting the 


projected curves so that they would both arrive 


2 Callendar, Proc. Roy. Soc. A120, 460 (1928). 
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Fic. 1. Mean molecular interaction energies in clusters 
formed in saturated vapors. 1—CH;COOH; 2.—CH:;Cl; 
3.—CH;0H; 4.—C.H; H; 5.—N;0; 6.—COs:; 7.—NH;; 
8.—H:0; 9.—O.; 10.—N2; 11.—Ne. Abscissae, T/Ty; 
ordinates, — W-10" erg. 


at their respective zero and unit values at the 
same temperature. If this extrapolation gives 
P(V—V;)/NkT=J, and (Vi/V.)=K at T., 
these values are inserted in (16). The values of 
¢o and ¢3 at T, are arbitrarily assumed to be those 
given by (13) for all substances without attempt- 
ing to adjust them to individual cases. This then 
gives 

(1+¢3,) =4.74K/(0.276—J). (17) 


This determines ¢3, which, in turn fixes the, 
whole curve of Eq. (16) as applied to the de- 
termination of aA(T). The best way to find T, 
then appears to be an extrapolation of the 
curve of aA(T) to its value corresponding to 
¢3y. As a check, the previous extrapolations of 
P(V—Vi)/NkT and (V;/V.) are available. It 
must be admitted that this process has been very 
arbitrary, but it is quite definite, and any other 
reasonable way of applying the formal relation 
(16) would lead to substantially the same results. 
Having thus secured ¢o as a function of T, W 
at once follows from (17) as mentioned above. 
The value of x2 having already been fixed at 
3.2kT., this results in an estimate of W2=W+x2, 
at 7.. These estimates are included in Table I. 
They compare fairly satisfactorily with those 


found from the imperfect gas equations. The . 


actual discrepancies might be reduced by a more 
careful choice of ¢3- and hence of $3, for in- 
dividual cases. 

Figure 1 shows curves of W against T/T,, for 
some of the substances appearing in Table I. At 


first sight it is surprising that the energy increases 
in magnitude with increasing temperature, but 
this is due to the increasing density of the 
saturated vapor towards T,. There is most de- 
cidedly no tendency for W to vanish either at T, 
or even at 7). 


REDETERMINATIONS AT THE CRITICAL POINTS 


In a few cases, the discrepancy between T;, 
and T, leads to significant changes in accepted 
critical point data. In these cases let V. continue 
to represent the mean volume at 7., as computed 
by the law of rectilinear diameters: 


1/V+1/Vi=2/V.. (18) 
Also let J be the extrapolated value at T, of 
P(V—V)/NkT.=J. (19) 


From (18) and (19) both true volumes V and V; 
at T. can be found. The true critical volume V, 
at TJ, can then be found by extending the recti- 
linear diameter to T,: 


2/Vy=a+0T, (20) 


where a and b are found from the data on V and 
V1 expressed in the form, near T,, 


1/V+1/Vi=a+0T. 


The actual results of these calculations are shown 
in Table III. In this table the extrapolated value 
J is given, and the corresponding value of ¢3,. 
V,. is the true liquid volume at 7,, and V, is the 
ordinarily accepted mean volume at TZ,, in 
cm*/mole. 


APPLICATION OF NERNST’Ss HEAT THEOREM 


Either by using the theoretical expression, ob- 
tained below, for the latent heat L, or simply by 
assuming Clapeyron’s relation, ZL can be found 
in terms of the temperature derivative of 


TABLE III. Critical point values at T, and Ty. 


Gas oy J Te °K | Ty °K Ve Vie Vy 
H20 24 0.114 | 647.0 | 653.5 51.6 42.1 47.8 
Ne 170. 0.252 44.4 46.9 41.8 30.7 44.2 

35 0.172 | 126.0 | 128.5 90.2 70.8 | 91.3 
CH;:0H 60.7 | 0.213 | 513.1 | 518.0 | 117.7 60.7 | 119.4 
37.1 | 0.181 | 309.6 | 313.7 96.9 65.5 | 98.7 
CO2z 26.9 | 0.133 | 304.2 | 304.6 94.8 72.1 94.9 
C:H;OH 24.4 | 0.114 | 516.2 | 518.0 | 166.6 | 128.9 | 167.9 
CHsCOOH | 22.3 | 0.114 | 594.7 | 598.0 | 171.2 | 122.6 | 172.2 
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(W/kT). Then assuming Nernst’s heat theorem 
this equation can be integrated to give W without 
the arbitrary factor of T, in terms of the data 
on L. This process is by no means precise, if only 
because the function can hardly be assumed con- 
tinuous down to 0°K, through the various phase 
changes, but it can be regarded as a rough check 
on the other methods of finding W. 
Thus, the total energy of the vapor assembly is 
E=NY wkT*(0/0T) In f(T), (21) 
s=1 
and if we use the expression (2) of Part II for 
f.(T) this reduces to 


E=N ¥ u.kT?(0/0T) In (1—8,) 


(22) 


Neglecting the temperature variation of 8,, and 
the contribution of the surface energies of the 
clusters in the vapor, this becomes 


E=(3)NRT us 


s=1 


+N 


s=1 


(23) 


In the liquid phase, again neglecting the surface 
energy, the total energy, to the same approxi- 
mation, is 


(24) 


so that the total latent heat is approximately 


L=E-E,+P(V—V;) 
=[P(V—V1)/NRT][SNRT/2 


—N(W-TaWw/aT)]. (25) 


The procedure now is to take experimental values 
of L and of P(V—V:)/NkT and to write 


W—TdW/dT =5NkT/2 
= F(T), say. (26) 


Then F(T) is expanded in the form C+DT?+--- 
and, using Nernst’s theorem, integrated at once 
to give 


W=C—DT?-::-. (27) 


The coefficients C and D for several substances 
were found in this way and are shown in Table 
IV. The values of W2=W-+ x2 at T, found from 
these values appear in Table I. Agreement is not 


Fic. 2. Mean numbers per cm* of heavy clusters main- 
tained over liquid surfaces in saturated vapors. Numbers 
identical with those in Fig. 1. The ordinate scale for Ne 
a reduced by one-half. Abscissae, T/T; ordinates, 
10°, 


very precise, but is much better in those cases 
where T, is low. 


SURFACE ENERGY AND SURFACE TENSION 


As described above, ¢o was found as a function 
of T. From Table II this means that by using (7) 
and (8) of Part III, the surface energy function 
aA(T) is known. As in (21) of Part III the true 
surface tension should be given by 


E(T) =(1/4.837) aA (T), (28) 
where », is the number of molecules per cm* in 


TABLE IV. A pplication of Nernst’s theorem to W=C—DT?. 
C in ergs, D in ergs/(°K)?. 


Gas —C-1018| D-1019 Gas -—C-10%3 | 
NH; 3.9 1.2 CH;OH 6.3 39.2 
H.0 7.30 7.81 | C:H;,OH 8.0 16.4 
Ne 0.23 0.30 | CH;Cl 2.75 4.75 
Ne 0.90 33.0 CH;COOH | 6.8 36.3 
O. 1.0 0.21 | C;H;OH 8.5 17.0 
COz 0.17 9.6 2(C2Hs)O 4.0 6.9 


the liquid. Actually E(T), as given by (28), is 
always greater than the measured surface ten- 
sion, remaining finite at T,. This has already been 
explained as due to the gravitational settling 
down of the larger clusters over the liquid surface. 
Thus suppose the mean size of these clusters is 
k molecules, and their mean number per cm? over 
the liquid is m,. Then the surface tension actually 
measurable will be 


(29) 


S(T) =E(T) +m {RW+aA (T)R3}. 
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It is certain that k must be at least of the order 
107, so the term in k? will be negligible. Then the 
product 7; can be found from the relation: 


—nkW=E(T)—S(T). (30) 


The interaction energy — W having already been 
determined, either as in Fig. 1 or in Table IV, 
this equation at once gives ;,k as a function of 
temperature. Assuming the mean value of k=10’, 
and W from Fig. 1, the results are shown in 
Fig. 2. 

It is specially remarkable that the numbers 
tend to increase with decreasing temperature, 
and that they are generally larger for those sub- 
stances showing marked difference between 7, 
and 7,. To fully appreciate these results it is 
necessary to abandon naive conceptions of a 
static interface between liquid and vapor, and, 
with Adam? realize that, from a microscopic 
point of view, the surface is in violent motion and 
almost never really ‘“‘flat.’”’ The present results 
very strongly suggest that the exchange across 
the boundary does not consist only of single 
molecules, but that the motion is so violent that 
whole clusters of molecules are also being torn 
from the surface in sufficient numbers to main- 
tain, against gravitational sedimentation, a per- 
manent layer over the surface of the liquid. 


PRECISION OF THE THEORY 


The theory lacks precision at two points. First 
in the arbitrary choice of ¢3.. To check this, the 
data on COz were analyzed, using values of ¢3. 


3N. K. Adam, The Physics and Chemistry of Surfaces 
(Oxford, second edition, 1938), pp. 6, 7. 


ranging from 2.3 to 5000; the values of km 
changed by a factor of five, while those of W by 
about 40 percent. This might thus be used to get 
a better agreement between the different parts of 
the theory, although too large a value of $3 
would be forbidden by Eqs. (9), (12) and (13). 

The other point is more fundamental and has 
already been emphasized by Mayer.‘ The cluster 
is here defined in terms of interaction energy 
between pairs, and it is necessary to fix arbi- 
trarily a limiting interaction distance beyond 
which the energy shall be neglected and a pair 
not counted as a cluster. Regions of less than 
normal density are not accounted for at all. 
Above T,, this objection must always be valid, 
there being no actual demarcation between 
clustered and unclustered phases. But below 7, 
the situation is surely improved by the existence 
of a definite distinction between the two phases. 
For clusters sufficiently large, the intrinsic pres- 
sure will operate and render unstable clusters of 
density intermediate between liquid and vapor. 
The limiting interaction distance thus ceases to 
be arbitrary. However, for small clusters this 
definite distinction between phases can still not 
be claimed ; in fact, our present results on W: for 
the saturated vapor, shown in Fig. 1, definitely 
indicate that the clusters formed in the saturated 
vapor have densities which increase towards T,, 
and which must, therefore, be intermediate be- 
tween liquid and vapor densities. In spite of these 
failings, the method has proved very convenient, 
and productive of interesting results, which do 
not appear to be without significance. 


4J. E. Mayer and S. F. Streeter, J. Chem. Phys. 7, 
1019-1025 (1939). : 
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Raman Spectra of Some Ethers Containing One or More Phenyl Groups 
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Raman frequencies, estimated intensities and depolarization factors are listed for n-propyl- 
phenyl, n-butylphenyl, ethylbenzyl, n-butylbenzyl and dibenzyl ethers. Intense frequencies 
characteristic of monosubstituted benzenes whose values were constant in the five ethers are 
618(5)0.8, 1001(10)0.1, 1029(7)0.3, 1156(4)0.8, 1176(4)0.8, 1588 and 1604 (6)0.7, 3058(9)0.3. 
Frequencies characteristic of the C-O—C bond may be those near 440, 900 and 1125 cm™. 


INTRODUCTION 


N a recent study of the simultaneous catalytic 

dehydration and dehydrogenation of benzyl 
alcohol,! it was desired to obtain information in 
regard to the composition of an intermediate 
fraction by use of Raman spectra. In this 
connection it was found necessary to determine 
the Raman spectrum of dibenzyl ether since no 
data for this compound were found in the 
literature. Since that time, depolarization factors 
for this ether and displacements, estimated 
intensities and depolarization factors for ethyl- 
benzyl, n-butylbenzyl, n-propylphenyl and n- 
butylphenyl ethers have been obtained. The five 
ethers are related in that each contains at least 
one phenyl group and may thus be regarded as 
a monosubstituted benzene. 


EXPERIMENTAL 


The experimental technique and the apparatus 
employed were essentially the same as described 
in a previous paper.? The spectra were obtained 
in the liquid state, the liquids being carefully 
purified by distillation in a 30-cm, continuous 
spiral packed, heated column. The slit width 
was 0.08 mm for the regular exposure and 0.12 
mm for the polarization exposures. Previous 
data apparently have been obtained only for 
dibenzyl ether.! Details regarding the samples 
(all Eastman products) are as follows: 


1 Douglas Meigs, Thesis, Armour Institute of Tech- 
nology, 1940. 

*Forrest F. Cleveland, M. J. Murray, Herschel H. 
(ioe and Julia Shackelford, J. Chem. Phys. 8, 153 


n-propylphenyl ether. b.p. 77° at 17 mm; exposure 
1.0 hr.; polarization exposures 2.5 hr. each 
n-butylphenyl ether. b.p. 95° at 17 mm; exposure 
1.3 hr.; polarization exposures 2.0 hr. each 
ethylbenzyl ether. b.p. 70° at 15 mm; exposure 
3.0 hr.; polarization exposures 2.0 hr. each 
n-butylbenzyl ether. b.p. 88° at 6-7 mm ; exposure 
1.8 hr.; polarization exposures, 2.8 hr. each 
dibenzyl ether. b.p. 176° at 20 mm; exposure 
1.0 hr.; polarization exposures, 1.0 hr. each. 


RESULTS 


The Raman displacements Ay in cm, esti- 
mated intensities J and depolarization factors p 
are given in Table I. Certain lines that are 
unusually sharp or broad are designated by the 
letters s and b, respectively. Data in regard to 
which there is some uncertainty are enclosed by 
parentheses. The frequencies joined with a brace 
represent cases in which the corresponding lines 
were unresolved on the polarization spectrogram 
and the depolarization factor given in such cases 
is the value read for the unresolved line. In a 
few instances it was not possible to list a numeri- 
cal value for the depolarization factor but only 
to designate it as polarized (P). 


DIscussION 
Phenyl group frequencies 


The molecules concerned in the present study 
are of interest not only as ethers but also as 


monosubstituted benzenes.* Kahovec and Reitz‘ 


3 For a review of the Raman spectra of benzene and its 
derivitives, see K. W. F. Kohlrausch, Physik. Zeits. 37, 
58-79 (1936). 

4L. Kahovec and A. W. Reitz, Akad. Wiss. Wien, Ber. 
145, IIb, 1033 (1936). 
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TABLE I. Raman spectra of ethers containing one or more phenyl groups. 


n-PROPYLPHENYL n-BUTYLPHENYL ETHYLBENZYL n-BUTYLBENZYL DIBENZYL — 
by I p | I p | I p | I | I p 
185 1 —| 171 1 0.9| 190 1 _ 
248 4 0.8 245 4 0.7 | 227 3b 0.9| 225 2 — | 200 3 0.9) 
319 0.6 | 291 3 0.8) 323 4 0.7} 285 2 — | 226 2b 0.8 | wio 
353 1 
413 4 0.7 
427 2 0.7 
450 3 454 1 
512 1 _ 510 1 0.9 
600 1 
614 5s 0.7 615 5s 0.7} 619 6s 0.8} 620 5s 0.9} 620 5s 0.7 | we 
753 4 0.7 | 754 3 0.9| 743 2 — | 749 1 740 2b 0.8 | y(CH)? 
775 1 — 
817 5 797° 3 - 803 4 0.4| 804 3 801 4 ws 
. O4 | 812 3 0.5 0.7 
832 3 827 én 822 5 0.4| 819 4 816 4 (CH) 
843 849 3 | 832 844 0.6 
878 2 <0.8 878 1 _— 
886 1 0.6 | 908 1 <0.8| 887 3b P | 902 1 — | 903 3 0.6 
947 2 0.8 
975 4 <0.7 974 1 _ 
999 10 0.1 999 10 0.1} 1003 10 0.1{| 1002 10 0.1} 1001 10 0.1 | w* 
1029 6 0.2 | 1029 8 0.2 | 1030 a 0.3 | 1029 7 0.3 | 1028 7 0.4] wn 
1068 1 — | 1068 3 —_ 
1098 — 


have listed eighteen frequencies which ‘more 
or less unambiguously” correspond to char- 
acteristic vibrations of monosubstituted benzenes 
(CeHs—X). The frequencies in Table I which 
correspond to these are set in boldface type and 
the designations of the frequencies as given by 
Kahovec and Reitz are written in the column 
at the right. While some of the selections are 
rather arbitrary, the frequencies agree, on the 


whole, quite well with those listed by Kahovec 
and Reitz. Depolarization factors for some of 
these phenyl group lines, are listed by Kohl- 
rausch* for X = NHe, CH3, Cl and Br. The mean 
of the numerical values of the depolarization 
factors listed by Kohlrausch for the various 
frequencies of the four compounds and the mean 
of the corresponding values obtained here for 
the ethers are compared in Table II. Excluding 
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RAMAN SPECTRA OF PHENYL ETHERS 


TABLE I.—Continued. 


n-PROPYLPHENYL n-BUTYLPHENYL ETHYLBENZYL n-BUTYLBENZYL DIBENZYL 
Av av p Av I p av I p ay I p 
11311 ie 0.6|1124 4 0.6|1125 1 
1155 3 0.9 | 1152 4 0.8|1159 4 0.8\1155 4 0.8|1158 4 0.9 | 6(CH) 
11744 0.6 |1173 0.81179 4 0.8|1176 4 0.8|1180 4 0.9 | 6(CH) 
1227 1 0.9 |1222 1 0.7 
1245 0.4 |1246 6 04/1210 7 0.2}1210 7 04/1210 6b 0.3| ws 
1281 2 — 11298 5b 0611273 1 0.8|1296 46 0.8 
1299 4 0.6 
1340} — 0.7 13611 
1393} — |1390 1 0.9 1385 1 _ 
14382 ~ 2) 14364 3(CH) 
1456 4 0.7 | 1450 1452 5b | 1454 1455 4 
1467 1 0.7 
1500 — |1497 —|1485 1 0.7|1481 4 0.9|1476 3 ws? 
1587 4 1590 4 1588 3 1588 3 4 08 
2872S 0.2 |2870 7 0.1|2847 4 2863 8b 2856 
2867 
29124 0.7 |2911 6 0.6 2911 0.6 
2937 «6 0.3 | 2937 5 0.2|2931 8 0.1| 2937 8 0.3|2932 3b 04 
2973.2 0.9 |2963 2 0.8|2980 4 0.7|2960 3 0.9|2976 1 sna 
3015} 3008 1 —|3008 1 —|3004 1 — | »(CH) 
3040 4 —|3045 4 3040 3 
3066 9 0.3 |3068 8 0.3|3062 9 0.3 | 3064 3061 
3165} —|3166 } —|3162 — | »(CH) 
3206 36 —|3202 36 —|32004 — 


the uncertain values for w;, w; and y(CH), the 
agreement is quite good, the average deviation 
being 0.08. 

The phenyl group frequencies corresponding 
to seven of the more intense lines are remarkably 
constant in a wide variety of monosubstituted 
benzenes. The mean values of these constant 
frequencies obtained from the list of Kahovec 
and Reitz (for X=NH2, OH, F, CHs, Cl, SH, 
Br and I) are compared with the corresponding 


mean values for the present ethers in Table III. 
Mean intensities and (for the ethers) mean 
depolarization factors are included. 


Frequencies near 3000 


The investigations of Ingold and associates® 
and of Langseth and Lord® indicate that the 


5 Angus, Bailey, Hale, Ingold, Leckie, Raisin, Thompson 
and Wilson, J. Chem. Soc. 971 (1936). 

6 A. Langseth and R. C. Lord, Jr., Kgl. Danske Vid. 
Sels. Math.-fys. Medd. 16, 1-85 (1938). 
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frequency 3062 should be assigned to the C—H 
valence frequency of the phenyl group and that 
the less intense frequencies at 3009 and 3164 are 
overtones. The 3062 line is resolved into two 
components in the benzyl ethers. 

The weak line near 3204 may be the first 
overtone of the 1600 frequency. It was observed 
only for the benzyl ethers. A strong line at 2911 
appeared for the three ethers containing propyl 


or butyl groups but was not observed for the 


other two. This frequency does not in general 
appear for methyl compounds and is weak in 
ethyl compounds, but is strong in all normal 
alkyl hydrocarbons containing chains of three 
or more carbon atoms. It thus appears to be 
more characteristic of the chain than of the 
CH; group, to which Fox and Martin’ have 
assigned the weak infra-red band at 2912. It 
may be that the frequency sometimes observed 
in the ethyl compounds is an overtone of the 
1450 frequency. 

The polarized frequency near 2860 (») and 


the depolarized frequency near 2945 (v,) have. 


been attributed to the —CH2— group.”*® The 
polarized frequency is higher in the two phenyl 
ethers than in the three benzyl ethers. In 
ethylbenzyl ether it is resolved into two compo- 
nents. As would be expected, it is strongest in 
the butyl ethers which contain the greatest 
number of —CHz— groups. The depolarized 
frequency has a value considerably higher than 
the value 2945 assigned to the ‘‘normal methylene 


TABLE II. Mean depolarization factors of phenyl group lines. 


Ko#L- Kout- 

FREQUENCY | RAUSCH | ETHERS || FREQUENCY | RAUSCH | ETHERS 
¥ 0.86 | 0.83 3 0.07 | 0.10 
wi0 0.78 | 0.73 @1 0.18 | 0.28 
w1 0.29 | (0.65) |}1155, 6(CH)} 0.88 | 0.84 
we 0.86 | 0.76 ws 0.21 0.34 
w4 0.12 | (0.50) 7, 8,9 0.85 | 0.68 
+(CH) dp | (0.50) 


asd diy? Fox and A. E. Martin, Proc. Roy. Soc. A175, 208 
8 (a) K. W. F. Kohlrausch, Der Smekal-Raman Effekt, 
p. 190; (b) Forrest F. Cleveland and M. J. Murray, J. Am. 

Chain. Soc. 62, 3185 (1940). 
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TABLE III. Constant frequencies of the phenyl group. 


FREQUENCY KAHOVEC AND REITZ ETHERS 
Ap I Ap I p 

we 615 4 618 5 0.8 
3 1001 11 1001 10 0.1 
11 1024 6 1029 7 0.3 
6(CH) 1157 4 1156 4 0.8 
7, 8,9 

v(CH) 3062 10 3058 9 0.3 


group.”’8* It is highest in the ether containing 
the ethyl group, lower in propylphenyl and 
dibenzyl ethers, and lowest in the two ethers 
containing the butyl group. This slight variation 
in frequency may be caused by the presence of 
different atoms or groups adjacent to the 
—CH.— group, to a neighboring unresolved 
methyl group frequency, or to shifts involving 
some resonance interaction. Fox and Martin 
have concluded from infra-red studies® that the 
symmetrical and unsymmetrical frequencies of 
the —CH.— group are 2853 and 2926, respec- 
tively. The latter is more intense than the former 
in the infra-red spectra. 


Frequencies of the ether group, R-O—R’ 


Pai,® using 1102, 921 and 416 as fundamentals 
for dimethyl ether, calculated a value of 102° 
for the C—O—C valence angle. This is in the 
neighborhood of the value 111+4° found by 
Sutton and Brockway” from an electron diffrac- 
tion investigation. Recently Crawford and 
Joyce" have assigned the frequencies 1122 and 
940 in dimethyl ether to vibrations involving a 
stretching of the C—O bond and the frequency 
412 to a C—O-—C deformation vibration. 
Frequencies observed for the present ethers 
which may correspond to these are the ones 
near 900, 440 (observed in only three cases) 
and 1125. 

® Gopala Pai, Current Sci. 2, 386 (1934). 


10 L. E. Sutton and L. O. Brockway, J. Am. Chem. Soc. 


57, 473 (1935). 
1B. L. Crawford, Jr. and L. Joyce, J. Chem. Phys. 7, 


307 (1939). 
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Line absorption spectra of thirteen sulfate and selenate 
chrome alums have been measured at 15—20°K. Pronounced 
differences in the spectra are attributed to structural 
differences in the crystals and are related to the size of 
the anions and the univalent cations. Anomalous behavior 
on change of temperature is shown by potassium sulfate 
and ammonium sulfate and selenate alums. Transitions 
disrupting the lattice and significantly altering the spec- 
trum occur in ammonium chromium sulfate at 81+1°K 
and in ammonium chromium selenate at 106+5°. Combi- 
nation of electronic frequencies of Cr**+* with internal 


oscillations of anions are weak. Combination of lattice 
vibrations with both upper and lower electronic states of 
Cr+** are found. Bands are observed in all the alums 
shifted about 330, 570 and 800 cm™ to higher frequencies 
than the electronic absorption lines. These are thought 
due to photo-excitation of a nonhomopolar complex, 
probably Cr(H2O).+*+*. Repetition of the electronic 
multiplet and associated vibrational structure at some of 
these higher frequencies is found in several of the alum 


spectra. 


INTRODUCTION 


ISCONTINUOUS absorption spectra of a 

large number of coordination compounds 
of trivalent chromium at liquid-air temperature 
have been recorded.’-* Sauer, especially, has 
studied the effect of isomorphous replacement 
in a group of eleven sulfate and selenate chrome 
alums at —190°C, and in a smaller number of 
these alums at —78° and +418° also. The 
discrete absorption is in general found within a 
range of a few hundred wave numbers in the red 
and orange and merges into a strong continuum 
with maxima at approximately 5700 and 4000A. 
With lowering temperature the continuum 
narrows, revealing additional lines and bands. 
The latter also become narrower and less diffuse 
and often resolve into finer structures. 

An investigation in this laboratory’ of potas- 
sium chromium sulfate alum was carried out in 
the temperature range —190° to —259°. The 
qualitative effects of changing temperature 
measured by Sauer are found in heightened 
degree in the lower range. The bands of con- 
tinuous absorption are so narrowed at the 
temperature of liquid hydrogen that a series of 
sharp lines is revealed in the blue which seems 
to show a relationship to the violet continuum 
similar to that of the red-orange discrete spec- 
trum to the yellow-green continuum. In addition 


1 Schnetzler, Ann. d. Physik 10, 373 (1931). 

? Snow and Rawlins, Proc. Camb. Phil. Soc. 28, £22 (1932). 
* Sauer, Ann. d. Physik 87, 197 (1928). 

* Spedding and Nutting, J. Chem. Phys. 2, 421 (1934). 


a gradual change in the relative intensity and 
structure of the strong ‘‘characteristic doublet”’ 
and its associated weaker lines near 6700A is 
noted which suggests the presence of a multiple 
lower state. Magnetic behavior of potassium 
chrome alum at temperatures down to about 1°K 
makes it seem unlikely that the lower electronic 
state of Crt+*+ is multiple. Alternatively, 
characteristic oscillation frequencies of the 
crystal lattice or of some ion or complex com- 
posing it might conceivably combine with the 
single electronic state of Crt+*+ to give an 
ensemble of states showing the observed temper- 
ature behavior. 

The spectra of other alums, notably the 
selenates, are richer in lines than potassium 
chromium sulfate, and the lines are superior in 
sharpness at the temperature of liquid air. With 
the increased definition and precision of meas- 
urement made possible at liquid hydrogen 
temperature it was hoped that constant fre- 
quency differences, repetition of line groups, or 
frequencies equal to infra-red or Raman fre- 
quencies might be found in the absorption 
spectra which would give a clue to the origin 
of the numerous absorption lines of these alums. 
In addition, further information on the effect of 
temperature on the structure of the characteristic 
doublets was desired. Accordingly, an investiga- 
tion was made of the visible absorption spectra 
of as many different chrome alums as could be 
prepared. The spectra of single crystals of the 
following thirteen chrome alums have been 
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TABLE I. Absorption lines of the chrome alums. Intensities with respect to the background (1) have been estimated roughly 
on the scale 1-10. The letter s indicates a line or band whose edges are sharp and well-defined; d a line or band whose edges 
are diffuse and less precisely measurable. The total width (w) of most of the lines has been measured and is reported in cm. 


I w I w I w I w 
(a) () 
20° 20° 
14,931 1s 15,178.5 2d 11 14,898 1s 15,240 4s 
937.5 10s 1.5 205 1d 10 906 1s 246 2d 
941.5 10s 1.1 255 3d 12 922.1 10s 1.7 252 3d 
944.0 3s 270 3d 12 927.9 10s 1.3 259 2d 
946.1 3s 292.5 1d 13 980 1s 17 279 1d 24 
978.5 1d 11 373 1d 17 15,005 1s 353 1d 27 
15,004 1d 12 431 1d 14 013 2s 420 1d 20 
028 1d 8 465 1d 10 020 1s 454 2d 12 
048.3 4s 1.6 501.8 8s 5 031 4s 486 8s 8 
052.2 4s 1.2 510.5 3d 5 037 4s 497 5s 7 
057 1s 549 2d 4 042 1s 533 2d 8 
065 1d 8 553 2d 4 049 1s 5 555 1d 
079 2d 6 574 1d 8 061 3s 638 1d 25 
087 2s 675 1d 9 068 5s 660 1d 19 
093.5 3s 1 682.5 2d 7 074 3s 697 4d 22 
098 2s 709 5d 20 080 4s 721 1d 9 
103 1s 731 4d 11 086 4s 741 4d 14 
106.5 4s 755 5d 13 092 4s 760 2d 8 
111 2s 781 1d 8 098 4s 777 1d 14 
113 1s 800 1d 5 104 2d 822 1d 28 
118 1d 835 id 9 113 4s 863 1d 21 
123 1d 849 2d 11 121 3s 910 1d 23 
128 1d 4 875 3d 17 128 3s 959 1d 42 
134.5 1s 1 923 2d 12 135 2s 16,052 1d 78 
142 1s 979 2d 19 164 4d 14 1d 16 
146 1s 16,042 . 1d 34 180 1s 279 1d 34 
162.5 1d 5 186 1s 
194 1s 
(b) RbCr(SeO.)2+12H20 233 2s 


20° 
(d) KCr(SeOx)2-12H20 


14,923.1 10s 1.9 | 15,134.7 3s 3 15° 
28.0 10s 1.4 147 1d 8 
930.4 3s 163 3d 13 
942 1d 187 2s 2 14,901.7 10s 2.2 | 15,092 4s 
952 1d 192 2s 3 906.6 10s 1.9 096 1s 
990 1d 228 1d 7 908.8 3s 105 1d 6 
995 1d 234.6 2s 3 926 1s 112 1d 
15,003 1d 239.8 2s 4 934 1s 122 1d 6 
id 260 3d 16 940 1s 135 1d 7 
011 1d 282 1d 13 953 1s 140 2d 3 
017 1d 359 1d 24 959 1s 163 1d 
022 id 364 2s 964 1s 168 1d 
033.2 4s 2.2 419 2d 9 971 1s 203 1d 
038.3 4s 2.3 455 3d 9 976 1s 205 1d 
044 1s 489.5 7s 7 981 1s 213 3s 
052 1d 8 497.2 5s 7 986 1s 218 4d 
064.7 3s 3 534.9 3s 4 994 1d 4 237 4d 10 
070.9 4s 2 540.7 3s 3 15,004 1d 260 1d 15 
076.0 3s 3 555 2d 5 012 2d 335 id 17 
081.7 3s 3 646 2d 27 018 3s 394 1d 10 
087.3 3s 3 672 1d 6 043 2d 433 1d 12 
092 1s 697 3d 22 050 3s 469 4d 10 
094.5 3s 2 719 2d 8 054 1s 479 4d 10 
100.2 3s 2 742 3d 14 061 1s 516 2d 
108.2 3s 2 764 id 7 068 3s 521 2d 
114.0 5s 3 783 1d 7 073 2d 532 2d 
119.5 4s 2 823 1d 20 077 1s 537 2d 
124.2 2s 3 862 1d 30 083 5s 623 1d 19 
128.2 3s 3 087 3s i 675 1d 18 
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SPECTRA OF CHROME ALUMS 
TABLE I.—Continued. 
I w I w I w cm! I w 
(e) (i) NHsCr(SO«)2- 12H20—Continued 
Clear crystal Shattered crystal 
14,851 id 20 15,152 2d 35 | 14,931.8 10s 1.1 | 15,276 5d 15 
903.5 10s 6 227 3d Sts 935.9 10s 2.2 366 id 15 
913.8 10s 6 327 ld 44 941.6 10s 1.4 401 1d 11 
955 id 26 392 id 25 944.5 10s 1.2 422 4d 10 
15,008 id 30 434 2d 17 947.5 4s 2 454 4s 
056 id 33 472 4d —- 30 955.5 2s 3 465 3d 7 
097 3d Ss 332 530 2d 33 | 15,007.5 3s 7 479 3d 7 
025 2s 493 3d 4 
042. 3s 13 
(f) *12H20 066 1s 531 1d 10 
Shattered crystal 082 2s 621 1d 20 
114 2s 1 
14,846.0 9s 2.2 | 15,164 1d 14 122 = 794 1d 13 
878.9 10s 2.6 8 1d 7 130 5s 8 7163 1d 39 
903 5d 6 202 1d 10 140 3. 817 1d 7 
917.6 5s 1.4 245° 3d 9 180 id 2 
921.5 10s 2.3 255 1d 11 
925.0 10s 1.8 277 id 10 
979 1s 378 1d 10 (j) 
987 1s 401 3d 10 
15,015 432 5s 7 
0 2 
093 id 6 > re Bs 14,916 2s 15,239 4s 6 
103 4s 3 590 id 12 923 2s 248 2d 
120 id 1 609 id 10 925 2s 255 5d 
128 id 6 648 1d 10 = 1.3 
146.5 1d 5 692 id 22 = 
943 2s 376 1d 
(g) 949 id 385 2d 
957 1d 411 1d 
1 4 
14,917 5d 9 15,399 ids 63 937 id 454 3d 
927 10s 3 455 id 24 904 id 462 d 
935 10s 4 484 4d 12 15,000 1d 488.5 8s 3 
940 10s 5 501 4d 10 009 1s 492 8s 4 
15,124 2d «=a 532 id 15 021 id 13 500.5 6d 3 
185 id 48 568 id 20 033 1s 503 6d 2 
259 4d 37 038 1s 520 1s 
045 1s 6 529 2s 
‘ 053 1s 41 1s 
Shattered crystal 063 2s 5 648 1d 8 
067 2s 4 667 1d 1 1 
14,8662 9s 24 | 15,134 2d 24 = 
9023 10s 23] ‘187 id 37 08s 
9287 9 15] 233 id 15 
101 3s 734 1s 5 
932.3 10s 2.1 265 3d 
938.4 9s 3.2 422 13 107 3s 
941.6 9s 2.9 462 477 30 115 2 ‘ in - 8 
15,069 ld 28 493 2 4 
1 1d 20 512 1d 7 9 
142 3s 843 
(i) 145 882 1d 18 
152 2 926 1d 29 
annual 159 2d 6 977 ld 36 
170 3d 9 
14,868.0 10s 1.1 | 15,196 3d 19 177.5 3d 7 
885 1d 5 212 1d 187 3d 
904.9 10s 1.2 235 2d 15 197.5 3d 5 
919 2s 7 250 1d 205 3d 6 
928 Is 4 264 5d 9 231 1d 
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TABLE I.—Continued. 


I w I w w I w 
(k) TICr(SO.)2:12H:0 (m) RbCr(SO,)2-12H20 
85° 20° 
14,745.5 1d 7 15,135 7s 7 14,814 1s 15,158.5 7d 5 
796 1d 4 147 5d 6 913 1s 164 1s 
822.5 1d 7 153 4d 5 922 1s 177 1d 6 
876 1d 8 183 5d 8 930 1s 191 5d 11 
891 1d 4 209 4d 11 942 1s 200 2d 
913 3s 4 274 8d 20 945.5 10s 1.4 214 4d 9 
920.5 1s 3 390 4d 8 947.4 10s 1.2 262 1d 
938.5 10s 1.8 432 4d 7 950.7 2s 275 8d 8 
941.6 10s 2.0 466 4d 14 954.0 3s 281 4d 4 
948 2s 499 Od 4 970 3s 3 288 4d 4 
953 2s 505 4d 7 997.5 1d 7 398 1d 10 
961 2s 4 533 4d 8 15,012 3s 4 433 2d 
969.5 2s 5 575 3d 21 020. 1d 9 443 2d 
987 2s 7 639 1d 11 031 1d 8 467.5 1d 7 
15,001 3d | 672 1d 13 039 1d 8 476 3d 6 
011 2d 5 709 6d 19 048.5 2d 7 507.5 8d 7 
021 1d 2 734 1d 7 055.5 1d 7 513 1d 4 
036 1d 10 757. 7d 11 062.5 6d 5 541 4d 9 
058 5s 8 804 1d 20 067.5 6d 5 553 1d 
069 2d 2 838 3d 8 078.5 3d 7 562 1d 
083 5d 6 869 3d Pat 089.9 7s 4.2 581 1d 15 
089 4d 6 921 1d 12 096.3 7s 2.2 627 1d 23 
102 4s 3 981 1d 15 102 4s 655 1d zi 
114 5s 4 107.1 7s 2.6 686 1d 23 
114 3s 716 3d 16 
421.7 7s 241 737 3d 6 
(1) TICr(SO,)2-12H:0 129 3d 2 764 4d 10 
“i 136 3d 3 781 id 14 
142.5 9d 5 846 2d 16 
14,911 1s 15,259 1d 2 151.5 7d 7 879 2d 22 
934 1s 285 8d 6 
939 4s 394 6d 10 14,875 1s 15,109 4s 
— 942.2 10s 1.9 433 5d 7 884 1s 117 1s 
944.8 10s 2.0 439 6d 5 889 1s 124 1s 
950 1s 2 469 5d 6 895.4 10s 1.4 134 2s 8 
956 1s 475 6d 5 907 1s 148 1s 
962.4 As 503 9s 4 918.9 10s 1.5 157 1s 
974 3s 4 508 6d 5 925 1s 168.5 2s 7 
983.5 1d 5 537 7d 12 929 1s 192 1s 8 
994 3s 5 556 2d 8 939 1s 216 1d 
15,009.1 4s 4 583 3d 19 944 1s 223 1d 
015.4 3s 4.2 624 1d 29 949 1s 229 1d 
025 2d 4 649 1d 18 956 1s 240.5 4d 5 
042 3d 8 679 1d 31 963 1s 264.5 2d 7 
050 1d 4 713 6d 19 970 1s 347 2d 9 
060.5 7s 3.9 733 5d 9 979 1s 382 3s 6 
064.4 7s 3.9 761 7d 7 987 1s 402 2s 11 
072.5 5s 5.6 776 1d 10 999 1s 419.5 3s 5 
085.6 7s 4.0 804 1d 13 15,007.5 1s 3 437 3s 6 
093.6 5s 2.4 816 1d 11 010 1s 2 449 2d 4 
099 1s 841 3d 21 017.5 1s 5 460 5s 4 
105.2 8s 3.1 872 3d 20 027 2s 6 471 4s 
111 1s 4 922 1d 18 034 2s 486.5 3s 7 
117.7 8s 4.2 937 1d 047 2s 510.5 2s 9 
124.5 2d 3 981 1d 22 052 2s 530 1d 
132 1d 2 16,030 1d 19 058 2s 4 553 1d 12 
139 8d 6 063 1d 28 061.5 2s 3 617 1d 15 
150 6d 8 104 1d 24 070 3s 643 1d 21 
156.5 6d 5 164 1d 26 077 2s 669 2d 20 
176 1d 202 1d 14 083 1d 714 2d 42 
187 7d 11 247 1d 30 088 2d 757 1d 12 
198 3s 4 290 1d 25 093 1d 819 1d 22 
212 6d 9 098 1d 850 1d 13 
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TABLE I.—Continued. 


cm~! I w cm"! I w 


I w cm"! I w 


(0) *12H20 


(q) CsCr(SO4)2-12H20 
20° 


14,573.9 10s 1.6 14,946 1d 30 
587 1s 15,000 1d 30 
621.5 1d 9 047 1d 16 
629 3s 084 2d 20 
653.5 1d 9 127 4s 7 
684.5 id 7 145 1s 8 
692.5 4s 5 155 1s 
701 1s 163.5 2s 3 
730.5 2s 7 176 1s 9 
738 3s 187 3s 4 
741 1s 214 2s 7 
756 1s 8 240 1d 16 
765 1s 7 269 1d 21 
789 1d 292 1d 3 
797.5 4s 3 302 1d 4 
804.5 1s 3 321 1d 13 
816 1d 11 351 1d 17 
832 1d 7 459 1d 6 
876 1d 11 508 1d 48 
898 1d 15 583 1d 20 
910 1d 


14,501.2 10s 1.5 | 14,830 1d 11 
507 1s 946 id 8 
517 1s 955 1d 9 
549.5 1s 5 974 1d 9 
557 1s 6 15,000 id 
566.5 1d 5 011 3d 7 
580 1d 11 042 1d 
593 1d 12 046 id 
617 1s 053 2s 
625.5 3s 7 060 id 
634 Is 072 1d 
646.5 2s 5 083 2s 
660 2d 9 099 1d 11 
669 2s 132 1d 29 
677 is 241 id 12 
689.5 1d 5 268 id 18 
719.5 1d 5 393 id 11 
730.5 3s 7 416 id 7 
745 1d 4 437 1d 24 
752 1d 513 1d 12 


14,479.3 10s 2.2 14,776 1d 
517 1s 787 1d 
523.5 5s 2.1 810 2d 
528 1s 817 3s 
535 1s 827 1d 10 
537 1s 841 1d 
540.5 2s 5 861 1d 10 
549 1s 927 3s 8 
553 3s 2 947 1d 5 
555.5 3s 3 960 3s 8 
565 1s 975 1d 8 
571.7 4s 995.1 9s 4 
575.1 4s 15,034.2 9s 4 
578 4s 040.1 10s 2.4 
582.5 4s 052.5 8s 4.6 
586.0 4s 059 4s 
600.7 9s 17 065 5s 6 
612.3 5s 3 086 1d 8 
618 1s 096 3d 
626.5 4s 5 112 1d 23 
632 5s 136 1d 
635.5 4s 5 147 1d 
643 1s 154 5d 8 
650 1s 161 4d 6 
655.5 3d 7 176 2d 9 
663 1d 184 2d 8 
668 1d 212 3d 27 
675 2d 260 1d 12 
679 2d 276 1d 
683 2s 298 1d 11 
693 4d 332 2d 25 
696 4d 371 2d 7 
699 4d 417 3d 30 
708.1 8s 3.4 469 1s 
713 2s 514 1d 12 
719 4s 544 1d 
727 1d 579 1d 7 
736 1d 622 1d 36 
742 1d 
755 2s 6 


photographed at 85° and 20° or 15°K: both the 
sulfate and selenate alums of potassium, rubid- 
ium, cesium and thallium, ammonia and methyl- 
amine, and the sulfate alum of hydroxylamine. 


EXPERIMENTAL PART 


Large single crystals of most of the alums 
were obtained by evaporation at 6°C of solutions 
containing equivalent amounts of the sulfates or 
selenates of chromium and the various univalent 
cations. Rubidium and cesium chrome alums 


were obtained as large crystals only after adding 
concentrated nitric acid to increase the otherwise 
slight solubility of these salts. The crystals gave 
no test for nitrate and their spectra were identical 
with the spectra of the small crystals separating 
from solutions in water alone. 

A guanidine sulfate alum has been reported.® 
The spectrum of the guanidine chromium sulfate 
which we prepared was widely different from 


that of any of the chrome alums. Crystallo- 


5 Canneri, Gazz. Chim. Ital. 55, 611 (1925). 
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graphic examination of the guanidine compound 
showed that it has hexagonal symmetry and 
thus that it is not an alum. 

Preparation of alums of many other substi- 
tuted amines was attempted unsuccessfully. 
Among them were both the sulfate and selenate 
alums of ethylamine, dimethylamine, propyl- 
amine, and isopropylamine; and sulfate alums of 
triethylamine, tetramethylamine, aniline, methyl 
guanidine and dimethyl guanidine. The univalent 
cations of these organic substances are evidently 
too large to be included in an alum lattice. The 
few cubic crystals sometimes obtained were 
attributable to ammonia or methylamine im- 
purities in the amines employed.® 

Hydrazine sulfate chrome alum and sodium 
sulfate and selenate chrome alums’~* have been 
described. We were able to prepare neither these 
alums nor silver sulfate nor hydroxylamine 
selenate chrome alum. 

The spectra were photographed at the temper- 
ature of liquid air and liquid hydrogen, em- 
ploying apparatus and technique described 
previously. The thickness of the single crystals 
used was from 3 to 7 mm. No investigation was 
made of the polarization of the absorption lines 
or of the directional dependence of their intensity. 


DESCRIPTION OF THE SPECTRA 


The position of the centers of the absorption 
lines, the width and character and the relative 


THE RED ABSORPTION SPECTRA OF THE CHROME ALUMS MCr (ROb-i210 


uM 

a & A L man 

- ».. moe. 

om wi led 

om 

aos 


Fig. 1. Graphical representation of the alum spectra. 
The solid circles at the extreme left indicate that the 
corresponding crystals were shattered. 


6 Wendekamm, Zeits. f. Krist. 85, 35 (1933). 

7™Sommer and Weise, Zeits. f. anorg. allgem. Chemie 94, 
51°(1916). 

8 Orlov," Ukrainskii Khem. Zhurnal 3, 115 (1928). 
ass Bhattacharya, J. Ind. Chem. Soc. 7, 765 
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intensity are given in Table I for all the alums 
at hydrogen temperatures, with the exception of 
KCr(SO,)2:12H,O, whose spectrum has been 
reported in detail from this Laboratory. Data at 
higher temperatures on TICr(SO,)2°12H2O are 
included because they indicate several excited 
levels in the lower state not found by Sauer, and 
data on the ammonium alums because of unusual 
temperature behavior. The positions reported 
for the strong, sharp lines in the region of the 
characteristic doublet are probably correct within 
a fraction of a cm while for the faint, diffuse, 
multiple lines and bands the error may well be 
one or more cm~!. The spectra are represented 
graphically in Fig. 1 and by photographs in 
Figs. 2, 3 and 4. 

At the long wave-length end of the spectrum 
of all the alums is a line or group of lines of 
singular sharpness and intensity, originating 
probably in pure electronic jumps within the 
3d shell of the Cr+++ ion. Extending to perhaps 
700 to 1300 cm on the violet side of the strong 
group is a series of weaker lines whose sharpness 
and resolution decrease very greatly with their 
separation from the strong group. The discrete 
absorption ultimately is masked beneath a 
continuous absorption band. The most striking 
changes on isomorphous replacement of univalent 
cations and of anions appear in the electronic 
spectrum. Three classes may be distinguished. 
The first class, which comprises most of the 
alums, has a strong doublet accompanied by a 
variable number of much fainter lines or satellites 
grouped close to 6700A. The spectrum of alums 
of the second class is shifted in its entirety about 
400 cm-! toward the red and in place of a 
characteristic electronic doublet there is a single 
strong, sharp line near 6900A. The third class 
shows about six strong lines near 6700A. 

To the first class belong the four alums of 
rubidium and thallium, the selenate alums of 
cesium and potassium, the sulfate alum of 
hydroxylamine and the form of all the ammonium 
alums that is stable at higher temperatures. 
The position and the separation of the compo- 
nents of the doublet is variable. On lowering the 
temperature below 85°K the doublets sharpen 
perceptibly, their satellites become much more 
pronounced, and many of the diffuse bands 
narrow and resolve into a number of weak, 
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CsCr(SeOs)2-12H20 


20°K 


CsCr(SeOs)2- 12420 


TICr(SeOs)2-12H20 


TICr(SeOs)2-12H20 


Fic. 2. 


RbCr(SeOs)2- 12120 


KCr(SeQs)2-12H20 20° 


KCr(SeO4)2-12H20 


—| 15,400 
15,100 


sharp, closely spaced lines (about 5 cm~!) around 
15,100 cm“. (See especially 
in Fig. 2.) The selenates of this class show 
remarkable similarity. Their spectra correspond 
practically line for line in respect to their relative 
position and intensity. In the sequence potas- 
sium, thallium, rubidium and cesium there is a 
displacement of the entire spectrum to the violet 
with each succeeding member. The amount of 
the displacement depends largely on the size of 
the cation and not upon its mass as the following 
compilation shows. Av is the shift of the mean 
frequency of the characteristic doublet. Re- 


kt Tit Rbt Cs* 
at wt. 39 204 85 133 
radius 1.33 1.48 1.48 1.69A 
Ap 0 21.0 21.5 35 cm 


A 


placing sulfate by selenate in thallous alums 
induces a red shift of the center of the doublet 
of 19 cm~; in rubidium alums, 20 cm~!. Com- 
parison of sulfate and selenate alums of this 
class produces little evidence of differences in 
their spectra which can be associated with charac- 
teristic oscillation frequencies of these anions. 

The two alums of methylamine and the sulfate 
alum of cesium constitute the second class. The 
spectra are characterized by their diffuseness at 
85°K. The weak absorption bands of the methyl- 
amine alums are, indeed, scarcely visible. At 20° 
they are, however, quite sharp. The attempt was 
made to resolve the isolated strong line of this 
class into a “‘characteristic doublet.”’ The spec- 
trum of the cesium alum was photographed at 
20° in the second order of the grating (dispersion 
1.8A per mm) but no indication of resolution 
was discovered. 
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TICr(SOs)2-12H20 


RbCr(SOs4)2-12H 
Fic. 3. 


KCr(SOs4)2-12H20 


20° 


—| 15,400 
100 


CsCr(SO«)2-12H20 20° 
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Potassium sulfate alum and the low tempera- 
ture form of the two ammonium alums make up 
the third class. Characteristic of this group is 
the large number (six or more) of strong sharp 
lines which appear near 6700A. In their change 
with temperature the ammonium alum spectra 
are as anomalous as the potassium chrome 
sulfate spectrum, but in a different way. The 
potassium alum has notably different spectra 
at 85° and at 14°K, and at intermediate temper- 
atures both spectra appear with prominence 
which varies with the temperature. The whole 
behavior is such as to suggest the presence of a 
close multiplet structure in the lower state. The 
ammonium alums also have “high” and “low” 
temperature spectra, but they do not coexist. 
It is often possible to keep an ammonium sulfate 
alum crystal indefinitely in liquid nitrogen or 
liquid hydrogen in its original perfect, highly 
transparent condition. In this case the spectra 
at the higher and lower temperatures differ only 
in the line breadths. On other occasions an 


| 


apparently perfect crystal retains its trans- 
parency for a few minutes after being placed in 
the cooling liquid and then almost instantane- 
ously becomes practically opaque, as though 
broken into an infinite number of tiny crystals 
by the disruption of the lattice. It is these 
nearly opaque crystals which will be called 
“shattered crystals” in what follows. Exposure 
times for the shattered crystals are of the order 
of one hundred times as long as for the perfect 
crystals and a thickness of 1 mm is ample. The 
spectra are quite unlike (Fig. 4). The high 
temperature form has two strong lines near 
6700A, the low temperature form has six lines 
much more widely spaced. , 
The influence of temperature upon the 
ammonium alums was investigated carefully. In 
liquid oxygen (90°) the sulfate alum gave the 
spectrum of the high temperature form (Table I 
(g)). The crystal was gradually cooled further 
by reducing the pressure on the oxygen until at 
81+1° the crystal suddenly became opaque, and 
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20°K 


*N HiCr(SeOy)2- 121120 


85° 


110° 


| 


*N HaCr(SO4)2- 12H20 


85° 


85° 


20° 


NH;3OHCr(SOs)2- 12H20 


Fic. 4. Spectra of ammonium and hydroxylamine alums. 


its spectrum was now that of the low temperature 
form (Table I (h)). The crystal was warmed to 
room temperature and re-immersed in liquid 
oxygen at one atmosphere. It again gave the 
spectrum of the high temperature form. The 
selenate alum was treated similarly, using 
liquefied natural gas (about 90 percent methane) 
as cooling liquid. In the methane at one atmos- 
phere pressure the crystal gave the high temper- 
ature spectrum (Table I (e)). On lowering the 
temperature the crystal became opaque at 
106+5° as determined roughly by a thermo- 
couple. Its spectrum was now that of the low 
temperature form (Table I (f)). After standing 
in the cooling bath until the temperature had 


A 
* Indicates shattered crystals. 


risen to about 140° the crystal gave the high 
temperature spectrum. These experiments clearly 
indicate the existence of transitions in the crystal 
lattice of the ammonium alums at low temper- 
atures. It has been possible to supercool the 
high temperature form of the sulfate alum to 
15°K, but in no case was it possible to supercool 
the selenate alum. Shattered crystals often regain 
most of their original transparency on warming 
to room temperature. The density of perfect and 
of shattered crystals is practically identical. 
Discrete absorption bands in the blue have 
been found at 20° in the following alums: the 
sulfates and sclenates of rubidium, thallium, 
cesium and methylamine, the sulfate of potas- 
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sium and the high temperature form of the 
sulfate of ammonium. The bands are diffuse and 
very weak. They could be photographed only 
when thick, perfect crystals were available. The 
bands extend into the regions of continuous 
absorption on either side. 


DISCUSSION 


The large difference in the spectra of alums of 
classes I and II is doubtless consequent upon 
small differences in the lattice structure. Lipson 
and Beevers'® have found two distinct structures 
among the aluminum sulfate alums, differing 
only in the position of the sulfate ion along a 
trigonal axis of the unit cell. One position is 
assumed in alums with the largest univalent 
cations, cesium and methylamine, while the 
second position of the sulfate is taken when the 
cations are of medium size, such as potassium, 
ammonium, thallium and rubidium. From their 
absorption spectra one may venture that cesium 
sulfate chrome alum has the one structure, while 
cesium selenate chrome alum has the other. A 
third, more widely different structure has been 
found by Lipson and Beevers for alums with 
small cations. This structure differs from the 
other two in the orientation of the anion along 
the trigonal axis. The spectrum of a chrome alum 
with this structure might well be distinctly 
unlike the others. Attempts were made, accord- 
ingly, to introduce the small ions Na* and Agt 
into a chrome alum lattice, but without success. 

The dielectric measurements of Guillien" are 
of interest in connection with the spectroscopic 
anomalies of potassium sulfate and ammonium 
alums. In sulfate alums of iron, chromium and 
aluminum with ammonium or potassium Guillien 
finds that the dielectric constant decreases with 
temperature down to about —100°, rises when 
the temperature is lowered to —187.5° and then 
very suddenly decreases again. The electric 
absorption and specific heat also have maxima 
near —190°. The phenomena usually show 
hysteresis and the maxima seem to be true A 
points. They are probably to be ascribed to 
cooperative relaxation of librational restraint 


” Lipson, Nature 135, 912 (1935); Proc. Roy. Soc. 
(London) A151, 347 (1935); Lipson and Beevers, Nature 
134, 327 (1934) ; Proc. Roy. Soc. (London), A148, 664 (1934). 

1 Guillien, Comptes rendus 209, 21 (1939). 


upon water molecules or water-containing com- 
plexes with rising temperature. In the ammonium 
chrome alums imposition of this librational 
restraint upon cooling evidently induces strains 
in the lattice which are often relieved in the 
shattering of the crystal. A more stable lattice 
is formed in which the chromium ions are acted 
on by fields of slightly different strength and 
(probably lower) symmetry as is expressed in 
the changed absorption spectrum. The lattice 
of the potassium alum is not similarly disrupted 
and the high and low temperature spectra may 
exist together. 

The nature of the electronic states giving rise 
to the characteristic doublet or corresponding 
lines in other Cr+++ compounds has not yet been 
satisfactorily described. Bowen™ gives the fol- 
lowing positions for terms arising from the 3d* 
configuration of gaseous Crt++, 

4F 3/2, 5/2, 7/2, 9/2 0, 237, 553, 949 cm 

4P 1/2, 3/2, 5/2 14,058, 14,177, 14,476 cm 

2G 7/2, 9/2 15,056, 15,405 cm 

2H 9/2, 11/2 21,067, 21,320 cm 


2D levels were not found but probably lie not far 
from */7. The interval ‘F—°G of the gaseous ion 
is nearly equal to the frequency of the sharp 
lines in the red spectrum of Cr++*+ compounds. 
Similarly, the energy difference between *F and 
277 levels of the free ion is nearly the same as 
the frequency of the blue line absorption spec- 
trum. Hence it has been proposed several times™ 
that the line spectra of the chrome alums do in 
fact represent *F—*G and ‘F?—/J/ transitions. 
However, owing to the drastic modification of 
coupling schemes and selection rules and the 
splitting and shifting of the energy levels of the 
free ion in the crystal fields, the above assign- 
ment of levels is probably not correct.'** 

Precise measurements of the magnetic sus- 
ceptibility have been made only on potassium 
chromium sulfate alum." The strict linearity of 
the susceptibility-temperature curve from. 1.3° 
to 290°K makes it probable that there is no 

2 Bowen, Phys. Rev. 52, 1153 (1937). 

13 Deutschbein, Ann. d. Physik 14, 712 (1932); Joos, 
Ergeb. exakt. naturwiss. 18, 78 (1939). 

13a While this article was in press, papers on the elec- 
tronic energy levels of the chrome alums have been pub- 
lished by Van Vleck (J. Chem. Phys. 8, 787 (1940)) and 
Finkelstein and Van Vleck (ibid. 8, 790 (1940)). 


4 de Haas, Gorter and van den Handel, Proc. Acad. Sci. 
Amsterdam 33, 676 (1930); 36, 158 (1933). 


s 
h 
e 
t 
li 
sl 
as 
ti 
E 
re 
pl 
vi 
eX! 
ple 
wi 
in 
of 
chr 
ion 
ion 
freq 
pou 
met 
trip 
freq 
and 
Rev. 
would 
which 
magn 
The u 
a mul 
kinds 


SPECTRA OF CHROME ALUMS 


excited electronic level in Cr+++ within 200 cm— 
of the ground level.“* Heat capacity data are 
not available to add further evidence on this 
point. The two or more strong, characteristic 
lines of the chrome alums originate then in 
transitions from a single lower state to two or 
more rather closely spaced upper electronic 
states.!5 The weak lines and the diffuse bands 
of higher frequency represent simultaneous 
photo-excitation of the electronic transition and 
oscillation of the lattice or certain of its con- 
stituents. The first (green) continuum may 
have a similar origin and also may include other 
electronic jumps and associated secondary struc- 
ture as yet undetected. The blue absorption 
lines seem to be the counterpart of the weak 
sharp lines around 15,100 cm-! and may be 
associated with a strong, pure electronic transi- 
tion which is hidden in the continuum. 

In absorption spectra of neodymium salts 
Ewald'* has found groups of weak, sharp lines 
repeating strong, sharp, pure electronic multi- 
plets at intervals corresponding to internal 
vibration frequencies of the anions. Photo- 
excitation of the characteristic vibrations takes 
place with moderate probability simultaneously 
with the electronic transitions of the Nd+*+* ion 
in spite of the interposition of neutral molecules 
of an Nd(H.O),+*++ complex. As in many 
neodymium compounds, the anions of the 
chrome alums are separated from the trivalent 
ion by an octahedral shell of water molecules. 

The fundamental frequencies of the sulfate 
ion are about 450, 620, 980 and 1100 cm-. The 
frequencies vary somewhat in different com- 
pounds, and removal of degeneracy in asym- 
metric lattice fields may result in doubling or 
tripling of the energy states.!? Characteristic 
frequencies of the selenate ion are 340, 410, 835 
and 870 cm~. None of the sulfate or selenate 


M4 Van Vleck, J. Chem. Phys. 7, 61, 72 (1939); Phys. 
Rev. 57, 426 (1940). 

% Crystallographic non-identity of all the Cr*** ions 
would lead to a difference of energy of the ground states, 
which, however, would not be revealed necessarily in the 
magnetic susceptibility-temperature curve of the alum. 
The upper states, too, would have different energy so that 
a multiplet could arise in a single transition in different 
kinds of chromium ions as well as in multiple transitions 
to a group of levels in a single kind of chromium ion. 

16 Ewald, Ann. d. Physik 34, 209 (1939). 

7 Hibben, The Raman Effect and its Chemical A ppli- 
cations (Reinhold, 1939), p. 387, 396. 
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frequencies appear prominently in the chrome 
alum spectra. The clearest indications of them 
are in the 15,200-15,600 cm bands of the 
selenates of K, Rb, Tl and Cs and the sulfates 
of Rb, Tl and NH,. Most of the bands in this 
frequency range are common to both selenates 
and sulfates. The ones which are not are broad, 
diffuse and weak. Bands present in selenates and 
absent in sulfates occur shifted 355 and 435 cm 
from the mean frequency of the characteristic 
doublet, while bands present in sulfates and 
absent in selenates are at 450 and 645 cm~. 
Absorption bands surely associated with the 
higher frequencies of the anion have not been 
found. 

The fundamental frequencies of the water 
molecule are 1600, 3600 and 3750 cm~. Because 
of proximity, combination of these frequencies 
with electronic frequencies of Cr++*+ might have 
large probability of occurrence. Corresponding 
absorption bands would lie in the continuum 
and indeed may contribute largely to the 
observed intensity. Other internal vibrations of 
homopolar complexes can occur only in am- 
monium, hydroxylamine and methylamine alums. 
The frequencies are so high (800-3400 cm~!)!8 
that corresponding absorption would also occur 
within or close to the continuum. In intensity 
it would probably be no greater than absorption 
due to internal oscillation of the anions. 

Among the observed weak lines certain of the 
strongest are of interest because they have nearly 
the same separation from the pure electronic 
lines in sulfate and selenate alums of all the 
classes we have distinguished. Such separations 
occur at about 330, 570 and 800 cm~. Sometimes 
the electronic spectrum and its subsidiary 
structure is repeated exactly at a frequency 
330 cm higher, but the diffuseness of bands 
near 570 and 800 cm is such that detailed 
matching of structures is scarcely possible. 
Lines near 195 cm are common to most of the 
spectra. The electronic spectrum is repeated in 
lines of moderate strength and good definition 
at about 110 cm~ in selenate alums of K, Rb, 
Tl and Cs; at 115 cm~ in sulfate alums of Rb 
and Tl; at 120 cm in cesium sulfate and the 


18 Hibben, reference 17, p. 263 ff. 
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Fic. 5. Partial energy level diagram of TlCr(SO.)2-12H,O 
at 20° and 85°K. Vibrational levels associated with both 
the lower and upper electronic states and the transitions 
between them are shown. 


methylamine alums; and at 140 in the hydroxyl- 
amine and ammonium alums. 

Variable low frequency bands (shifted from 
the electronic spectrum by amounts ranging from 
a few wave numbers up to perhaps 100 or 
150 cm-') presumably arise in mutual oscillation 
of massive or loosely bound fragments of the 
lattice, while the higher frequency bands of 
constant separation from the electronic spectrum 
(330, 570 cm=, etc.) may be due to the super- 
position of characteristic librations of a tightly 
bound electrostatic complex (most probably 
Cr(H2O).***) upon the electronic spectrum and 
its attendant low frequency vibrational bands. 
Raman shifts assignable to internal oscillations 
in aquo-complexes have not been observed but 
do appear at frequencies ranging up to 1000 cm= 
or more in metal complexes with ammonia, 
pyridine and ethylene diamine.!® Potassium 
aluminum sulfate has Raman shifts of 189 and 
322 cm“, which probably correspond to the 195 
and 330 cm~ bands we have found. Hellwege”? 


19 Hibben, reference 17, p. 462. 
» Hellwege, Ann. d. Physik 34, 521 (1939). 
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has discovered several reflection maxima in the 
infra-red spectrum of three alums between 56 
and 20u or between 180 and 500 cm. From 
his work it is not at present possible to determine 
absorption maxima for direct comparison with 
our results. Furthermore, since Hellwege’s experi- 
ments were performed at room temperature, 
some of his reflections may represent combina- 
tions of acoustical frequencies (excitable by 
thermal vibrations but not by light) and optical 
frequencies, while in our own experiments at 
20°K acoustical frequencies except of very low 
energy do not come in question. 

A few temperature sensitive lines often appear 
on the long wave-length side of the electronic 
spectrum in the red. Separations of about 115, 
200 and 330 cm~ are found in several alums for 
which Sauer gives data at —78° and +18°. The 
lines arise in transitions from thermally excited 
vibrational states associated with the ground 
state of Cr+++. Because of the low intensity 
and the natural diffuseness of all the absorption 
lines at —78° and higher temperatures the 
separations cannot be very accurately defined, 


_ but their near identity with photo-excited 


vibrations associated with the upper electronic 
state at 20°K is obvious. Similar effects have 
been observed in the absorption and fluorescence 
of Cr+++ dispersed in various foreign lattices.” 
In TICr(SO,)2-12H.O at 85°K we find vibra- 
tional levels associated with the lower state 
19.5, 27, 49, 64, 117.5, 144 and 195.5 cm~ from 
the mean frequency of the doublet. At 20° we 
find upper vibrational states at 19.9, 50, 65.6, 
117, 142 and 195.5 cm—. Figure 5 is a partial 
energy diagram of this alum showing the 
vibrational levels of the lower and upper states 
and the observed transitions among them. 
Three vibrations in the ground state of 
NH;OHCr(SO,)2:12H,O at 20°K have been 
found: »,;=11.4, v2=30, v3=44. Multiples and 
combinations of these frequencies added to the 
two upper electronic states of the characteristic 
doublet, =14,895.4 and e2.=14,918.9, account 
in a satisfactory manner for all of the first 22 
lines of the absorption spectrum (Fig. 6). The 
doublet separation, 23.5 cm@, is repeated 
prominently at 139, 273 and 346 cm™, and in 
the fine structure associated with these lines 
there is distinguishable to a greater or less 
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SPECTRA OF CHROME ALUMS 


extent the vibrational pattern close to the 


doublet. 

Repetition of sequences is observable in all 
the alum spectra. Examples are given graphically 
in Fig. 7. In Fig. 7-D the ten lines comprising 
the (electronic) doublet of TIlCr(SO,)2-12H,O 
and the first eight (vibrational) lines on the 
short wave-length side of the doublet are 
repeated after the 143 cm~ line. Other repeti- 
tions, associated with 118, 196 etc., cannot be 
pointed out with as much assurance of their 
correctness owing to overlapping and to the 
width and diffuseness of the lines. In the corre- 
sponding rubidium alum there is perfect corre- 
spondence among lines appearing in both spectra, 
but matching lines are not always found, pre- 
sumably because the proper combination of 
crystal thickness and exposure time for adequate 
absorption intensity was not arrived at. This is 
especially true immediately to the violet of the 
doublet, where several lines are missing in the 
rubidium alum. The possibility of unequivocal 
establishment of repeating sequences at higher 
frequencies is thus diminished. 
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Fic. 6. Partial energy level diagram of 
NH;OHCr(SO,)2-12H20 at 20°K. 
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Fic. 7. Repetition of line sequences in some of the alum 
spectra. A. NH4Cr(SO,4)2-12H2O (shattered crystal). B. 
NH,Cr(SO,)2-12H2O (clear crystal). C. TiCr(SeO,)2-- 
12H.,0. D. TiCr(SO,)2-12H2O. The intensity of the ab- 
sorption lines is indicated by the height of the lines in the 
= Lines not crossed at the bottom signify absorption 
ines whose centers only were measured. The length of the 
horizontal lines is a measure of the breadth of the absorp- 
tion lines. Horizontal lines upturned slightly at the ends 
signify diffuse absorption bands, whose position and 
breadth is not determinable with high precision. 


It will be noticed in Fig. 7 that perfect repeti- 
tion of the electronic doublet and its closely 
associated vibrational lines is not always achieved 
at the higher frequencies. This may be attribut- 
able in part to inaccuracy in measurements due 
to diffuseness or low intensity of the absorption 
bands. It may also be attributable to an actual 
small shift in the relative positions of the 
vibrational energy levels close to an excited 
electronic level produced by an additional, 
superimposed vibrational jump of, for example, 
195 or 330 cm. A further possibility of ac- 
counting for the observed shift lies in the 
expectation that vibrational frequencies of the 
lattice or its constituents will be somewhat 
dependent upon the electronic state of the Cr++*+ 
ion. The magnitude of the vibrational frequencies 
associated with an electronic jump in a Crt+* 
ion is a function of the masses of the vibrating 
units and the force constants of the bonds 
which link them. The masses are invariant but 
the force constants may be slightly different for 
the various electronic states of Cr+*++ and this 
would lead to small differences in the spacing of 
the vibrational levels superposed upon the 
several electronic states. 
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Using an apparatus similar to the one used by Lozier, the negative ion formation in I; vapor 
was investigated as a function of the electron velocity (0-7 ev). The following reaction is 
responsible for negative ion formation at these velocities: I,+e—~I~+I. Measurements of the 
kinetic energy of the ions for different electron energies show that all the ions are formed by 
the same reaction. The variation of the yield as a function of the velocity of the electrons can 
be understood with the help of a potential-curve diagram, by considering transitions from the 
ground state of the I, molecule to the ground state of the I,~ molecule. The probability that an 
electron colliding with an I, molecule forms a negative ion varies between 10~* and 107%, 
With the aid of an energy cycle one can calculate the electron-affinity of the iodine atom, if one 
knows the dissociation energy of the neutral molecule, the kinetic energy of the incoming 
electron, and the kinetic energy of the ions. The value of the electron affinity obtained by this 
method is 3.0+0.2 ev, which is in agreement with other determinations. 


INTRODUCTION 


HEREAS the formation of positive ions, 

either by electron collision or irradiation 
of light, has been the subject of numerous 
publications! in which one has been able to 
determine the energy necessary to form a positive 
ion and the probability with which the particular 
process occurs, there are only few data available 
in regard to the formation of negative ions.? 
Even for such gases as Cle, Iz, Bre, where the 
electron affinity? of the atom is well known, and 
where the process leading to negative ion forma- 
tion has been investigated, (I.+e—I-+I),’ the 
yield of this reaction and the variation of the 
yield with the electron velocity are not well 
known. The formation of negative ions in iodine 
vapor, which is the subject of this paper, has 
been investigated twice before. One paper is by 
Mohler‘ and the other is by Hey and Leipunski.*® 
The large difference between the values of 
Mohler and Hey and Leipunski in regard to the 
yield does not seem to be a real one. If one 
calculates the yield (given by the expression 


(1) 


where o; and o, are, respectively, the cross 


1 Handbuch der Physik, Vol. 23 (1933). 
wn” Negative Ions (Cambridge, University Press, 
3 Hogness and Harkness, Phys. Rev. 32, 784 (1928). 
4 Mohler, Phys. Rev. 26, 614 (1925). 
5 Hey and Leipunski, Zeits. f. Physik 66, 669 (1930). 


sections for ion formation and total electron 
scattering) using the same electron scattering 
cross section throughout, and if one further 
makes some correction in the evaluation of Hey 
and Leipunski’s data, one obtains values which 
are of the same order of magnitude. The fact 
remains, nevertheless, that the data regarding 
the variation of the yield with the electron 
velocity do not agree and that in both cases it 
is difficult to give a qualitative theoretical 
interpretation of the observed maxima and 
minima. The purpose of this investigation is, 
therefore, to reexamine the variation of the ion 
formation in I, vapor as a function of the electron 
velocity and to measure, at the same time, the 
kinetic energy of the ions. An attempt is then 
made to interpret the results from a theoretical 
point of view. 


DESCRIPTION OF THE APPARATUS 


The apparatus used in this investigation is 
quite similar to the one used by Lozier® with 


v s 


Fic. 1. Schematic representation of the apparatus. 
® Lozier, Phys. Rev. 46, 268 (1934). 
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NEGATIVE ION FORMATION 


which he measured the negative ion formation 
Og and CO. 

Electrons which are emitted by a thoriated 
tungsten filament (F) pass three grids (Pi, Ps 
and P3;) and travel down the axis of the tube 
where they are collected on the plate C (see 
Fig. 1). A solenoid (S) establishing a magnetic 
field parallel to the electron beam fits over the 
tube. When the magnetic field is strong enough 
(about 150 gauss), electrons will not move far 
away from the axis even when they collide with 
gas molecules. On the other hand, the ions formed 
by collision of the electrons with the gas mole- 
cules will hardly be affected by the magnetic 
field and wil therefore move in all directions. 
Some of them will get through the vanes (V)— 
there are twelve of them evenly. distributed over 
a circle—and reach the cylinder Py. The purpose 
of the vanes is a twofold one: (1) To measure the 
kinetic energy of the formed ions. To do this a 
variable retarding potential is applied to the 
outside cylinder. The vanes shield the innermost 
part of the tube sufficiently so that the velocity 
of the electrons remains unchanged. (2) To catch 
those electrons which by multiple reflection or 
by scattering at right angles to their original 
direction might have traveled rather far away 
from the axis. 

The filament was heated by an alternating 
current and was kept at a positive potential with 
respect to the rest of the apparatus when the 
heating current was on. When the heating 
current was off, a negative potential was applied 
to the filament to give the electrons the desired 
velocity, so that the electrons actually emitted 
started from an equal potential surface. The 
second grid was usually kept at a small positive 
potential to overcome part of the space charge 
effects. P3, the vanes and the disk with the small 
opening at the lower end of the vanes were 
always kept at the same potential (at ground) 
so that the electrons were moving in a field- 
free space when they collided with the gas 
molecules. A small positive potential was applied 
to C, to make sure that the total electron current 
actually reached C. 

The electron current measured at C was of the 
order of magnitude of 10-'—10~-* ampere except 
for measurements below 1 volt, where the 
current varied between 10-* and 10-5 ampere. 
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The current to Ps, was measured by a highly 
sensitive galvanometer (lower limit 3.10-" am- 
pere). 

For the few measurements that were made in 
oxygen, no particular precautions were taken to 
get it pure. As a matter of fact, air was used and 
it was assumed that the negative ion formation 
in air is the same as in oxygen. (Experiments by 
Bradbury’ furnish proof for this assumption.) 
For the measurements with Is, purest chemical 


CROSS 
SECTION 


10 ORDINATE UNITS = 


23 4 5 6 7EV 


Fic. 2. Capture cross section plotted against 
energy of electrons in ev. 


iodine available was used; before it was let into 
the apparatus, the iodine was once more sub- 
jected to a vacuum distillation to remove the air 
adsorbed on the crystals. The pressure range 
varied between 6.10-° and 3.10-* mm Hg. In 
both cases (O2 and I) the gas was streaming 
through the apparatus. The pressure of the I, gas 
was measured by a quartz fiber manometer. 
Small changes of the temperature in the cooling 
flask around the iodine container produced rather 
noticeable pressure changes. It should be pointed 
out here that the pressure could be kept constant 
only within certain limits. 


A serious difficulty in the type of apparatus used is the 
uncertainty in the total path traveled by the electrons. Two 
other types of apparatus were constructed, which, however, 
for various reasons, failed to yield satisfactory results. 

In one of these, the electrons originated from a filament 
situated in the axis of a solenoid, passed first through a slit 
in a concentric cylinder and then were bent back by the 


7 Bradbury, Phys. Rev. 44, 883 (1933). 
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Fic. 3. (a) Maximum kinetic energy of negative ions 
lotted against energy of electrons in ev. (b) Ordinate: 
on current-to P, (arbitrary units). Abscissa: Retarding 

potential applied at P,. 


magnetic field to the outside of the same cylinder. Ions 
formed in a field-free space outside of this cylinder were 
collected on a second concentric cylindrical plate of larger 
diameter. The apparatus failed because the ion current was 
masked by the current caused by electrons which succeeded 
in reaching the outer cylinder, probably because of multiple 
collisions or reflection. 

The uncertainty in regard to the value of the mean free 
path is a more serious difficulty. To overcome this point, we 
constructed an apparatus the main part of which is a thick- 
walled iron cylinder with a slot parallel to the axis of the 
cylinder. Inside the iron cylinder is a similarly slit thin- 
walled brass cylinder of slightly smaller diameter. Across 
the slot of the iron cylinder is a magnetic field of about 250 
gauss, strong enough to prevent electrons getting to the 
outside of the iron cylinder. Electrons emitted by a 
tungsten filament and focused by a simple electron lens 
move along the axis inside of the iron cylinder where they 
collide with the molecules. Those electrons scattered in the 
process of collision are collected by the brass cylinder 
whereas the negative ions can easily penetrate the slot and 
are measured by a plate outside of the iron cylinder. The 
ratio of the scattered electron current to the negative ion 
current multiplied by a certain factor due to the geometry 
of the apparatus, will give immediately the absolute yield 
of ion formation without involving a knowledge of the mean 
free path. In order to obtain good results, the stray field 
inside of the magnetic cylinder must be as small as possible. 
It was this particular difficulty which made it impossible 
for the time being to obtain satisfactory results. 


RESULTS 


Figure 2 shows the change in capture cross 
section with electron velocity for I;.* The cross 
section is defined here in the following way: 


I; 
“A No. of molecules/cm* Xd Xf 


(2) 


where d is the distance the electrons travel and 
f is determined by the geometry of the appa- 
ratus. Measurements were only made up to 
about 7.0 ev. At higher velocities positive and 
negative ions are probably formed simultaneously 
and the construction of the apparatus did not 
permit a good analysis of these processes. 

To make sure that the current measured at P, 
was actually caused by ions formed in the gas 
phase and to nothing else, several tests were 
made. 

(1) If one measures the current to P, as a 
function of the magnetic field (everything being 
kept constant) with or without a scatterer in the 
tube, one ought to find that the current to P, 
drops rapidly with increasing field. The speed 
with which it drops depends somewhat on the 
gas pressure, but sooner or later the current to 
P, must become 0, if no ions are formed. On the 
other hand, if a gas is present that forms positive 
or negative ions, the current ought to drop also 
rapidly in the beginning, but then settle down to 
a constant value over a long range. This was 
shown to be the case. 

(2) If we assume that the negative ions are 
formed according to the following reaction 
I,+e-I-+I, it is clear that the maximum 
kinetic energy of the ions for different electron 
velocities must have a definite value which one 
can calculate if one knows the electron affinity 
and the dissociation energy of the neutral mole- 


TABLE I. Maximum yield of negative ions. X2+E—~X-+X 
process leading to negative ion formation. 


Gas OxYGEN IODINE CHLORINE 


Lozier 9.5-10-4 

Bradbury 3-10-% 
Bailey 1.5-10-% 
Hey 5-104 

Mohler 2-10-3 

Present Results 1-10-3 8-10-* 


* This curve is the average of three different runs at 
different pressures. 
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Fic. 4. (a) and (c) Ordinate: o; ——— 

(arbitrary units). Abscissa: Energy of 
electrons in ev. (b) and (d) Ordinate: 
Volts. Abscissa: Internuclear distance 

in angstrom units. 
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(c) 
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(d) 


'23 4 5 6 7etv 


cule. All our measurements show that the ions 
have a kinetic energy which corresponds to the 
dissociation of the molecular ion into a negative 
atomic ion and an unexcited* neutral atom. 
(See Fig. 3(a).) If one extrapolates the straight 
part of the curve in Fig. 3(a) where it cuts the 
x axis, one can find the electron affinity of the 
iodine atom. The value obtained in this manner 
is 3.0+0.2 which agrees with other determina- 
tions by Mayer® and Glockler.® At this point the 
question concerning the formation of ions at the 
filament might arise. It is well known that the I, 
molecule dissociates into atoms on a hot tungsten 
filament above a certain temperature. Negative 
ions are then formed easily. These ions would 
then gain the same energy as the electrons, 
might make collision with neutral I, molecules 
and reach Py. The possibility that something 
like this happens cannot be completely excluded 
offhand. Nevertheless the geometry of the appa- 
ratus makes it quite improbable that many of 
those ions would reach P,. If so, their maximum 
kinetic energy ought to be equal to the total 
energy they have gained. For example, with 
V;=—6 one ought to find some particles reach- 


* The possibility that some of the collisions will yield a 
gaa atomic ion plus an excited atom will be discussed 
ater. 

8 Mayer and Sutton, J. Chem. Phys. 3, 20 (1935). 

® Glockler and Galvin, J. Chem. Phys. 3, 771 (1935). 


ing P, which possess 6-volt kinetic energy or a 
little less. Figure 3(b) shows that this was never 
the case. 

(3) Our results for oxygen (in the range of 
0-10 ev) agree well with the experiments done 
by Lozier. 


DISCUSSION 


a. Absolute value of yield 


To calculate the absolute value of the yield as 
defined in Eq. (1), we must first know the cross 
section for ion formation, o; given by the ex- 
pression (2). Nevertheless, the value obtained in 
this manner can only be an approximate one, as 
we do not know the exact value for d* (the 
actual path length of the electrons), because of 
the fact that all the electrons that have once 
made collision with a molecule will move on 
helixes of different sizes instead of straight lines. 
This remark applies to Mohler’s data also and 
in a similar way to Hey and Leipunski’s results. 
In order to compare our results with the data of 
Mohler and Hey and Leipunski, we made use 
of data given by Fisk" on the cross section of 
total electron scattering in chlorine, assuming 


* The value of d which was chosen for the evaluation of 
Eq. (2) was the distance between P; and the lower end of 
the vanes. One thus obtains a maximum value for o;. 

10 Fisk, Phys. Rev. 51, 25 (1937). 
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that the cross section for iodine is certainly not 
smaller than that of chlorine. Q=1500 cm~ (at 
1 mm Hg) was used as an average value for the 
cross section o;. If one uses this value to calculate 
the maximum yield, one gets values which are of 
the same order of magnitude for all three investi- 
gations, as one may see from Table I. For com- 
parison, Table I also gives the values for Cl. and 
for O2 when the ion formation is caused by the 
same process: (X2+e—-X-+X). In Os: this 
process can only take place with electrons of a 
certain amount of energy as the electron affinity 
is about 2.8 ev smaller than the dissociation of 
the neutral molecule. One sees, therefore, that 
for all processes of negative ion formation which 
proceed according to the above-mentioned equa- 
tion, the yield is about the same. On the other 
hand, the accuracy of the absolute values is 
small. 


b. Variation of the yield with the electron 
velocity 


The results of the different authors regarding 
the variation of the yield with the electron 
velocity cannot be made to agree satisfactorily. 
Without discussing the details of the apparatus 
used by Hey and Leipunski, we shall merely 
remark that their arrangement does not seem to 
give good results for the processes due to very 
slow electrons, mainly because the effect of the 
earth’s magnetic field on these slow electrons 
cannot be compensated for in a satisfactory way. 
In regard to Mohler’s paper, one must say that 
because of the high pressure of iodine vapor 
(10-*—10-? mm Hg) used in his investigation, 
multiple scattering sets in which might change 
the shape of the yield curve to a considerable 
extent. On the other hand, it seems to us that 
the curves measured with our apparatus could 
hardly be false. We, therefore, limit ourselves 
to a theoretical discussion of the result repro- 
duced in Fig. 2. 

The reaction I,+e—-I-+1+E xin) is energeti- 
cally possible for electrons with zero velocity 
because the dissociation energy of the iodine 
molecule is considerably smaller than the elec- 
tron affinity of the iodine atom. The difference 
between the electron affinity (A) and the disso- 
ciation energy (D) must be taken up as kinetic 


energy of the iodine ion and the iodine atom. 
For electrons with zero velocity the kinetic 
energy of the ion is 7,=D—A/2.* If we tempo- 
rarily disregard the secondary maxima at 1.75 
and 2.5 ev (and possibly at 3.5 ev), we get a 
curve like the one shown in Fig. 4(a). Electrons 
with zero velocity give indeed a considerable 
yield of ion formation according to the equation 
I,+e—I-+I, but the yield increases with the 
velocity of the electrons up to 0.4 ev (+0.15 ev). 
We have no reason to doubt the reality of that 
maximum. For electrons with a velocity higher 
than 0.4 ev the yield drops sharply but still 
gives a measurable yield up to 7 ev (the maxi- 
mum electron velocity used for the investigation 
in iodine vapor). Bradbury," using the filter 
method due to Loeb,” found a maximum of ion 
formation in chlorine at 0.75 ev. His measure- 
ments do not extend to higher electron velocities. 

One can explain this curve in a simple way 
with the help of the Franck-Condon principle, 
if one assumes the potential curve of I~ to be 
as shown in Fig. 4(b). The dotted part of the 
curve is quite uncertain. Following considera- 
tions due to Mulliken, one must conclude that 
the I;- molecule has a definite dissociation 
energy, but the actual depth of that minimum 
and its accurate position are unimportant for 
our considerations. On the other hand, the steep 
part of the curve and the point where this 
curve cuts the curve for the normal I; molecule 
are essential for our interpretation. The maxi- 
mum at 0.4 ev is then caused by a transition 
from I.+e—I,- corresponding to the perpen- 
dicular distance from the minimum of the I, 
curve to the I;- curve. The rather extended de- 
crease of the curve following the maximum at 
0.4 ev indicates that the repulsive part of the 
I;- potential curve must be very steep in that 
region so that the I,- curve runs for a long time 
inside of the shaded area which corresponds to 
the zero-point vibration of the I, molecule. We 
are here confronted with potential curves which 
are similar to those which are responsible for the 
appearance of long continuous absorption spectra 


* D=1.54 ev; W. G. Brown, Phys. Rev. 38, 709 (1931). 
1 Bradbury, J. Chem. Phys. 2, 835 (1934). 
22 Loeb, Phys. Rev. 48, 684 (1935). 
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of molecules. (The fact that the iodine vapor was 
usually at a temperature of about 60°C because 
of the heating effect of the solenoid, naturally 
brings a small number of iodine molecules into 
their first vibrational state; but in first approxi- 
mation we can neglect that fact.) 

We now turn to the discussion of the maxima 
at 1.75 and 2.5 ev (and possibly at 3.5 ev, our 
results do not indicate definitely whether or not 
there is a maximum). (See Fig. 4(c).) Both these 
values (1.75 and 2.5 ev) agree quantitatively 
with the perpendicular distance from the ground 
state of the normal I; molecule to the first 
excited states, as one can see from the potential 
curve diagram (Fig. 4(d)) which is known from 
spectroscopic work.'* The transition from the 
ground state of I, to the strong repulsive curve 
of the iodine molecule is about 3.5 ev." 

It was obvious to presume that the maxima 
at 1.75 and 2.5 ev were caused by inelastic 
collisions of the electrons with the iodine mole- 
cules. Two possibilities exist: (1) Because of the 
inelastic impacts of the electrons the velocity 
distribution of the primary electrons is changed 
to such an extent that besides electrons with a 
velocity of X volts there is at the same time an 
appreciable amount of electrons present with a 
velocity of X—1.75 and X—2.5 volts. As the 
distance from the filament to the electron col- 
lector is of the same order of magnitude as the 
mean free path, one should have an appreciable 
amount of these electrons providing the proba- 
bility of excitation is very great in the region of 
1.75 and 2.5 volts. As a matter of fact, this 
might not be impossible. But if these secondary 
maxima are caused by electrons which already 
have been slowed down, they should not lie just 
at the minimum excitation potential, but should 
be shifted by an amount of 0.4 ev; because 
electrons of 0.4 ev have a higher probability to 
form negative ions than electrons of zero velocity. 
This reasoning is only justified if we accept 
definitely the first maximum as given. But even 
if one doubts this, one has difficulties in in- 
terpreting the measurements of the velocity 
_ distribution of the ions which are formed. If, for 


8 Brown, Phys. Rev. 38, 1187 (1931). 
4 Turner, Zeits. f. Physik 65, 464 (1930). 
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example, the maximum at 2.5 ev is caused by 
electrons of a velocity of X —2.5 volts, then the 
number of these slowed-down electrons should 
increase, with increasing electron velocity accord- 
ing to the excitation function. A measurement 
of the velocity distribution of the ions formed 
with electrons of 3 ev should, therefore, contain 
two distinct groups of ions of different velocities, 
one with a maximum velocity of 2.35 volts and 
another cne with a maximum velocity of 1.7 
volts. (Fig. 3(b) shows the results for ions 
formed with electrons of 3 ev.) Because of the 
geometry of the apparatus, an exact evaluation 
of the actual energy distribution of the ions is 
not possible. Nevertheless it is possible to decide 
whether or not there are different groups of a 
definite velocity present. The differentiated 
curves will indicate these groups even more 
distinctly. The monotonic shape of these curves 
seems to be a strong argument against the just- 
discussed possibility. We therefore come to the 
conclusion that the secondary maxima are caused 
by a process which occurs in one elementary act. 
The possibility that this can occur does exist, 
as the steep potential curve of the iodine ion 
molecule cuts the potential curves of the excited 
neutral molecule at points which lie practically 
perpendicular above the ground state of the 
neutral molecule, as it is indicated in Fig. 4(d). 
This would mean therefore that the same electron 
which excites or dissociates the iodine molecule 
is caught by one of the iodine atoms. The mono- 
tonic shape of the decrease of the yield with 
increasing electron velocity will be interrupted 
at those places where the potential curve of I,- 
cuts the potential curve of the excited iodine 
molecule. It is important that these cutting 
points lie directly above the ground state. At 
2.5 ev and above the reaction I,+I-—I-+I* 
might occur, where I* is an excited iodine atom. 
The yield for this reaction will be considerably 
smaller. 

Our interpretation could be checked experi- 
mentally if the dependence of the yield could be 
studied as a function of the pressure over a wide 
pressure range. Following the interpretation 


which we prefer, one ought to find that the rela- 


tive strength of the maxima is independent of the 
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pressure, whereas according to the first-men- 
tioned and rejected interpretation, they ought to 
increase quadratically with the pressure. The 
apparatus used in this investigation unfortu- 
nately did not allow us to make measurements 
at very low pressure. 


I am greatly indebted to Professor James 
Franck, under whose direction this investigation 
was started, for his continued help and advice. 
It is a pleasure to acknowledge many helpful 
discussions with Professor J. E. Mayer and 
Professor J. A. Bearden. 
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The Differential Diffusion Constant of Calcium Chloride in Aqueous Solution 


E. A. HoLLINGsHEAD AND A. R. GoRDON 
Chemistry Department, University of Toronto, Toronto, Ontario, Canada 


(Received October 25, 1940) 


The differential diffusion constant for calcium chloride in aqueous solution has been measured 
for concentrations ranging from 0.0025 to 0.5 mole per liter and for temperatures from 15° to 
35°C. The relation connecting the diffusion constant at a finite concentration with that at 
infinite dilution, k=o(1+cd In f/dc)(mo/n), which holds for moderate concentrations for 
several 1—1 electrolytes, is valid for this salt only up to 0.01 mole per liter; the corresponding 
value of ky is about 1.5 percent less than that computed from the ionic mobilities by the Nernst 
relation. The activation energy of diffusion for a given concentration decreases with rising . 
temperature as is customary; for a given temperature, however, it decreases with increasing 
concentration in dilute solutions, and finally becomes independent of concentration. 


REVIOUS papers! from this laboratory have 

reported the differential diffusion constants 

for hydrochloric and sulphuric acids in aqueous 

solution as functions of temperature and con- 

centration. The present work was undertaken 

to obtain similar data for a typical salt of the 
2—1 valence type viz. calcium chloride. 

As in the earlier measurements, diffusion cells 
of the Northrop-McBain type were employed, 
the technique of the measurements and the 
method of analysis of the solutions by inter- 
ferometer being the same as were used in the 
previous work. Stock calcium chloride solutions 
(made up from the twice recrystallized salt and 
conductivity water) were analyzed by precipita- 
tion of silver chloride; all solutions used in the 
diffusion measurements as well as all reference 
solutions for the interferometer were made by 
gravimetric dilution of these stock solutions; 
weight concentrations were converted to volume 
concentrations by means of the density data in 


1 (a) W. A. James, E. A. Hollingshead and A. R. Gordon, 
J. Chem. Phys. 7, 89 (1939); (b) W. A. James and A. R. 
Gordon, ibid. 7, 963 (1939); (c) E. A. Hollingshead and 
A. R. Gordon, ibid. 8, 423 (1940). 


the International Critical Tables. As in the 
previous measurements, the cells were recali- 
brated after every three or four calcium chloride 
runs by diffusing 0.1 KCI solution in water at 
25°C; the value of & used in the calibration'* 
was 1.842 X 10-5 cm?/sec. 

Table I contains the results; the second and 
third columns give the average initial concen- 
trations of the ‘‘inner’’ and ‘‘outer” solutions 
expressed in moles per liter. The next column 
gives the square root of the mean concentration 
of the experiment—the concentration for which 
the effective diffusion constant k (computed by 
Eq. (2) of reference 1a) is equal to the true 
differential diffusion constant.!*? The last column 
gives the average k, obtained in each case from 
measurements in four different cells; the indi- 
vidual values showed in general a mean absolute 
deviation from the mean of somewhat less than 
0.5 percent. 

Figure 1 shows the results graphically, the 


continuous curves of the figure corresponding to . 


the smoothed entries of Table II. Adequate 
2 A. R. Gordon, J. Chem. Phys. 5, 522 (1937). 
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DIFFERENTIAL DIFFUSION CONSTANT 


thermodynamic data to test the relation previ- 
ously proposed? for 1—1 salts 


k=ko(1+cd In f/dc)- (n0/n) 


are only available for calcium chloride for 25°C 
and for c=0.1. The dotted curve of the figure 
corresponds to this relation, the thermodynamic 
data being those of Shedlovsky and MacInnes? 
and the viscosities from the International 
Critical Tables. The value selected for ko is 
1.312 X10-5, slightly less than that computed by 
the Nernst relation from the ionic mobilities’ 
viz. 1.335X10-*. The figure shows that there is 
appreciable deviation for c>0.01 ; this is perhaps 
not surprising when it is remembered that the 
ionic mobilities of this salt, even for low concen- 
trations, do not obey the rules valid for 1—1 
electrolytes. 

From Table II it is possible to obtain values 
of the activation energy of diffusion; at the 
lowest concentration (c=0.0025) the mean 
activation energies for the temperature differ- 
ences 15—25° and 25-35° are 4.72 and 4.37 kcal., 


respectively; with rising concentration the 
activation energy falls until when c20.04, it 
TABLE I. 
T°C co’ co” ((co’ +c0")/2)8 108k 
15 0.00565 0 0.053 0.904 
15 .01025 0 .072 .905 
15 -0500 0 .158 
15 .1522 0.0515 319 .870 
15 3001 .2023 874 
15 .5023 4010 .672 
25 .00513 0 .051 1.201 
25 .01024 0 .072 1.174 
25 .02060 0 101 1.161 
25 .0483 0 155 1.154 
25 .0510 .318 1.129 
25 .3020 .1982 .500 1.138 
25 .5043 4034 .674 1.156 
35 .00559 0 .053 1.538 
35 .01045 0 .072 1.488 
35 .0488 0 .156 1.449 
35 1491 .0368 .305 1.423 
35 .2999 .2002 -500 1.431 
35 4962 .4000 .669 1.456 


’ T. Shedlovsky and D. A. MacInnes, J. Am. Chem. Soc. 
59, 503 (1937). 

‘Int. Crit. Tab. Vol. 5, 14 (1929). 

5D. A. MacInnes, J. Frank. Inst. 225, 661 (1938). 

°L. G. Longsworth and D. A. MacInnes, J. Am. Chem. 
Soc. 60, 3070 71938). 
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TABLE II. 10° k for aqueous solutions of CaCle. 


a 15° 25° 35° 
0.05 0.913 1.204 1.530 
893 1.162 1.475 
20 876 1.135 1.433 
30 871 1.128 1.422 
50 878 1.138 1.432 
70 898 1.162 1.462 


becomes practically independent of concentra- 
tion at 4.4 and 4.2 kcal., respectively. Thus the 
activation energy falls with rising temperature 
as in the other cases studied."’ Its dependence 
on concentration, however, is different from that 
found for the acids; for hydrochloric, the energy 
of activation did not depend on concentration 
up to 0.2 normal, but then increased, while for 
sulphuric the activation energy was practically 
constant for solutions varying from 0.01 to 1.0 
mole per liter. This would seem to indicate that 
the change in the diffusional energy of activation 
with concentration as well as its absolute value, 
is, for a given solvent and temperature, a specific 
property of the diffusing substance. 


7H. S. Taylor, J. Chem. Phys. 6, 331 (1938). 
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Using a potential function of the form 


Ce C3 


Ro 


the vibrational wave number and energy of dissociation of the molecules HgHe, HgNe, HgA, 
HgKr, HgXe, Hge, (O2z)2 and (NO)s in the normal state are calculated. It is found that, if the 
single value p=0.28A is chosen, satisfactory agreement is obtained with experimental data 
except in the case of (NO)». The latter structure is believed to owe its binding partly to first- 


order attractive forces. 


INTRODUCTION 


PECTROSCOPIC research as well as results 
from many other sources strongly indicate 
the existence of quasi-molecules, such as HgA or 
(O2)2, which are not molecules in the ordinary 
sense because, at least in the gaseous state and 
under normal conditions, their lifetime is ex- 
ceedingly small and their concentration minute. 
These molecules are considered to owe their 
existence to the so-called van der Waals forces 
and are therefore spoken of as van der Waals 
molecules. They show very flat potential minima 
(depth A=0.01+0.1 ev) at relatively large 
interatomic distances 

The subject of van der Waals forces has been 
thoroughly dealt with recently in a paper by 
London! and in Margenau’s? review article. 

The present investigation is concerned with the 
eight molecules HgHe, HgNe, HgA, HgkKr, 
HgXe, Hge, (O2)2 and (NO). in the gaseous 
phase and the lowest energy state. The two last- 
named molecules can be safely treated as con- 
sisting of two bodies inasmuch as the two atoms 
in each normal O.[NO] molecule are fairly 
tightly bound (energy of dissociation slightly 
above 5 ev*) and their separation of 1.21[1.15 JA 
is much smaller than the intermolecular distance, 
Ro say, of (Oz)2[ (NO). ]. Disregarding the contri- 
butions of permanent quadrupole and higher 


* Part of a dissertation presented to the Faculty of the 
Graduate School of Yale University in candidacy for the 
degree of Doctor of Philosophy. 

1F, London, Trans. Faraday Soc. 33, 8 (1937). 

2H. Margenau, Rev. Mod. Phys. 11, 1 (1939). 

3G. Herzberg, Molecular Spectra and Molecular Struc- 
ture, I (New York, 1939), Table 36. 


moments of Hg, O2 and NO, whose magnitude is 
not known at present but is probably small, the 
dipole-dipole terms in the alignment and induc- 
tion effects are utterly negligible as compared to 
the dispersion effect in the case of nitric oxide, 
whose permanent dipole moment is somewhat 
less than 0.2(10)-!8 e.s.u.,4 and they are zero for 
the other molecules, whose constituents carry no 
permanent dipole moments. Thus only the dis- 
persion forces will be considered in the attractive 
part of the potential. 

In addition to the second-order attractive 
forces, the first-order contribution must be taken 
into account. In London’s' origina! calculation of 
the depth of the potential minimum the repulsive 
potential was approximated by a vertical line 
drawn at a distance equal to the sum of the gas 
kinetic radii of the interacting atoms. The depth 
of the potential at the minimum was also com- 
puted by Ekstein and Magat,® who improved on 
London’s calculations by employing a term DR~ 
to represent the repulsive potential; and the 
same type of potential function was used by 
Wohl’ and by Salow.* We shall, following Born 
and Mayer’s® treatment of the alkali-halide 
crystals, write the repulsive potential in the 

‘Cc. P. a and K. B. McAlpine, J. Chem. Phys. 1, 
60 (1933); E. Watson, G. Gundo Rao and K. L. 
Ramaswamy, Proc. Roy. Soc. London A143, 558 (1934). 

5F. London, Zeits. f. physik. Chemie Bil, 222 (1930) 
[Hg:, HgA, HgKr}. 

®H. Ekstein and M. Magat, Comptes rendus 199, 264 


(1934) [Hg]. 
—e Zeits. f. physik. Chemie B14, 36 (1931) 
‘H. Salow, Dissertation (Berlin, 1935), pp. 40-48 


[(O2)2]. 
? M. Born and J. E. Mayer, Zeits. f. Physik 75, 1 (1932). 
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THEORY OF SOME VAN DER WAALS MOLECULES 


form Ae-*/*, where A and p are positive con- 
stants. Detailed quantum-mechanical calcula- 
tions” indicate that such an expression may be 
expected to be a fair approximation in a certain 
region about the potential minimum. 


That the first-order potential energy will always lead to 
repulsion is by no means evident as far as the molecules 
(Oz)2 and (NO)2 are concerned, their (paramagnetic !) 
constituents possessing unsaturated spin moments (normal 
state: *& and II, respectively). In the absence of theoretical 
information as to the size of the first-order attractive forces, 
however, they will be considered to be negligible. As regards 
the validity of this assumption, the reader is referred to the 
last paragraph in the section on the molecule (NO)>. 


THE POTENTIAL ENERGY FuNcTION V(R) 


Denote by R the (variable) distance between 
the centers of positive charge of the two normal- 
state atoms or groups of atoms forming a van 
der Waals molecule. Their mutual potential 
energy, V(R), will be the sum of the repulsive 
potential, Ae-*/*, and the attractive dispersion 
potential, AE® say. 

For twoelectrically neutral charge distributions 
I and JI representing normal atoms or molecules 
the dispersion-force potential is of the form 


Ci 
AE® = — . (1) 
R2i+2) 
The coefficient of the leading term, —c,/R®, 


which represents the dipole-dipole interaction, is 
given by London’s general expression,® which for 
two atoms or groups of atoms in the normal state 
reads 


3 4h 


3274 mo? + 9/4) 


where .e and mp are charge and mass of the 
electron, the summation is extended over all 
energy states k (for atom J) and / (for atom JJ), 
and the »’s and f’s are the frequencies and 
oscillator strengths associated with the transi- 
tions from the ground state to the excited state 
labeled k or J and are available from the most 
complete experimental Sellmeier-type dispersion 
formulas for substance J or JJ. Approximate 


0L. Pauling, Zeits. f. Krist. 67, = (1928); H. Briick, 
Zeits. f. Physik 4) ‘ag (1928): P . Gombas, Zeits ft. 
Physik 93, 378 (19 (1934 
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values for the coefficients of the higher powers of 
R-~ in the expansion (1) are found by means of a 
perturbation calculation, using Margenau’s"™ 
harmonic oscillator model, the result being 


45 1 a! hv 
63 =— — 
8 fi 2hv?+hr™ 
fl +hv! 
105 1 
fi 7 3hv' 


3 ag! al! 


2 fl fl + 


(3) 


(hv") 


2 
4725 1 (hy')8 
32 4hy' 
ag! aol! ao! hy! 


+2— 
fl fll ft 3hv! 


Here v’, f/ and ao’(j=T, IJ) are the characteristic 
frequency, oscillator strength and static molecular 
polarizability of atom j as given by experimental 
one-term Sellmeier-type dispersion formulas. For 
a medium with a one-term dispersion formula 
there exists the relation 


(hv!) 
(5) 


+hy! 


e? 


i= 


(j=1, IN). (6) 


(v')? 


The coefficients of higher powers in R~* than R-” 
in (1) are neglected. 


The derivation of (3), (4), and (5) is similar to 
Margenau’s," the difference being that the two constituents 
I and II of the molecule will here be allowed to be different 
and that higher terms will be included in the classical 
Coulombian interaction energy. For two neutral charge 


11H. Margenau, J. Chem. Phys. 6, 896 (1938). 
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TABLE I. List of the principal numerical data. (p=0.28A.) 


WAVE ZERO- 
NUMBER POINT 
a C2 Ro Rm p* A ENERGY | ENERGY OF DissocIATIONt 
MOLECULE |(EvV XA‘) (Ev XA8) (EV XA!) (A) (A) (A) (Ev) (cm~) (1073 Ev) (1073 Ev) 
Hg He 9.12 33.85 130.4 623.9 3.28 3.1 0.3627 689.7| 46.6 2.88 (€0) | 3.01 (Do) 
2.53 (e1) 3.36 (Di) 
Hg Ne 17.62 65.37 251.6 1202 3.09 3.0 0.3459 1237 31.5 1.95 (eo) | 11.22 (Do) 
1.88 (e1) | 11.30 (Di) 
Hg A 66.09 257.8 1055 5409 3.39 3.35 0.3998 6601 | 35.9 3007a) | 2.20 (eo) | 27.5 (Do) 
2.16 (a) | 27.5 (Di) 63(d)® 
Hg Kr 95.80 380.2 1624 8714 3.46 3.4 0.4056 10275 31.7 30%) | 1.96 (eo) | 38.1 (Do) 35@7=) 
1.93 (a) | 38.1 (D1) 75(A)® 
Hg Xe 150.6 610.1 2759 15687 3.66 3.65 0.4442 22375 28.1 1.74 (e0) | 41.9 (Do) 
1.72 (ea) | 41.9 (D1) 
Hg 143.3 781.2 4049 26430 3.20 3.2 0.3512 16133 | 35.2 3605) | 2.18 (eo) | 86.0 (Do) 66+91°) 
2.16 | 86.0 (D1) 
84(A) (6) 
(O2)2 34.69 70.54 232.5 883.4 3.48 3.4 0.4397 2806 | 34.4 48%) 2.13 (€0) | 10.77 (Do) 5.6%) 
2.04 (a) | 10.86 (Di) 10%) 
p. 62] 
10(Do 
(NO): 34.66 83.28 268.9 1094 3.68 3.6 0.4662 3857 | 30.5 1.89 (0) 7.90 (Do) —_ p. 1792] 
1.80 (a) | 7.99 (Di) 3 


t References e and g are to footnotes to Table II. 


distributions each of which has its positive charge localized 
at one point this energy, up to and including the term in 
R“, is 


U!(R) = Upp + Upot+(Uoo+ Upa)+(Ueat+ Ups), (7) 


where, in a form convenient for the computation, 
e 
3 
st 


4 s Pt s Pt 
— 152,75 .?+2(r,- S201)? ]; 
1 
+5 (22+) 7261) J; 
2 
wt 8 Re = — 6 — (2r,- 72.5 


+7 (2,3 (4r, 9261) 


R points from J to IJ, (xs, Vs, 22) (Es, ne, £1) ] is the 
position vector of the sth [/th] electron e of multipole 
I(IT) referred to the respective center as origin, and the 
z’s and ¢’s are taken along R. The subscripts p, 9, 2 and ® 
affixed to the U’s refer to the (instantaneous) dipole, 
quadrupole, octupole and fourth-order pole moments whose 
interaction is represented by these components of (7). 
When the matrix elements of (7) are calculated with 
the use of oscillator eigenfunctions, employing different 
stiffness constants (41°m/h)vi=6i(7=TI, IJ) for the two 


constituents of the molecule, and then substituted into the 
formula for the second-order perturbation energy of the 
normal state and finally the value of 8’ is expressed by 
means of (6), there results 


co 
R® R& RO RI’ 


complete up to and including the term in R™!*, with 2, cs 
and ¢& given by (3), (4) and (5), whereas the coefficient 


c= vad ao (8) 


already derived by London,‘ is the special case of the 
accurate expression (2) for substances with one-term 
dispersion formulas. 

The numerical factors in the energy denominators of 
the coefficients (8), (3), (4), and (5) indicate in an obvious 
way the combination of poles which has given rise to the 
corresponding term. Note that, owing to the inclusion of the 
dipole-octupole interactions, the coefficient (4) of the term 
in R-” in the special case of two alike atoms or molecules 
is 7/3 times the corresponding coefficient as found by 
Margenau," who did not consider the interactions involving 
poles higher than quadrupoles. It will be seen that the 
effect of the term in R~!2, here included for compuaianess, 
is small. 


In Table I will be found the values of the 
coefficients c; to c, given by (2) to (5) as com- 
puted from the following experimental dispersion 
formulas: the three-term formulas of Wolfsohn” 
for mercury, Larsén™ for argon, Ladenburg and 


2G. Wolfsohn, Zeits. f. Physik 83, 234 (1933). 
13 T. Larsén, Zeits. f. Physik 88, 389 (1934). 
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THEORY OF SOME VAN DER WAALS MOLECULES 


Wolfsohn" for oxygen; the two-term formulas of 
Kronjager’ for krypton and xenon, and Koch!* 
for nitric oxide ; the one-term formulas of Quarder, 
Herzfeld and Wolf’ for argon and of C. and M. 
Cuthbertson!* for mercury and oxygen, krypton 
and xenon, helium and neon and nitric oxide. 

The potential energy function of a van der 
Waals molecule is thus 


R® RO RY 


It will have a minimum at some interatomic 
distance R,, if and only if 


45 46%). Rmlp 
R,4 = ® (10) 
=A(p) 
and 
4 Co 5 C3 C4 
3K,‘ my (11) 
p< p”. 
12 co 55 26 7 
7 Ry? 21 7 


A more detailed study shows that (since p and 
the c’s are all positive) the equality sign in (11) 
cannot be admitted as it would lead to a point of 
inflection instead of a minimum, also that, pro- 
viding (10) and (11) hold, (9) has its maximum 
at some point Ry <R,, and has no other (finite) 
R where the tangent is horizontal. Since the 
values of the parameter p that we will have 
occasion to use are about ? times the maximum 
value p* defined by (11) and hence small enough 
to permit the maximum of V(R) to rise fairly 
high above the R axis, the potential function has 
the shape shown in Fig. 1. It is sufficient for our 
purpose and will be assumed that V(R) is trust- 
worthy over the classical range of vibration for 
the lowest vibrational level and thence up to 
infinite separation. The depth A of the potential 
42 (193 ee and G. Wolfsohn, Zeits. f. Physik 79, 

* W. Kronjager, Zeits. f. Physik 98, 17 (1935). 

oa Koch, Arkiv f. Matematik 18, No. 3 (1924). 

B. Quarder, Ann. d. Physik [4], 74, 255 (1924); K. F. 

Heraeld and K. L. Wolf, ibid. [4], 76, 1 (1925). 
18 C. and M. Cuthbertson: Proc. Roy. Soc. London A83, 


151 (1909) [O., Hg]; A84, 13 (1910) Kr, Xe]; A135, 40 
tNOL [He, Ne]; Trans. Roy. Soc. London A213, 1 (1913) 
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(9) at the minimum is 


C2 8p 
+2 (1-—) 
Rn 
10 12 
R,,"* | 


and the second derivative at the minimum is 


2 7p Co 9p 
+4 (1-*) 
R.,* Rn 


C3 lip C4 13p 
+5—{ 1-—— +6—(1-— | (13) 
Rn R,,* Ra 


p being bounded by (11) and the subtrahends in 
the parentheses expressing the effect of the 
repulsive part of the potential energy. 


Rn 


I 


V"(Rn) = 


THE ZERO-POINT VIBRATIONAL ENERGY 


Let vo stand for the classical frequency of 
vibrations with infinitely small amplitude of a 
diatomic molecule about the equilibrium posi- 
tion, which then coincides with the separation 
R,» corresponding to the minimum of the po- 
tential function. As V’’(R») =(2mv0)?, where is 
the molecule’s reduced mass, the wave number 
associated with yo is 


1 (—) 
c 2nxc 


with V’’(Rm) given by (13). 
In zeroth approximation the potential function 
(9) can be replaced by a Hooke’s law potential. 


(14) 


>R 


Ry Rn 


Fic. 1. The potential function V(R). 
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The pure vibration energy of the molecule in the 
lowest energy state is then known to be 


thvp=she- Vo (15) 


with the small amplitude wave number 7 
from (14). 

The potential function (9) is not a parabola, 
however; and it will be shown in Fig. 3 for the 
case of the quasi-molecule HgHe that the Morse 
potential function M(R) which best fits our 
function V(R) represents its asymmetry over a 
considerable range of the variable R around the 
minimum value R, very well. One may thus 
expect the zero-point vibrational energy, «1 say, 
on the Morse potential to be a better approxi- 
mation than the oscillator energy (15), in terms 
of which it may be written 


1 €0 
4A 


where A is the depth (12) of the potential mini- 
mum and the subtrahend in the parenthesis 
expresses the effect on the vibrational energy of 
the asymmetry of V(R) about R,,,—an effect that 
is almost negligible for all the molecules here 
considered except HgHe. 


DETERMINATION OF Ry» AND p 


The energy of dissociation of a van der Waals 
molecule in the state of lowest electronic and 
vibrational energies is 


A—e=Do (17) 
in zeroth approximation, and 
A- = D 1 (18) 


in first approximation, A standing for the depth 
(12) of the potential at the minimum. 

The quantities vo, €0, €1, Do, D1, as well as the 
potential curve (9, 10) and the Morse curve best 
fitting it, depend on the coefficients ci to Cu, 
which have already been expressed and computed 
as functions of the constants of experimental 
dispersion formulas; on the molecule’s reduced 
mass u; on the separation R,,, by means of which 
the constant A(p) could be eliminated in the 
potential ; and on p. 

The asymmetry of the potential curve V(R), 
already responsible for the slight depression of 


the vibrational level ¢; as against €9, makes itself 
felt, moreover, in that the distance R,», corre- 
sponding to the minimum of the potential func- 
tion V(R), cannot be expected exactly to equal 
the equilibrium separation, Rosay, of the molecule 
in its normal state, though the relative difference 
should be small. To correlate R,, and Ro, a 
graphic procedure, suggested by Kane,’® will be 
employed. He identified the midpoint of the 
classical range of oscillation of the lowest vibra- 
tional level with the equilibrium distance Ro. 
(See Fig. 1.) So we will have to choose R,, such 
that, taking some value of p satisfying the 
minimum condition (11), that midpoint comes 
just at the equilibrium distance Ro. We shall feel 
contented, however, when the difference between 
midpoint and Ro does not exceed one percent 
of Ro. 

The equilibritim distance Ry ought to be de- 
duced from the spectroscopic constant for rotation 
in the lowest vibrational state, By=h/(8m*uRo’). 
Unfortunately, however, the rotational structure 
of the band spectra of the van der Waals 
molecules under discussion is missing or (in the 
case of Hg) unresolved, and hence the Bo 
constants are not known. To find approximate 
values for Ro for the various molecules, recourse 
will be had to crystal and related data. It is to be 
hoped that these Rp values are not much in error. 

The values adopted for Rp are listed in Table I. 
For the five mercury-rare gas compounds, we 
shall take the arithmetic mean of the distances 
between two nearest neighbors at or near 0°K in 
the crystal lattices of mercury, 2.98A as found by 
extrapolation with the aid of density measure- 
ments” of the value given by Neuburger,” and of 
the inert gases; for helium this is*® 3.57A (at 
1.45°K), whereas the values used for the re- 
maining rare gases differ by less than one 
percent from those employed by Kane.'® Mercury 
being a metal, the true equilibrium distance of 
two normal Hg atoms forming an isolated Hg» 

1G. Kane, J. Chem. Phys. 7, 603 (1939). [Rm and Ro, 
p. 607; Ro values, first row of Table V and p. 605, where 
reference d ought to read: W. H. Keesom and H. H. Mooy, 


Proc. Roy. Acad. Amst. 33, 447 (1930), being Comm. 
Leiden No. 209b.] 

20A. Sapper and W. Biltz, Zeits. f. anorg. allgem. 
Chemie 198, 184 (1931). 

21M. C. Neuburger, Zeits. f. anorg. allgem. Chemie 212, 
40 (1933). 


2 W. H. Keesom and K. W. Taconis, Physica 5, 161 . 


(1938); being Comm. Leiden No. 250e. 
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molecule ought to be somewhat larger than the 
crystal value, 2.98A; and we shall adopt the 
value 3.29A, also employed by Kuhn,” which is 
an average of the atomic diameter of Hg as 
determined from the density of liquid mercury 
assuming closest tetrahedral packing (3.25A) and 
from the amalgams™* Ag-Hg (3.09A) and Pb-Hg 
(3.24 + 3.28A). For oxygen, we shall use Vegard’s™ 
x-ray value, 3.48A, for the intermolecular dis- 
tance of two O,2 molecules forming a pair in 
y-oxygen at about 48°K. As to nitric oxide, for 
which no x-ray measurements are available, Ro 
will be set equal to the arithmetic mean, 3.68A,* 
of Vegard’s value for oxygen and the smallest 
distance* between the centers of two neighboring 
Ne molecules at 21°K in a-nitrogen. 

Except for having an upper bound in accord- 
ance with the minimum condition (11), the 
positive quantity p is a disposable parameter. 
The energy of dissociation and the vibrational 
wave number are known spectroscopically for the 
molecules HgA, HgKr and Hg: in the lowest 
vibrational level of the electronic ground state. 
For each of the three molecules the computed 
data come out in good agreement with the two 
spectroscopic values if, and only if, the adjustable 
parameter is taken very nearly equal to 0.28A. 
The same value of p will be adopted for the five 
remaining molecules; for, as regards energies of 
dissociation and vibrational wave numbers, no 
experimental data are available excepting esti- 
mates from nonspectroscopical sources of the 
dissociation energies of (Oz). and (NO)s. 

The dependence of the zero-point energy (15) 
and energy of dissociation (17) on the value 
chosen for the parameter p is shown in Fig. 2, 
which is plotted for HgA but is representative of 
the functions ¢€(p) and Do(p) for the other 
molecules as well. It is seen that the energy of 


3H. Kuhn, Proc. Roy. Soc. London A158, 230 (1937). 
(1928) M. Goldschmidt, Zeits. f. physik. Chemie 133, 397 

*L. Vegard, Nature 136, 720 (1935); Zeits. f. Physik 
98, 1 (1935). Also: W. H. Keesom and K. W. Taconis, 
Physica 3, 141 (1936); being Comm. Leiden No. 240d. 

_*It so happens that 3.68A is just the product of the 
diameter of an NO molecule as deduced from the critical 
temperature and pressure, by the correction factor 1.30. 
For the five inert gases, the departure from 1.30 of the 
ratio of the x-ray value of the distance of two nearest 
neighbors (at or near 0°K) to the (Perit, Terit) atomic di- 
ameter does not exceed 4.2 percent. 

**M. C. Neuburger, Zeits. f. Krist. 93, 1 (1936). (On 

p. 15, replace his reference number 105 by 109.) 
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dissociation is not very sensitive to a change of p 
over a considerable part of its range, whereas the 
zero-point energy (or the vibrational wave num- 
ber) depends much more strongly on the value of 
that parameter.** 


HgA anp HgKr 


The spectra of the molecules HgA and HgKr 
are well known from Oldenberg’s?” measurements 
in fluorescence, and recently the spectrum of 
HgA has also been obtained in absorption*® and 
in an electrical discharge.2® Oldenberg?” sug- 
gested that these bands are caused by transitions 
from some fixed vibrational state of the excited 
molecule Hg’A or Hg’Kr to all the vibrational 
levels of the normal molecule and from this 
interpretation derived the experimental values 
given in Table I for wave number and dissociation 
energy of these two molecules. 

As an example for the actual computation of 
the vibrational wave number and energy of 
dissociation of a normal van der Waals molecule, 
consider the molecule HgA. Refer to the 
corresponding row in Table I. The equilibrium 
distance is taken as Ry =3.39A. Try Rn=3.35A. 
Using the known c’s, the zero-point energy €o 
from (15) and energy of dissociation (17) are 
computed and plotted in Fig. 2 as a function of 
the parameter p, whose upper limit p* as given by 
(11) will also be found in Table I. Choosing 
p=0.28A, the calculated values as listed are 
obtained for the vibrational wave number (14) as 
well as the zero-point and dissociation energies in 
zeroth (15, 17) and first (16, 18) approximations. 
Next the constant A (p) given by (10) is computed 
and the potential energy function V(R) from (9) 
plotted ; it is similar in shape to the V(R) curve 
for the molecule HgHe shown in Fig. 3. Since the 
midpoint of the classical range of oscillations for 
the lowest vibrational level €9 or €; comes at 


** In Fig. 2, R», rather than Ro is kept constant (having 
been chosen such that the midpoint of the level ¢o or «: 
agrees with Ry for values of p not too far off 0.28A). Also 
the approximation of the zero-point energy by (15) [or 
even (16) ] is unreliable for p values so small that the ratio 
€o/A becomes comparable to unity. The curves are, 
— approximate only except in the neighborhood of 


27 (a) O. Oldenberg, Zeits. f. Physik 55, 1 (1929); (b) 

H. Kuhn and O. Oldenberg, Phys. Rev. 41, 72 (1932). 
28H. Kuhn, Proc. Roy. Soc. London A158, 212 (1937). 
29 W. M. Preston, Phys. Rev. 51, 298 (1937). 
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Fic. 2. eo(p) and Do(p) for HgA. (Ro=3.39A; 
Rn =3.35A; p* =0.3998A.) 


about 3.37A, in sufficient agreement with Ro, the 
choice for R,, has been appropriate. The value 
0.28A has been selected for the parameter p in 
order to obtain approximate agreement between 
Oldenberg’s spectroscopic data for the vibra- 
tional wave number and energy of dissociation 
and the values here calculated. 

The procedure followed is analogous for the 
other molecules treated in this paper. As to the 
results for the molecule HgKr, see the corre- 
sponding row in Table I. It also lists the results 
of previous theoretical investigations. 


Mercury vapor is practically monatomic.” 
That two Hg atoms in the normal state (64S) can 
form a van der Waals molecule Hg2(!Z,*) follows, 
however, from the band spectrum to which 
gaseous mercury gives rise. This spectrum con- 
sists of a multitude of continua and more or less 
diffuse bands and has been discussed in detail in a 
recent excellent monograph of Finkelnburg.*! As 
to the spectroscopic values for the energy of 
dissociation and vibrational frequency, reference 
is made to a paper by Kuhn,” who, incidentally, 
found 0.02 percent for the fraction by weight 7 
of Hgs in mercury vapor at a pressure of about 
8 mm Hg and at 448°K. 

The numerical data pertaining to the molecule 
Hg: are assembled in Table I. 


30 A. Kundt and E. Warburg, Ber. 8, 945 (1875); Ann. 
d. Physik [2] 157, 353 (1876). 

31W. Finkelnburg, Kontinuierliche Spekiren (Struktur 
und Eigenschaften der Materie, Vol. 20 (Berlin, 1938)). 


HELLER 


HgHe, HgNe HgXe 


Unlike HgA and HgKr, the spectra of HgHe, 
HgNe and HgXe, observed in fluorescence by 
Oldenberg,?’ and in the case of HgHe in an 
electrical discharge by Preston,?® are continua 
which show no structure (except for HgNe on the 
short wave-length side of the mercury line 
2537). The absence of a diffuse structure on the 
long wave-length side in the HgXe spectrum may 
have been caused by strong scattering?” from the 
quartz bulb used. On Oldenberg’s?”* hypothesis, 
according to which such a structure represents 
the vibrational quanta of the molecule in its 
lowest electronic state, the experimental fact that 
the structure is missing in the HgHe spectrum is 
what should be expected from the present calcu- 
lation inasmuch as, in view of Fig. 3 and the 
numerical data in Table I, the ratio of zero-point 
vibrational energy and depth of the potential 
minimum is so big (€o/A=48.9%; e1/A=42.9%) 
as to suggest that the first excited vibrational level 
would not be stable. But the absence of a 
structure on the long wave-length side in the 
spectrum of HgNe cannot be accounted for in 
this way. 


(O2)2 


The molecule (O2)2 is perhaps the most inter- 
esting of all known van der Waals molecules.® 
That (Oz)2 molecules exist in gaseous oxygen, and 
in the liquid and solid phases as well, follows 
conclusively from the continuous absorption 
bands, some of which are responsible for the pale 
blue color of compressed, liquid and solid oxygen. 
Particulars regarding the absorption spectrum 
caused by (Oz)2 can be found in Finkelnburg’s® 
monograph and in the recent dissertations of 
Guillien* and Herman.* 


% It has no connection with a hypothetical compound, 
Ou, called oxozone, with great chemical activity, which C. 
Harries (Chem.-Zeitung 31, 609 (1907); Liebig’s Ann. 374, 
293 (1910); 390, 235 (1912); Zeits. f. Elektrochemie 17, 
631 (1911); 18, 130 (1912); Ber. 45, 936 (1912)) claimed 
was contained in an appreciable quantity in gaseous and 
liquid ozone. This substance never has been isolated; and 
E. H. Riesenfeld and G. M. Schwab (Ber. 55, 2088 (1922)), 
who succeeded in preparing pure gaseous, liquid and solid 
ozone, have stated that there is no evidence for the existence 
of oxozone. 

%R. Guillien, Ann. de physique [11] 8, 555 (1937); 
being his Thése (Paris, 1937). 

*L. Herman, Ann. de physique [11] 11, 548 (1939); 
being his Thése (Paris, 1939). 
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THEORY OF SOME VAN DER WAALS MOLECULES 


The spectroscopic evidence is corroborated by 
the results of a good many other investigations, 
magnetic and otherwise, concerning the physical 
properties of oxygen. Some of these are listed in 
Table II. 


(NO). 


Table II contains also some information re- 
garding magnetic data which suggest the exist- 
ence of the dimer (NO)> of the ‘‘odd’’ molecule 
NO. There is, in addition, some evidence from 
chemical sources that it is the molecule (NO)> 
that serves as the intermediary product in the 
oxidation of nitric oxide. 


This reaction, of considerable importance at one stage in 
the preparation of nitric acid from nitric oxide and also 
in the theory of the lead chamber process, is of the third 
order and has a small negative temperature coefficient.*® 
As recent experiments of Porret** proved the reaction rate 
constant to be practically independent of the pressure of an 
added foreign gas (nitrogen), Kassel’s*? statistical theory of 
trimolecular mechanisms does not seem applicable to this 
reaction. If, instead, the reaction proceeds in two bimo- 


4 ev 


° 


-6 


Fic. 3. Potential energy and zero-point energy for 
HgHe. (Ro=3. 28A; Rm=3.1A; p=0.28A.) : V(R); 
Morse curve M(R). 


% M. Bodenstein and L. Wachenheim, Zeits. f. Elek- 
trochemie 24, 183 (1918) [273+363°K]; M. Bodenstein 
and Frl. Lindner, Zeits. f. physik. Chemie 100, 87 (1922) 
Hass oak: F. Matthes, Dissertation (Berlin, 1933) 
143+273°K J; E. Briner, W. Pfeiffer and G. Malet, J. 
chim. phys. 21, 25 (1924) [90 +348°K]. 

D. Porret, ’Dissertation (Berlin, 1937), 1-6. 
“1. &. Kassel, J. Phys. Chem. 34, 1777 (1930); also: 
Homogeneous Gas Reactions (New York, 1932), 
ection 
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lecular steps, the intermediary substance may be either 
NO; or (NO)s. The existence of NO; is assured both 
spectroscopically** and on the basis of chemical kinetics,*® 
though it has not yet been isolated;@ but the complete 
absence*! of its characteristic lines in the absorption 
spectrum of a mixture of nitric oxide and oxygen would 
appear to rule it out as the intermediate in the reaction 
under consideration (unless one postulate‘? ad hoc the 
existence of an isomer of NO3;). The belief, now held by 
Bodenstein,** that the intermediate is the dimer (NO)s is 
strongly supported*** » by the close analogy and parallelism 
in the chemical behavior of the reactions of nitric oxide 
with oxygen, bromine*** and chlorine.**> 


From the temperature dependence of the reac- 
tion rate constant as deduced from Bodenstein’s 
older measurements of the rate of the oxidation 
of nitric oxide, the energy of dissociation of an 
(NO)2 molecule has been estimated by Kassel*’ to 
be 56(10)-* ev and by Christiansen*® to be 
35(10)-* ev. This estimate is several times larger 
than the value of about 8 millivolts} computed 
by the present author, suggesting that first-order 
attractive forces arising from the tendency for 
saturation of the spin moments of two normal 
(711) NO molecules may be present, as was 
already pointed out by Briegleb“ in an attempt 
to explain the discrepancy between the experi- 
mental and computed second virial coefficient of 
nitric oxide (at 273°K). (NO)s would not then be 


a pure van der Waals molecule. It is interesting 


38 P. E. M. Berthelot, Comptes rendus 92, 82 (1881); the 
references in footnotes 2 and 3 of Sprenger, reference 39; 
also P. Hautefeuille and J. Chappuis, Comptes rendus 92, 
946 (1882), and J. Chappuis, Ann. Ec. Norm. Sup. [2], lL, 
137 (1882); E. J. Jones and O. R. Wulf, J. Chem. Phys. 
5, 873 (1937). 

*°G. Sprenger, Zeits. f. Elektrochemie 37, 674 (1931), 
and the papers cited in his footnote 1. (Decomposition of 
ozone in the presence of nitrogen pentoxide.) See also: G. 
M. Calhoun and R. H. Crist, J. Chem. Phys. 5, 301 (1937); 
7, 719 (1939); and G. v. Elbe and B. Lewis, J. Am. Chem. 
Soc. 59, 2022 (1937); 61, 1350 (1939). 

40 A. Klemenc and W. Neumann, Zeits. f. anorg. allgem. 
Chemie 232, 216 (1937). 

“aH. J. Schumacher and G. Sprenger, Zeits. f. angew. 
Chemie 42, 697 (1929). 

“2FP, Swarts, 5¢ Conseil de Chimie, Instituts Solvay (Paris, 
1935), p. 37. 

43M. Bodenstein, (a) reference 42, P. 1 ff.; (b) Helv. 
Chim. Acta 18, 743 (1935); (c) Zeits. f . physik. Chemie 
A175, 294 (1935). 

4 (a) W. Krauss, Zeits. f. Dig os Chemie A175, 295 
(1935); (b) W. Krauss and Saracini, ibid. A178, 245 
(1936). The announced simultaneous discussion of the three 
reactions of nitric oxide has not been <n as yet. 

4 J. A. Christiansen, reference 42, p. 31. 

+ Just as for the other molecules, p=0. 28A has been as- 
sumed. No matter what particular value is taken for the 
parameter p, however, the calculated energy of dissociation 
would never exceed 13 millivolts. 

 G. Briegleb, Zeits. f. physik. Chemie B23, 127 (1933). 
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TABLE II. Some experimental evidence for the existence of (Oz)2 and (NO)s. (r denotes the fraction by weight of the dimer 
in the undiluted substance.) 


PHENOMENON PHASE AUTHOR REMARKS 
Oxygen 
1. ae ay from reduced equation Gas E. C. Wiersma and C. J. Gorter* 
of state 
2. Anomalous Faraday effect Gas R. Serber” 
H. R. Woltjer, C. W. Coppoolse 
3. Deviations from Curie’s law of Gee and E. C. Wiersma‘; 
paramagnetism E. C. Wiersma, W. J. de Haas and 
W. H. Capel 
at 77°K; r(extrap.) 
=0.03% at 293°K 
and 1 atmos. 
4. Cotton-Mouton effect Liquid P. Lainé’® 
5. Kerr effect © Liquid R. Guillien fa p. 619/20, 643] , 
6. Concentration dependence of Liquid R. Guillien [%3, p. 578] 7=40% at 77°K 
extinction coefficient 
7. Heat of mixing of liquid oxygen Liquid P. Lainé® D=10 mv 
and nitrogen 
8. Structure of y-oxygen Solid L. Vegard*5 Ro=3.48A 
9. Magnetic susceptibility Solid A. Perrier and H. Kamerlingh Onnesi 
Nitric Oxide 
1. Faraday effect Lasid } H. Bizette and B. Tsaii 7=93% at 110°K 
Liquid H. Bizette and B. Tsai* 7r=97% at 110°K 
2. Magnetic susceptibility Solid E. Lips! 


: Physica 12, 316 (1932); poee Comm. Leiden, Suppl. No. 73c. 
b Phys. Rev. 41, 489 (1932). 
on _— Roy. Acad. Amst. 32, 1329 (1929); being Comm. Leiden No. 
4 Proc. Roy. Acad. Amst. 34, 494 (1931); being Comm. Leiden No. 


ej. Am. Chem. Soc. 46, 2027 (1924). 
t J. de physique [7] 6, 14S, 62S (1935). 


to note that Wohl’s’ second virial coefficient of 
oxygen (from 121 to 373°K) is in excellent 
agreement with experiment, and that the experi- 
mental and theoretical values in Table I for the 
energy of dissociation of the molecule (O2)2 are 
fairly concordant, although no first-order forces 
had been assumed for the attraction between 
two normal (*2) Oz molecules. 


DISCUSSION 


The main numerical results have been as- 
sembled in Table I. The three last columns, 
giving the vibrational wave number, zero-point 
energy, and energy of dissociation, are subdi- 
vided, showing the results of the present investi- 
gation on the left-hand side and the data obtained 
by other authors to the right. Experimental 
results are indicated by figures printed in italics. 
The results derived from a potential energy 
function, on the other hand, are differentiated, 
except for the wave numbers, by means of the 


& Biréfringence magnétique (Paris, 1936), p. 35, 44. 
b Physica 3, 895 (1936); being Comm. Cate ‘No. 242d. 
— Roy. Acad. Amst. 16, 894 (1913); being Comm. Leiden No. 
1 
i Comptes rendus 204, 1638, 1965 (1937). 
k Comptes rendus 206, 1288 (1938). 
! Helv. Phys. Acta 8, 247, 380 (1935). 


symbols (€) or (Do), (€1) or (D1), or (A), ac- 
cording as they were computed with the aid of 
the oscillator zero-point energy ¢€o, the Morse 
zero-point energy €1, or are depths of the potential 
at the minimum. 

It is believed that the repulsive law Ae-*’’, 
successfully employed by Born and Mayer? for 
ionic crystals, has proved satisfactory for the van 
der Waals molecules here under discussion. 

If the term in R-™ in the attractive part of the 
potential (9, 10) were replaced by zero, the error 
thus introduced would be only 2.1 percent or less 
in the depth A of the potential minimum and less 
than 5.3 percent in the zero-point energy €o (or €1) 
or wave number jo for the eight molecules ‘con- 
sidered (assuming p=0.28A). However, the con- 
tribution of the terms in R-*® and R-, which 
were first considered by Margenau“’ and in which 
the dipole-octupole interaction has been included 
in this paper, is appreciable at R= Rn. 


47H. Margenau, Phys. Rev. 38, 747 (1931). 
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It is of interest to observe how a change in the 
polarizability of the constituents is counteracted 
by a corresponding change in size and especially 
in reduced mass of the van der Waals molecule 
such that the vibrational wave number 7» is 
nearly the same for the eight molecules, varying 
from 28.1 cm~! (HgXe) to 46.6 cm~! (HgHe). A 
similar result has already been obtained by 
Salow® for the molecules OoZ (Z=Os, He, Ne, A, 
Ne, CO2), his ranging from 39 cm—! (O2Ne) to 
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95 cm~! (O2(CO»)). The spectroscopic values for 
normal-state diatomic molecules with a chemical 
bond, on the other hand, extend all the way from 
about 42 cm~! up to about 4400 
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The Entropy of Formic Acid. The Heat Capacity from 15 to 300°K. Heats of Fusion 
and Vaporization 
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The heat capacity of solid and liquid formic acid has been measured from 15 to 300°K. 
The melting point is 281.40°K (0°C =273.10°K). The heat of fusion is 3031 cal. mole~! and 
the heat of vaporization at 298.10°K is 4754 cal. mole. (1 mole = 46.0260 grams.) The vapor 
pressure at 298.10°K is 4.31 cm of mercury. From the third law of thermodynamics and the 
calorimetric data the entropy has been calculated. To the value obtained from the calorimetric 
measurements an entropy of R/2 In 2, due to the random orientation of hydrogen bonds in 
the solid which has been predicted by Pauling, has been added. This gives 31.51 cal. deg.~! 
mole for the entropy of the liquid at 298.10°K and 47.46 for the entropy of 46.0260 grams 
of gas, P=4.31 cm, in equilibrium with the liquid at 298.10°K. These are the values that 


should. be used in thermodynamic calculations. 


HIS paper presents the results of a low 
temperature calorimetric investigation of 
formic acid, leading to values for the entropy of 
the liquid and the gas. The value of the entropy 
obtained is compared with those deduced from 
free energy and heat changes. 

The measurements were made in a calorimeter 
of laboratory designation Gold Calorimeter IV. 
The resistance thermometer was calibrated 
during the measurements by comparison with 
standard copper-constantan thermocouple W-21. 
The calorimeter and the calibration of the 
thermocouple have been described by Kemp and 
Giauque.! The apparatus and methods of 
measurement and calculation have been de- 


* Fellow of the Lalor Foundation, pei Now at the 
Massachusetts Institute of Technolog 
' Kemp and Giauque, J. Am. rg ‘Soc. 59, 79 (1937). 


scribed in previous papers? from this laboratory. 
The ice point was taken as 273.10°K, and the 
defined calorie, 4.1833 international joules, was 
used throughout. 


PREPARATION AND PuRITY OF Formic AcIpD 


Since the commercially available formic acid 
contained an objectionable amount of acetic 
acid, which might have proved difficult to 
remove, it was decided to prepare the sample by 
the reaction of a formate with a mineral acid. 
Phosphoric acid was chosen since it has a 
negligible vapor pressure and preliminary experi- 
ments had shown that it does not rapidly 


2 (a) Giauque and Wiebe, J. Am. Chem. Soc. 101 
(1928). (b) Giauque and Johnston, J. Am. Chem. ond 51, 
(c) Giauque and Egan, J. Chem. Phys. 5, 
45 (193 
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TABLE I. Melting point of formic acid 0°C =273.10°K. 


RESISTANCE 
TIME % MELTED THERMOMETER —THERMOCOUPLE 
4-25-39 
9:26 p.M. Heated into melting point 
11:25 P.M. 8 281.365 281.36 
4-26-39 
12:31 a.m. . 8 281.360 281.36 
10:36 A.M. 8 281.358 281.35 
11:37 a.m. Supplied heat 
12:22 P.M. 19 281.385 281.38 
2:04 P.M. 19 281.386 281.38 
3:21 P.M. 19 281.387 281.38 
5:28 p.m. Supplied heat 
6:01 P.M. 56 281.392 281.39 
8:15 P.M. 56 281.387 281.39 
10:10 P.M. 56 ~ 281.391 281.39 


11:16 P.M. 56 281.388 281.39 
Accepted value 281.40+.05 


decompose formic acid as does concentrated 


sulfuric acid. 85 percent phosphoric acid was 


- heated to 250°C to drive off some of the water 


and then cooled to room temperature. This acid 
was then added to ammonium formate and the 
formic acid was distilled off. During the distilla- 
tion the temperature of the flask was kept below 


50°C. The formic acid so prepared contained an 
appreciable amount of water which was satis- 
factorily removed by distilling the sample into a 
bulb containing anhydrous calcium sulfate 
(Drierite) which had been pumped under a high 
vacuum at 250°C. No decomposition of the 
formic acid such as was observed in preliminary 
experiments with phosphorus pentoxide or 
concentrated sulfuric acid was found with the 
anhydrous calcium sulfate. The sample was then 
distilled at room temperature through a small 
fractionating column and the middle fraction 


- selected for the measurements. 


The preparation and purification of the formic 
acid was performed in an all-glass system from 
which air had been evacuated. All stopcocks 
were lubricated with Apiezon grease. We ob- 
served no evidence of reaction of the Apiezon 
grease with formic acid such as has been reported 
by Ramsperger and Porter* for ordinary stopcock 
grease. To avoid the slow decomposition of 
liquid formic acid to which Coolidge! has called 
attention, the sample was kept constantly frozen 


} TABLE II. The heat capacity of formic acid. 0°C=273.10°K, molecular weight 46.0260. 


APPROX. Cp APPROX. Cp 
AT CAL. DEG.~! MOLE™! AT CAL. DEG.~! MOLE! 
15.05 0.479 163.49 11.02 
17.27 ae 0.673 170.63 7.0 11.23 
19.61 2.5 0.936 177.74 7.3 11.46 
22.28 2.7 1.235 184.94 7.0 11.69 
25.15 3.0 1.540 192.12 7.3 11.88 
28.28 3.2 1.908 199.30 7.0 12.10 
31.43 3.1 2.308 206.21 6.8 12.33 
34.63 3.3 2.740 213.07 7.1 12.57 
38.31 4.0 3.196 220.31 6.8 12.79 
42.59 4.6 3.670 227.85 7A 13.06 
47.39 5.1 4.163 234.84 6.8 13.31 
52.64 5.4 4.707 241.89 6.6 13.53 
57.98 Le §.212 248.78 6.4 13.70 
63.16 5.2 5.652 255.57 6.2 13.94 
68.56 5.6 6.114 262.12 5.9 14.17 
74.32 6.0 6.592 268.56 5.7 14.40 
80.47 6.3 7.071 274.92 5.4 15.24* 
86.88 6.6 7.522 281.40 Melting point 
93.82 te 7.963 275.28 4.2 23.39 
100.79 6.7 8.368 279.68 4.1 23.49 
107,72 8.747 284.16 4.0 23.55 
114.83 7.0 9.099 286.42 4.0 23.53 
121.63 6.6 9.422 286.46 4.0 23.63 - 
128.51 7.4 9.701 289.32 44 23.62 
135.68 7.0 10.02 290.76 4.0 23.63 
141.99 5.5 10.25 294.90 3.6 23.66 
148.09 6.6 10.47 299.79 3.5 23.68 
155.65 8.4 10.74 


* High because of premelting. 


3’ Ramsperger and Porter, J. Am. Chem. Soc. 48, 1267 (1926). 
sy J. Am. Chem. Soc. 50, 2166 (1928). 
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Fic. 1. Heat capacity in calories per degree per mole of formic acid. 


except during transfer and during the measure- 
ments above the melting point. 

The purity of the formic acid was estimated 
from a series of short heat capacity measure- 
ments immediately below- the melting point. 
The amount of liquid-soluble solid-insoluble 
impurity was six parts per 100,000 on a molal 
basis. The same value may be obtained from 
the change of the melting point with fraction 
melted. 


THE MELTING POINT 


A summary of the melting point observations 
is given in Table I. The temperatures measured 
by the resistance thermometer are given to 
0.001° only because of the high relative accuracy. 
In calculating the accepted value of the melting 
point, correction has been made for the lowering 


of the melting point because of the impurity 
present. Of previous determinations of the 
melting point the most recent is that of Coolidge’ 
who obtained the value 8.25°C as compared to 
our value of 8.30°C (281.40°K). 


THe Heat Capacity oF Formic Acip 


The amount of formic acid present during the 
measurements of the heat capacities and the 
heats of fusion was determined by distilling the 
sample into an evacuated glass bulb which was 
weighed. All weights have been corrected to 
vacuum. 156.041 g or 3.3903 moles of formic acid 
were used for the measurements. The heat 
capacity observations were corrected for vapori- 
zation into the small volume of the calorimeter 


5 Coolidge, J. Am. Chem. Soc. 52, 1874 (1930). 


166 J. W. STOUT AND L. H. FISHER 


TABLE III. The heat capacity of formic acid. Values from a smooth curve through the observations. 
0°C=273.10°K, molecular weight 46.0260. 


DEVIATION DEVIATION 
Cp G. L. AND P. Cp G. L. AND P. 
CAL. DEG.~! MOLE™! —THIS REs. % CAL. DEG.~! MOLE! —THIs REs. % 
15 0.474 170 11.21 
20 .980 180 11.51 +91 
25 1.527 190 11.82 +10.6 
30 2.124 200 12.13 +11.5 
35 2.781 210 12.46 
40 3.384 220 12.81 
45 3.921 230 13.14 
50 4.439 240 13.47 +14.5 
60 5.381 250 13.78 +15.5 
70 6.241 +13.7 260 14.10 +16.2 
80 7.033 +10.0 270 14.43 +16.8 
90 alee +6.7 280 14.75 
100 8.325 281.40 Melting point 
110 8.862 Liquid 
120 9.338 275 23.38 
130 9.779 280 23.42 
140 10.17 290 23.61 —0.5 
150 10.55 300 23.68 
160 10.89 
TABLE IV. Heat of fusion of formic acid. Melting point 281.40°K, molecular weight 46.0260. 
Heat INPut CpdT PREMELTING AH Fusion 
Run Temp. RANGE, °K CAL. MOLE~! CAL. MOLE™! CAL. MOLE! CAL. MOLE“! 
1 279.50-286.67 3214.3 186.4 4.7 3032.6 
2 278.11-283.71 3156.3 127.8 2a 3031.2 
3 278.36-284.04 3159.3 132.6 2.9 3029.6 
This research, mean 3031 +3 
Berthelot! 2524 
Coolidge? 3035 


1 Berthelot, Ann. Chim. Phys. [5] 4, 74 (1875). 
2 Cf. reference 5 of text. 


and line. In order to insure thermal equilibrium 
about 2 cm pressure of helium was admitted to 
the calorimeter at 250°K. No evidence of poor 
thermal conductivity of the solid, such as has 
been mentioned by Gibson, Latimer and Parks® 
was observed. Except for the runs immediately 
below the melting point, thermal equilibrium 
in the calorimeter was attained within five 
minutes after the addition of energy. 

The heat capacity measurements on solid and 
liquid formic acid are presented in Table II and 
represented graphically in Fig. 1. Two points 
on the heat capacity of the supercooled liquid 
are included. 

Values read from a smooth curve through the 
observations are given in Table III. These 
values are believed accurate to 0.2 percent above 
35°K, to 1 percent at 20°K and to 5 percent at 


6 Gibson, Latimer and Parks, J. Am. Chem. Soc. 42, 
1533 (1920). 


15°K. In Table III the earlier measurements of 
Gibson, Latimer and Parks* are compared with 
our results. 


THE HEAT OF FUSION 


The results of the measurements of the heat 
of fusion are given in Table IV. The runs were 
started several degrees below the melting point 
and ended at a temperature above. Corrections 
were made for the heat capacity of the solid 
before melting and the heat capacity of the 
liquid afterwards. Small corrections were also 
applied for the amount of premelting in the solid 
at the start of the run. For comparison the 
values of the heat of fusion obtained by previous 
investigators are also listed in Table IV. 


THE HEAT OF VAPORIZATION 


Since the vapor pressure of formic acid at 
298°K is only some 4.3 cm, it was not possible 


to 
pc 
pr 
m. 
co 
va 
ev 
ba 
gle 
th 
| joi 
thi 
col 
ap 
| aci 
rin 
tio 
cal 
wit 
ten 
bet 
(lal 
val 
| ure 
4¢ 
vat 
the 
S 
| use 
| atte 
dou 
the 
44 
| capi 
4 ar 
capi 
vapi 
heat 
| sligt 
diffe 
| men 
| the 
| obta 
| the 
| com 
D 
| after 


ENTROPY OF 


to collect the sample vaporized in an absorption 
bulb open to the atmosphere as is ordinarily 
done with a gas vaporizing at its normal boiling 
point. The procedure adopted was to utilize the 
pressure drop through a glass capillary to 
maintain the pressure in the calorimeter at that 
corresponding to the desired temperature of 
vaporization. The formic acid was condensed in 
evacuated glass bulbs which were cooled in a 
bath of carbon dioxide and methyl alcohol. The 
glass bulb, which was attached to the line from 
the calorimeter by means of a ground glass 
joint, was removed and weighed before and after 
the heat of vaporization. All weights have been 
corrected to vacuum. Corrections have been 
applied for the difference in the amount of formic 
acid present in the gas space above the calo- 
rimeter before and after the heats of vaporiza- 
tion. The electrical energy introduced into the 
calorimeter was corrected for the heat exchanged 
with the surroundings and for the change in 
temperature of the calorimeter and contents 
between the start and end of the energy input 
(labeled {C.aidT in Table V). In correcting the 
values of the heat of vaporization at the meas- 
ured average temperatures to 298.10°K the value 
4 cal. mole deg.—! for the change in heat of 
vaporization with temperature, derived from 
the measurements of Coolidge® was used. 

Several different rates of evaporation were 
used in order to determine if slowness in the 
attainment of equilibrium between single and 
double molecules in the gas appreciably affected 
the value obtained. Capillary A, 0.1 cm I.D. and 
44 cm long was used for runs 1, 2 and 3 and 
capillary B, 0.1 cm I.D. and 69 cm long for runs 
4 and 5. Runs 6 and 7 were taken using the two 
capillaries in parallel. The data on the heat of 
vaporization are presented in Table V. The 
heats of vaporization at the slowest rate are 
slightly larger than the others but since the 
difference is of the order of the expected experi- 
mental error we have adopted the mean of all 
the measurements as the best value. The value 
obtained from a smooth curve drawn through 
the observations of Coolidge® is included for 
comparison in Table V. 

During the equilibrium periods before and 
after each heat of vaporization the vapor 
pressure of the liquid formic acid in the calo- 
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rimeter was measured. The readings were made 
with a cathetometer used as a comparison 
instrument with a standard meter bar mounted 
in a case with the mercury manometer. The 
pressure readings were corrected to international 
cm of mercury. The average value of the pressure 
of liquid formic acid was 4.31 cm at 298.10°. 
This is to be compared with the value 4.29 cm 
obtained from a smooth curve through the 
measurements of Coolidge’ and that of 4.16 in 
the International Critical Tables.’ Coolidge* has 
called attention to the marked lowering in the 
vapor pressure of formic acid produced by small 
amounts of water. 


THE ENTROPY FROM CALORIMETRIC DATA 


The entropy was calculated by graphical 
integration of fj7C,d In T plus the entropies of 
fusion and vaporization. A summary of the 
entropy calculation is given in Table VI. 

Gibson, Latimer and Parks® obtained 24.30 
cal. deg.—! mole“ for the entropy of the solid at 
the melting point compared to our value of 18.68. 
About 4 cal. deg.-! mole of the discrepancy is 
caused by the extrapolation from the lowest 
measured point of Gibson, Latimer and Parks 
at 70.8°K to the absolute zero. 

Pauling® has pointed out that the entropy 
obtained from low temperature calorimetric 
measurements on formic acid should be too low 
by 4RIn2 per mole of monomer caused by a 
random orientation of the two hydrogen bonds 
present in a dimer molecule in the solid. A 
similar random orientation of hydrogen bonds 
has been confirmed experimentally in the case of 
H,0° and D,O.'" If the predicted random 
orientation in formic acid occurs, as seems very 
likely, an entropy of 0.69 cal. deg.-' mole 
should be added to the values listed in Table VI 
for the liquid and the gas. 

The entropy of the gas listed in Table VI is 
the sum of the entropies of the monomer and 
the dimer present in 46.0260 g of gas under its 
saturation vapor pressure, 4.31 cm, at 298.10°K, 
the entropy of each species being calculated at 
its partial pressure. Careful measurements of 
the equilibrium between the monomer and the 

TInt. Crit. Tab. Vol. III, p. 215. 

§ Pauling, J. Am. Chem. Soc. 57, 2680 (1935). 


® Giauque and Stout, J. Am. Chem. Soc. 58, 1144 (1936). 
10 Long and Kemp, J. Am. Chem. Soc. 58, 1829 (1936). 
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TABLE V. Heat of vaporization of formic acid, molecular weight 46.0260. 
TIME OF ELEc. HEAT 
MOLES ENERGY ENERGY EXCHANGE J 428.10 
RuN T,°K Evap. INPUT, MIN. MOLE™! CAL. MOLE! CAL. MOLE! CAL. MOLE™! CAL. MOLE! 
1 297.99 0.21482 56 4817.9 4.1 —75.9 4746.1 4747 
2 298.17 -21801 56 4737.1 3.9 +12.8 4753.8 4754 
3 298.15 .22031 56 4685.9 2.4 67.2 4755.5 4755 
4 298.10 21945 94 4732.1 14.8 15.9 4762.8 4763 
5 298.15 .22181 94 4678.6 15.6 66.3 4760.5 4760 
6 297.74 .22229 42 4610.0 21.2 108.9 4746.1 4748 
7 297.69 .21715 42 4713.0 35.8 0.0 4748.8 4750 
Mean 475445 
Coolidge’ 4805 ] 
ar 
dimer in the gas have been reported by Ramps- HCOOH (/) = H2(g) +CO2(g) (2) ty 
berger and Porter? and by Coolidge.* These ; in 
measurements do not, however, permit an by Bredig, Carter and Enderli.” Parks and 8 
estimation of AS of the dissociation reaction Huffman" employing a heat of combustion of ail 
more closely than about 1 cal. deg.-! mole-!. A __ liquid formic acid based on the measurements of 
more accurate treatment of the dissociation Berthelot and Matignon™ and Jahn’*® have ps 
equilibrium and the answer to the question of calculated AS%293.10=30.3 cal. deg.-' mole for i 
random orientation of hydrogen bonds in solid reaction (1) and 47.9 cal. deg. mole for (2). 
formic acid must be postponed until values of Using the S%bos.10 values, H2O(/) = 16.75,'* CO(g) 
the vibrational frequencies of at least the gaseous 9=47.32,"" H2(g) =31.23,!% and CO2(g) =51.07," 
monomer molecule have been obtained. one obtains for S%93.19 of liquid formic acid be 
The equilibrium 33.8 cal. deg.-! mole! from reaction (1) and | 
34.4 from reaction (2). The lack of agreement - 
OOK) (1) with the value 30.82+0.69=31.51 obtained 
has been studied by Branch" and the equilibrium from our measurements combined with the at 
assumption of random orientation of hydrogen 
molecular seid 46.0260. bonds in the solid is probably due to inaccuracy f 4 
in the heats of combustion. - 
0-15°K, Debye extrapolation 0.164 nee 
Fusion 308128140 rouBredig, Carter and Enderli, Monats. F. Chem. 54, F(e 
Parks and Huffman, The Free Energies of Some Organic 
ompounds emical Catalog Co., 
+0.1 aso and Matignon, Comptes rendus 114, 1145 or 
S of 46.0260 g of actual gas (P=4.31 15 h Wied. A 37, 408 (1889) 
cm) in equilibrium with the liquid 16 tossini elie ct and Jonston, J. Research Nat. Bur. 
at 298.10°K a cal. deg. 1 Stand. 24, 369 (1940). ’ 
17 Clayton and Giauque, J. Am. Chem. Soc. 54, 2610 
(1932); 55, 5071 (4933). SO 1 
Leremicwere rd 18 Giauque, J. Am. Chem. Soc. 52, 4816 (1930). 
11 Branch, J. Am. Chem. Soc. 37, 2316 (1915). 19 Giauque and Egan, J. Chem. Phys. 5, 45 (1937). 
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A symmetrical method of finding the Thiele semi-invariants \, in the partition function 
for a binary solid solution is described. In addition to the first four \’s already known, As and Ags 


are calculated. 


N studying the equilibrium properties of a 
binary solid solution consisting of N4A atoms 
and NzB atoms and inhabiting a lattice of the 
type AB (a lattice which can be divided into two 
interpenetrating sublattices, say a sublattice and 
8 sublattice), we are faced with the construction 
of the partition function* F(6, T) 


xX 
(Vaat Ven—2Vaz)X 


Xexp — (1) 
kT 


where Vas, Vee and Vaz are the interactions 
between AA, BB and AB pairs of nearest 
neighbors, and f(N,3N0@, Na—3N6@, X) is the 
number of ways of arranging N, particles on a 
lattice of the type AB consisting of N sites in 
total, with }N@ particles on the a@ sublattice, 
particles on the sublattice, and 
further with the particles forming X pairs of 
nearest neighbors or first neighbors. At T= ~, 
F(@, T) is simply 


f(N, Ns—3N0, X) 
xX 


so we can put log F(@, T) as 


or 


3N a 
log )( ) (2) 
2! 


3N6 3! 
where a stands for — and 
for By an ingenious method 
Kirkwood has been able to obtain i, \2 and 
A3s.' He gave for the special case Na=3N, 


* The interactions between more distant neighbors and 
the factor in the partition function caused by vibrations 
and excitations are left out for simplicity. 

Kirkwood, J. Chem. Phys. 6, 70 (1938). 


X(1—6)?(1—26)?, where z is a coordination 
number of the lattice denoting the number of 
first neighbors of a site. A subsequent paper by 
Bethe and Kirkwood? gave 


where Zaa’ is another coordination number de- 
noting the number of first neighbors common to 
the sites a and a’. 

In a paper by the writer,’ where Aj, A2 and A; 
are given for the general case Na+}N, it is 
unfortunately remarked that \, and the higher 
\’s cannot be found by continuing Kirkwood’s 
method. This is due to the fact that the intro- 
duction of coordination numbers other than z 
was then not considered. In the present note, 
Xa, As and Ag will be given for the general case. An 
attempt to calculate them is worth while, for, at 
the critical point of the setting in of the long 
distance order, the higher moments such as As 
and \¢ do not seem to converge quickly, as can be 
seen below. 

Since 


a? 
x 2! 3! 


we can write 


F(6, T)=(X f) 
where 
Mi=(X), M2=(X*), Ms=(X), 


2 Bethe and Kirkwood, J. Chem. Phys. 7, 578 (1939). 
169 


3 Chang, Proc. Roy. Soc. 173, 48 (1939). 
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with (u) denoting the a priori average of u over 
the }:xf ways of arrangements. In the first 
paper by Kirkwood, X is written as Yoavata, 
where vy, denotes the number of A atoms in the 
site a of the a sublattice and is accordingly either 
zero or unity, and m, denotes the number of A 
atoms among the first neighbors of the site a, the 
summation being taken over all a’s. The problem 
therefore reduces to the calculation of ((}>vama)"). 
In the paper by Bethe and Kirkwood, X is 
written as 


(5) 


where £2 denotes the number of A atoms in the 
site a of the a sublattice, m the number of A 
atoms in the site 6 of the 6 sublattice, and Xs is 
an element of a neighbor matrix which is unity 
if sites a and 6 are first neighbors and zero 
otherwise. The authors have not shown how to 
calculate but simply referred 
to Kirkwood’s paper for the calculation of 


((Xvata)"). Needless to say, the values found for 
and are necessarily the 
same. However, we shall describe below how to 
calculate ((SoAatam)"), for two reasons. First, 
such a calculation enables us to keep the sym- 
metry between the two sublattices always in 
sight, and second, it can be easily adapted (or 
generalized) to include interactions between 
more distant neighbors,* or to cases with other 
types of superlattice formations. For example, if 
we study the superlattice formations of the type 
AB; with the help of neighbor matrices similar 
to Xs in nature and restrict ourselves to the task 
of finding and we find that this can 
easily be done, the labor required being much less 
than that required in applying Bethe’s method.‘ 
As a latent heat phenomena instead of a dis- 
continuity of the specific heat is the main feature 
there and the amount of the latent heat does not 
depend so sensitively on dz, As, etc., we shall 
omit here discussions on superlattice formations 
of the type ABs. 


Consider ((>-Aafanv)?) for definiteness. A binomial expansion gives 


where a and a’, as well as b and b’, denote different sites. We first replace (fa: + +m: ++) by (fa: + *)(m***), 
since the a priori average of a product of ¢’s and 7’s is the product of the a priori averages of ¢’s and 
n's. Next, we replace by 61, (fata) by 62, (m?)=(m) by 61’, (non’) by 62’, since is the 
same for different a sites, ({afa’) is the same for whatever pairs of a and a’ sites, etc. After the replace- 
ment, (6) becomes 
2 

a,b a,a’,b a, b, b’ a, a’, b, b’ 
Hence we need to find 6;, 62, 61’, 62’, and (or > Aas), Each 
of the last four being a summation over all the suffices that appear in it. In finding ((S[Aafan)")’ 
similar replacements can be done, and we are left with quantities of two kinds. Quantities of the first 
kind are 61, 02, and 61’, 02’, where is *fa(n—1)), a, a’, being any different sites 
among the a sublattice. Quantities of the second kind are summations of the products of a succession 


of 
The different 6’s can be found as follows. In the first place, 6; is nothing but 6. From the relation 


we get 


. — are is provided by the subsequent paper, which includes second-neighbor interactions in the calculation of 
1, Ag and A3. 
* Peierls, Proc. Roy. Soc. 154, 207 (1936). 
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From the relation 


1)(;N—2)6s, 
(7) 
In a similar way we can show 6,;'=6’, ---, and 


we get 


In general 


The different summations of the product of two \.s’s can be found as follows. In the first place we 
have from the definition of \. and 2, 


and hence }'asAa=3Nz. Now we have the relation 
(Aad Aare) —D Daw, 
a,a’,b a,b a* a,b 


where in the summation over a*, a* may correspond to any site in the a sublattice, not necessarily 
different from the site a. Hence 


a,a’,b 
This is also the value of S°a,o,»'AaAav due to the symmetry between the a and 8 sublattices. Further we 
have the relation 


a, a’, b,b’ a,b, b’ a* a,b, b’ 
from which we get 
a,a’,b,b’ 


In a similar way, we can get the summation of the products of three \.s’s. Because of the symmetry 
between the a’s and the 0’s, there are only four significantly different products, i.e., }>AabAavAav’’, 
DL ANd They are given byt 


DX =D, (AadAad? Aa») — AadAaw = 3N2(sz—1)?, 
X =D (AabAad? Liat — = 
+N2(52*—62+2). (11) 
By using the same method we can express the summations of the products of four or five \q’s in 


terms of z and a quantity y defined as (2/N)>-\aAavAa'tAa'w Which may be considered as a new 
coordination number. Continuing this method further, we express the summations of products of six 


*It is understood from now on that all summations without subscripts under the summation sign are summations 
— all the subscripts in the summand, and summations with subscripts under the summation sign are taken over these 
subscripts only. 

} For clarity, we repeat that in the summation over a* (or o*), a* (or b*) may correspond to any site in the @ (or 8) 
sublattice, not necessarily different from a or a’ or a” (or b, b’, ---), while in summations over a, a’, «++ (or }, b’, ---), 
“the a’s (or b’s) are necessarily different sites. 
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or seven \qs’s in terms of z, y and two quantities y: and 72, defined as (2/.N) S-avAav-Aad’AabAa’b’Aa’d’’ 
and (2/N) S-Yasdav/Aa'b’Aa’b’Aa''bAa’s>- The writing out of expressions for all such summations is 
cumbersome and will be omitted. Evidently we can continue the method still further to find summa- 
tions of products of eight or more \.»’s with the introduction of other coordination numbers of the 
same nature as y, 71 and 72, but this will not be done here. 

With the 6’s and the summations of products of three, four, five, or six \.,’s completely determined, 
we obtain M,, M2, --- up to Ms. Expanding the M’s in powers of N-', substituting the expanded 
series into (4), and comparing the F(6, T) so obtained with (2), we find 
\1=43N200’, As 
—6)(1—0’)(1— 60+ 66?) (1 —66’ +66’) —6)2(1—8’)?, 

(1 —0)(1—6’) (1 —26)(1 —26’)(1 —126+126?) (1 —120’ +120’) 
— 6)2(1 
Ae = 3.200’ (1 — 6) (1 — 6’) {1 —300(1 — 6) (1 — 20)? — 300’(1 — 6’) (1 —20’)?} 
+15 N(13y+172)676’?(1 —6)?(1— 0’)? +30N(71—30y — 282) 670’2(1 — 6)?(1 — 0’)? 
X {0(1— 0) — 6’) } +60N(-y2—471+66y +422) — 6)9(1— 


(12) 


where terms of the order smaller than N are not written out. It is seen that the \’s are linear in the 
coordination numbers 2, y, yi and y2, meaning that they do not contain powers of y, 2, y1 and ye, nor 
their products. Whether this holds for the higher \’s remains to be seen, but it seems very plausible. 

The connection between the coordination numbers y, y: and yz and the number 2,. of first neighbors 
common to the sites a and a’ is as follows. Since aa is nothing but }>»AaAa», we have 


whereupon 
= FN {De —2(8—1)}, 


y= —2(2—1). (13) 


or 


In a similar way we can prove 


Da 3(2—2) Le Bear +22(2—1)(z—2). 


These equations furnish us a way of obtaining the numerical values of these coordination numbers. 
We have, for simple cubic systems, z=6, y=24, y1=0 and y2=264, and for body-centered cubic 
systems, 2=8, y=96, y1=144 and y2=1776. 
The calculation of the critical point for the setting in of a long distance order and the accompanying 
eg AC, of the specific heat can be carried out in the usual manner. For the special case 
=4N, the results are 


(14) 


4 1 
a at the critical point T, or 1+- 242 
32? 


4y 16 211 
~)=+ — + +}, (15) 


.34 3 6y 52 Sy2 4071 26y 28 
AC, in k per atom= +(- (16) 
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In the expansion we have considered y, 7; and y2 as of the order z, though strictly speaking we should 
not do so, as y, y1 and y2 are much larger than z. 

The fifth and the sixth terms in the above expansions are directly affected by what values A; and Ax. 
take. It is easily seen that owing to the large value of y, y: and y2 the convergence of the above series 
is slow. For example, the third term in the bracket in (16) is 2/(3z*), amounting to 0.0104 for body- 
centered cubic systems, while the fifth and the sixth terms for the same lattice structure are 0.0173 
and 0.0233 respectively. Thus, for a thorough investigation, a knowledge of the higher )’s is necessary. 
Neglecting the dotted terms in (15) and (16), we have 


—a at T,=0.838, AC, in k per atom =1.92 for simple cubic systems, 
—a at T.=0.598, AC, in k per atom =1.87 for body-centered cubic systems, 


which are not significantly different from values obtained earlier.! 

In any effort to coordinate the change of configurations with the change of the vibrations of the 
atoms in the binary solid solution,® it is necessary to construct the partition function with a factor Y 
corresponding to the vibrations of the solution, Y depending in general on 6, X and 7. Thus in the 
construction of the partition function it is necessary to know explicitly the function f(V, X). 
This function can be found from the }’s as follows. 

Replacing the summation (1) by the maximum term, and equating its logarithm to (2), we have 


log f(N, X)+Xa=1 aN (17) 
where X is a function of 6, 6’ and a@ given by 
(8/AX) log f(N, X)+a=0. (18) 


An alternative way of looking at the Eqs. (17) and (18) is to consider (17) as an equation between f, 
6, &’ and X, with a considered as a function of 6, 6’ and X determined by (18). This provides us at 
once a method of determining f. After eliminating a from (17) and (18), we obtain 


2N 2N /\\1—n( — 
where 


1 1 
——{(1 {(1 ~00"(1 —0)(1—0"), 
12 42 


1—6’)4—6’4}, 2 


1 
30 6\s 3s 2? 


1 


The higher C’s can be obtained when }; and the higher \’s are determined. A similar calculation of f 
from an expression of F(@:7) deduced from the equilibrium properties of the solid solution calculated 
by Bethe’s method has been done in reference 3. If we expand the formula for log f given there in 


® Chang, Proc. Camb. Phil. Soc. 35, 265 (1939). 
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powers of (X —}Nz66’)" and put the expansion in the form (19), the coefficients Co, --- up to C¢ are 
exactly the same as those given by (20), except that all terms containing y or 71 or yz are no longer 
present. In fact, if we put the expression for log F(#@:7) deduced from the equilibrium properties 
obtained by Bethe’s method in the form (2), the \’s are exactly the same as those given by (12), with 
y, Y1 and 2 put equal to zero. This accounts for the agreement of the present formula for log f and 
the formula given there when we remove all terms containing y or 1 or 72. It is almost needless to add 
that if we expand the values of a, and AC, obtained by Bethe’s method in powers of 2~', we obtain 
(15) and (16) after removing all terms in them which contain y or y; or ¥2. 

In conclusion, the writer wishes to express his thanks to Dr. J. S. Wang of Southwestern Associated 


University and Mr. S. K. Kao of this university for their interest and helpful discussions. 
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Second- Neighbor Interactions and the Critical Behavior of Binary Solid Solutions 


T. S. CHANG 
Department of Physics, Central University, Chungking, China 


(Received December 5, 1940) 


The Thiele semi-invariants 71, y2 and 3 of a binary solid solution are calculated with the 
assumption that there are second-neighbor interactions as well as first-neighbor interactions. 
The results give a decrease of the critical temperature 7. at which the superlattice sets in, 
an increase of the local order o at T., an increase of the discontinuity of the specific heat, thus 
in general agreement with the results obtained by applying Bethe’s method to the same 


problem. 


HE configurational partition function 
F(6@, T) of a binary solid solution consisting 
of N4A atoms and N—N,B atoms and inhabit- 
ing a lattice of the type AB (a lattice consisting 
of two interpenetrating sublattices, referred to 
as a and B sublattices) for configurations with 
3N0A atoms on one of the sublattices, say a 
sublattice, can be put into the form! 
-—) 


where 6’ stands for 2N4/N-—6 and y, are to be 
determined. When all interactions between the 
atoms are those between the first neighbors, the 
first six y, have been determined and are exactly 
(Vaat Ves—2 Vas)” times the A, given in 
reference 1, Vas, Vas and Vag denoting the 
interaction between AA, BB and AB pairs of 
first neighbors. As shown there, the convergence 
of the series }-yn(”!)-'(—RkT)— at the critical 
-_1See the preceding paper entitled “A note to Bethe- 


Kirkwood’s partition function of a binary solid solution,” 
J. Chem. Phys. 9, 169 (1941). 


point is slow, and for an accurate determination 
of the behavior of the setting-in of a long distance- 
order at the critical point, it is necessary to find 
7, Ys, etc. Before we undertake to find y;, ys, 
it is advisable to know what effects the intro- 
duction of second- or third-neighbor interactions 
has upon the lower y’s, say 71, y2 and 73. If the 
effect is large, it will be unwise to calculate 77, ys 
with second- or third-neighbor interactions 
neglected, since they will depend on the second- 
or third-neighbor interactions to a greater extent. 

In the following, we introduce second-neighbor 


(2) 
interactions oe Veep and = and calculate 71, 


2 and 7; following a method given in reference 1. 
The same problem has been studied by using 
Bethe’s second approximation.** The result 


given there is that, if vat Vee~3Ves 4 is also 


positive, just as Vaat+Veeg—2V4z, the critical 
temperature 7, is lowered, the local order o at 
T., say o,, is increased, and the discontinuity 


2H. A. Bethe, Proc. Roy. Soc. 150, 552 (1935). 
3T.S. Chang, Proc. Roy. Soc. 161, 546 (1937). 
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CRITICAL BEHAVIOR OF BINARY SOLUTIONS 


AC, in the specific heat considerably increased. 
The mere knowledge of 7:1, y2 and 3 for cases 
with or without second-neighbor interactions 
cannot give us the precise amount of changesin 7, 
or o, or AC,, since the series >y,(n!)—"(—kT)- 
converges slowly in both cases, meaning that the 
higher y’s are equally important in determining 
T., ¢- and AC,. But as can be seen below, the 
results calculated with y1, y2 and ys and with 
the higher y’s neglected points clearly to a 
decrease of 7., an increase of o, and an increase 
of AC,, which is in agreement with the results 
obtained by Bethe’s method. 

Denoting the number of AA pairs of first 
neighbors by X, the number of AA pairs of 
second neighbors by Y, —(Vaa+ Ves—2Vaz)/kT 


by a, and by ta, we 


may write F(@, T) as 


1N\/ 4N 
)( 
2! 


where M,=((X+ty)), Me2=((X+tY)), 
(u) denoting the a priori average of u over all 
configurations with }N@A atoms on the a 
sublattice. Introducing {4 as the number of A 
atoms in the site a of the a sublattice and m as 
the number of A atoms in the site } of the 6 
sublattice, we have 


X= Lab aN, (3) 
Y=i{> 08 Maa’ } (4) 


where \@ is a first-neighbor matrix which is 
defined to be unity if the sites a and 3 are first 
neighbors and zero otherwise and paa’ (Or pes’) is 
a second-neighbor matrix: which is defined to be 
unity if the sites a and a’ (or b and 0’) are second 
neighbors and zero otherwise. Further we 
introduce the coordination numbers 21, 22, 2112 
and 2222 defined by 


21= Dia Ld Aad, 


22> Kaa’ > Mob’, 


(S) 
2112={ — a,a’,b AabAa’bMaa’, 
112 N 


2 
(x) 


EE EEE 
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where a, a’ and a” as well as } and b’ denote 
different sites. Then it is possible to calculate 
M,, Mzand M; exactly in the same way as given 
in reference 1, i.e., we first replace (f--+7---) 
by (¢-++)(n-++), and then (fq) by 61, Safa’) by 
62, by 83, by 41’, (none) by 42’, etc., 
obtaining finally only quantities such as the 6’s 
and the summations of products of \.s’s and y’s. 
The 6’s are given in reference 1, and regarding 
the summations of products of Aq’s and p’s, we 
meet the following kinds in the calculation of 
M,, and M3, i.e., 


(i) Summations of products of one \q and one yp, 

(ii) Summations of products of two Aw’s and 
one 

(iii) Summations of products of one Aw and 
two 

(iv) Summations of products of two or three Aq’s, 

(v) Summations of products of two or three y’s. 


By following a method sketched in the reference 
1, we can express summations of products of 
the first and the third types in terms of 2; and 
22, those of the second type in terms of 21, 22 and 
Z112, those of the fourth type in terms of 2;, and 
those of the fifth type in terms of z2 and 222. 
After obtaining these summations and also the 
6’s, we obtain M,, M2 and M; in terms of 6, 6’ 
and the coordination numbers 21, 22, 2112 and 
2222. 

If we proceed further to find the higher M’s, 
we must introduce other coordination numbers. 
It can be seen that quite generally, each co- 
ordination number is connected with a summa- 
tion of the products of Aw’s and y’s, say the 
summation of 


in which, if we consider \;; as something con- 
necting the subscripts i and 7 and y,; as some- 
thing connecting the subscripts k and /, we can 
find two or more chains of q's and y’s, with no 
Xa» Or w in common, connecting any one of the 
subscripts to any other. Each of the subscripts 
occurs therefore at least twice, and the whole 
summation is of the order N/2. We can thus 
define the corresponding coordination numbers 
as 2/N times the respective summations. For 
example, in finding M4, we meet the coordination 
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numbers 


21122 


21111 =— DL 
N 


kT | 
(kT)? 
¥3 
(kT)* 


where terms of an order smaller than N are not 
written out. It can be noted that all the y’s are 
linear in the coordination numbers just as the 
in reference 1. 

The numerical values of such coordination 
numbers for simple cubic systems and body- 
centered cubic systems are as follows: 
for simple cubic systems 

2,:=6, Z.=12, 2112 = 24, Zo00=0; 
for body-centered systems 
2:=8, 22=6, 212=24, Ze2=0. 


2 
> 


It is needless to say that the number of co- 
ordination numbers increases very rapidly when 
we go into the higher M’s. 

Substituting the found values of Mi, M2 and 
M; into (2), expanding in powers of N-, and 
then comparing (2) with (1), we get 


= +4 — 0)? +0'2(1— 0’)? }}, 


—26)(1 — 26’) + ] 
+1 —0)2(1 — 0’) 2(1 +3 (6) 


Since Z222 is zero for both cases, we shall leave 
out the term in —y3/(Rk7)* containing 2222. 
Because of the large value of 22 and 2112, the 
effect of increasing t from 0 to some value such 
as 0.1 is already noticeable. For example, for 
simple cubic systems, at 0’=%, Ys is 
changed by 1.44 times. 

The calculation of the critical point JT, and 
the discontinuity AC, in the specific heat can be 
carried out in the usual manner. For the particu- 
lar case N4=3N, we have 


1 22; —2ef? 2(21 +22t*)? 


aat 7.ora,=— 


3 


4 
1+ 


21 +20? 162,;+ 182;10t (21 +22t*)? 
6 feel, (8) 
(z;— Zot)? —Zel)4 


AC, in k per atom =| 1+ 


The local order o can be defined by X = }.Nai(1—c) 
for this special case, X being the equilibrium 
value of X for the binary solution. Now X is 
given by (0/da) log F(@,T), where in the 
differentiation ta is kept constant. Hence we 
get at T., where =}, 


1 1 1 2112 


8 2 


and accordingly, 


a 12112 


| 4 8 21 


1 22112¢ 
1+ (9) 
(21 — Zot) (31 — Zot)? 21(21 — Zot) 


Thus, if ¢ is positive, T. is decreased, o-. 
increased, and AC, increased, but the increase of 
AC, does not seem to be large enough to bring 
the calculated results into agreement with 
experiments. The above expansions are in series 
of (s1;—2et)-!, and this makes the increase or 
decrease at once obvious. We can, of course, 
expand in z;"', and then the increase or decrease 
will be shifted to a lower term. 

Finally, the writer wishes to thank Mr. S. K. 
Kao for helpful discussions. 
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Granulation, Phase Change, and Microstructure 


Kinetics of Phase Change. III* 


MELVIN AvRAMI 
Columbia University, New York, New York 


(Received August 6, 1940) 


The theory of the preceding papers is generalized and 
the notation simplified. A cluster of molecules in a stable 
phase surrounded by an unstable phase is itself unstable 
until a critical size is reached, though for statistical reasons 
a distribution of such clusters may exist. Beyond the 
critical size, the cluster tends to grow steadily. The 
designation “nuclei” or “‘grains” is used according as the 
clusters are below or above the critical size. It is shown that 
a comprehensive description of the phenomena of phase 
change may be summarized in Phase Change, Grain Num- 
ber and Microstructure Formulas or Diagrams, giving, re- 
spectively, the transformed volume, grain, and micro- 


structure densities as a function of time, temperature, and 
other variables. To facilitate the deduction of formulas for 
these densities the related densities of the ‘‘extended”’ — 
grain population are introduced. The extended population 
is that system of interpenetrating volumes that would 
obtain if the grains granulated and grew through each 
other without mutual interference. The extended densities 
are much more readily derivable from an analysis of the 
fundamental processes of granulation and growth. It is 
shown that, under very general circumstances, the den- 
sities of the actual grain population may be expressed 
simply in terms of the extended population. 


A. NUCLEI AND GRAINS 


N the previous papers' we have developed a 
treatment of the transformation mechanism 

which we may now generalize. At the same time, 
it is desirable to simplify some of the previous 
notations. To recapitulate: It appears that the 
grains of a new phase generally originate at 
certain centers or heterogeneities in the old 
phase. In these it is desirable to distinguish 
between two stages. For even though, in the 
supercooled (or superheated) condition, the 
internal free energy of the new phase is lower 
than in the old, the additional opposing free 
energy change, involved in the creation of an 
interface, makes a cluster of the new phase below 
a certain critical size (function of the thermo- 
dynamic conditions) incapable of steady growth. 
Only if a cluster larger than the critical size 
occurs, e.g., through spontaneous or stimulated 
fluctuations, will it grow steadily. 

Now, in previous investigations,! we have used 
the name “‘germ nuclei” to refer to such activable 
regions before they reach the critical size, and 
“growth nuclei” for those which have reached 
and passed the critical size. It is possible, and 
perhaps desirable for psychological reasons, to 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

1(a) J. Chem. Phys. 7, 1103 (1939); (b) ibid. 8, 212 
(1940); referred to as K.P.C. I and II, respectively. 


simplify the terminology by calling the former 
simply “nuclei” and the latter “grains’’ (or ‘‘new 
grains” if necessary for specificity). 

This we shall do in the future, confining the 
term ‘‘nucleus’’ to a cluster of the new phase smaller 
than the critical size, the term ‘‘grain’’ to a cluster 
larger than the critical size. The region (originally 
a nucleus) around which a grain has grown may 
be referred to as the ‘‘grain center.’’ A nucleus or 
grain may have foreign matter (inclusion) as a 
core, or it may form on a foreign surface and 
thus not be entirely surrounded by the old phase, 
but, in any case, the specifications for the critical 
size may be worked out.'™ It has been suggested 
that the process of passing from the nucleus to 
the grain stage shall be termed ‘‘granulation.’’t 
The “granulation rate’ is the number of nuclei, 
per unit time and unit volume, which ‘‘pop,”’ 
i.e., pass through the critical size, and, becoming 
grains, embark upon a career of steady growth. 

Several recent investigators? have discussed 
the necessary existence, in equilibrium, of clusters 
of one phase in another before the temperature 


+ This replaces the term ‘‘nucleation” which was used 
in K.P.C. I in conformity with the current somewhat vague 
terminology. It seems logical, to confine the term nucleation 
to the process of becoming a nucleus. 

2 In particular J. Frenkel, J. Chem. Phys. 7, 538 (1939); 
W. Band, ibid. 7, 324 (1939); J. E. Mayer and S. F. 
Streeter, ibid. 7, 1019 (1939). A connection has also been 
established with the general statistical mechanical treat- 
ment given by Mayer and his collaborators. 
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for macroscopic transition is reached.{{ These 
“‘heterophase clusters” are nuclei in our sense of 
the word. In the case of a substance free from 
foreign particles, they are the only nuclei present. 
Their typical size-frequency density (i.e., number 
of given size per unit volume), at a certain tem- 
perature 7, above the transition point, is indi- 
cated by the solid curve in Fig. 1. Below the 
transition point, the critical size diminishes from 
infinity to some finite value. Thus, if we quench 
suddenly from 7; to a supercool temperature 7», 
for which the critical nucleus size is v,, all clusters 
larger than v,, present at 7), are effective grains 
and start growing steadily. Those below the 
critical size shift their frequency statistically and, 
it can be shown, generally increase in number, so 
that at a later time the frequency is given by a 
curve like the dotted one in Fig. 1. The general 
kinetics of the processes involved will be analyzed 
in a later paper. For the present, the underlying 
physical situation may be described as follows: 
Any particular nucleus below the critical size has 
a higher tendency to dissolve than to grow. But 
when we consider the entire group of a particular 
size, we find that the tendency of its majority 
to dissolve is more than balanced by the minority 
tendency of the larger group of next smaller size 
to grow. The granulation rate, at any time, is the 
number of clusters per unit volume which pass 
through the “bottleneck” v, per unit time. 
Sometimes, the distribution with which we 


start was not in equilibrium at the pretreatment | 


tt For simplicity of speech, and in conformity with the 
usual experimental conditions, we consider pressure kept 
constant and temperature as the controlling variable. Also, 
we take the new phase as the one stable at lower tem- 
peratures, and therefore speak of supercooling. The 
generalization can be made without difficulty. 
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temperature. Also, it may be non-uniform, so 
that the resulting density varies in space. The 
considerations in this paper will, however, apply 
generally to these and other industrially im- 
portant cases where the transforming substance 
contains foreign heterogeneities which catalyze 
the reaction by providing adsorbing-activating 
surfaces on which nuclei and grains form. These, 
like the heterophase clusters, may be present 
from the beginning of the reaction, and are 
therefore determined by its “‘prehistory’’; or they 
may be introduced physically, or formed chem- 
ically, during the reaction, and are therefore func- 
tionals of its history. 


B. VoLUME, GRAIN, AND FREQUENCY 
DENSITIES 


We may arrange observations on the kinetics 
of transformation under the following headings. 
First we consider the volume density V of the new 
phase (i.e., fraction of unit volume of space 
occupied), as a function of time ¢, position P, 
initial temperature To, and as a functional of the 
prehistory and of the subsequent type of thermal 
treatment 7(t, To, P). For given prehistory and 
type of history (e.g., isothermal, linear cooling, 
“seeding” at a certain time, etc.), and for the 
unit volume around a given point, a three- 
dimensional diagram relating V, t, and JT) may be 
constructed. We shall call this the ‘‘ Volume Dia- 
gram” or “Phase Change Diagram” in contrast 
with the familiar Phase Equilibrium Diagrams 
to which it reduces in the cross section at infinite 
time (system with isothermal history). The ana- 
lytical representation of 


V=V(t, To, P) 


will be called the “Phase Change Formula.”’ It is 
to be remarked that V corresponds to the vari- 
able, ‘“‘degree of advancement”’ of a reaction, 
used by De Donder-Rysselbergh* and Guggen- 
heim‘ in various thermodynamic studies. _ 

A second important type of relation is that 
giving the grain density N (i.e., number of grain 
centers per unit volume)* as a function of t, To, 

3 De Donder-Rysselbergh, Affinity (Stanford University 
Press, 1936). 

‘ Guggenheim J. Chem. Phys. 7, 103 (1939). ; 

*This is a reversal of the notation in the two earlier 


papers, where N’ was used to designate grain density and 
N was used for nucleus density. Since we shall deal 
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KINETICS OF PHASE CHANGE 


P, T(t, To, P), and the prehistory. Again, for 
given prehistory, type of history, and location, 
we may construct a three-dimensional ‘‘Grain- 
Number Diagram” between N, t, and To. The 
analytical representation of 


N=N(t, To, P) 


will be called the ‘‘Grain-Number Formula.” 

The grain-number diagram or formula is of 
great importance in the theory of microstructures 
resulting from given treatment. For a given type 
of crystal habit, a quantitative description of 
microstructure may be given by the volume- 
frequency density F(v)dv, i.e., the number, per 
unit volume, of centers of grains (or nuclei) with 
volume between v and v+dv. This, too, is a 
function of time, position, initial temperature, 
etc. Fixing the prehistory, history, and location, 
we still have one more variable here, and it would 
be necessary to use a four-dimensional plot to 
represent the frequency density as a function of 
both time and temperature. If, however, we fix 
the initial temperature, we may set up a three- 
dimensional diagram between F, v, and ¢t. This 
may be called the microstructure-time diagram. 
A more important diagram will be that giving 
the volume-frequency density, after infinite (or 
very long) time, as a function of initial tem- 
perature. This we shall designate the ‘‘ Micro- 
structure Diagram.” The general formula for F 
will be called the ‘Microstructure Formula.” By 
definition, the relations 


v= f Fede (1) 
and 
V= f vF(v)dv (2) 


must subsist. 

In some cases (crystalline aggregates), the 
orientation 7 (abbreviation for the two direction- 
specifying coordinates) of the individual grains 
with respect to some fixed direction in space may 
be distinguished. In others, where the grains are 


primarily with grain density, it is desirable to make this 
simplification even at the risk of some difficulty in com- 
paring with I and II. 

+The ‘‘nucleus density” has the obvious correlative 
definition fer (v)dv. Formula (2) is correct only if the unit 
volume is large compared with the grain sizes or the space 


gradient of grain density; see Eq. (5) below. 
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isotropic and possess complete spherical sym- 
metry (drops and bubbles), this distinction is 
lost. In the anisotropic grain case, a more detailed 
description of the system may be given by taking 
into account the distribution in orientation of 
the grains. This leads to an orientation-dependent 
grain density 


N(n, t, To, P)dn, 
and an orientation-dependent frequency density 
F(n, t, To, P)dn. 


The previously given NV and F may be regarded 
as orientation-averaged densities, and it is with 
these that we shall work in the following unless 
otherwise specified. We note that all the sub- 
sequent analysis may be carried through equally 
well with the orientation-dependent densities. 
The deduction of formulas for V, N, and F is 
considerably facilitated by the conception of 
related ‘‘extended densities’”’ which we shall now 
develop. As a transformation proceeds, two phe- 
nomena occur: (1) Some nuclei are reached and 
enveloped by growing grains and become inoper- 
ative for production of new grains. (2) The 
growing grains impinge upon one another and 
interact; the possible forms of this interaction 
will be discussed later. Those grains which would 
have appeared up to any given time, if their 
nuclei did not lie in occupied regions, will be 
referred to as “phantom grains.” ! The grains 
(phantoms included), taken in the fullest extent 
they would have, had their growth been unim- 
peded and unmodified by impingement upon one 
another, will be referred to as ‘‘extended grains.” 
To visualize the extended-grain population, we 
may imagine a representative geometrical system 
of interpenetrating volumes obeying exactly the 
same laws as the physical system except that the 
volumes granulate and grow through each other. 
For brevity, we shall use the adjective ‘‘ex- 
tended,” and the superior symbol ~, to designate 
any attribute of the representative geometrical 
system.{t Thus, the “extended volume-frequency 
density” F(é)dé is the number per unit volume 
of (phantom-included) grains or nuclei of ex- 
tended volume between 4 and i+dé. The 


“extended grain density” N is the total number 


tt This is a simplification of the notation in the previous 
papers in which the subscript x was employed. 
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larger than @, in the unit volume, i.e., 
N= f F (6) dé. (1’) 


Similarly, the extended volume density V is the 
fraction of unit volume of space occupied by the 
extended population. By definition then, over a 
volume large compared with the size of an 
average grain 


Evidently V can become larger than unity. 

(1’) and (2’) enable us to compute N and 7 
when F is given. If, on the other hand, N and the 
law of growth of the individual extended grain 
is easily found (a frequent case), V and F may be 
calculated therefrom. V, N, and F may then be 
calculated, so that a knowledge of N and the law 
of growth of the individual grain is all that is 
required for the determination of the kinetics of 
a granular aggregate (sections C and D). 

The introduction of the extended densities 
considerably facilitates the theoretical treatment 
of transformation kinetics. The problem can now 
be split into two parts. First, an analysis of the 
microkinetic processes of granulation and growth 
leads to the determination of the extended 
densities as a function of position, time, initial 
temperature, and as a functional of the prehistory 
and type of history.* Then, by the formulas to 
be developed in this paper, the observable 
macrokinetic quantities, V, N, and F, of the 
actual grains of the new phase may be deduced. 


C. AVERAGE DENSITIES AND THEIR RELATIONS 


It is necessary at this point to refine our 
concepts somewhat. So far, we have applied the 
term ‘‘density”’ to numbers or quantities found 


* We shall show elsewhere that information about pre- 
history and type of history corresponds to boundary con- 
ditions for the system of differential equations governing 
the nucleation and growth processes of the extended 
population. This system is also representable by a partial 
differential equation of the same form as that which appears 
in the theory of heat conduction in a heterogeneous bar or 
electric propagation in a non-uniform cable. When this is 
set up and solved under appropriate boundary conditions, 
it yields the extended densities. On the other hand, the 
latter may often be found directly by elementary methods, 
as in K.P.C. I and II, where assumptions holding in special 
cases are made. 
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in unit volume, without pausing to consider the 
dependence upon the absolute size of the unit 
volume. We obtain a useful, statistically suitable, 
definition of density for point distributions in the 
following manner: Subdivide the space contain- 
ing the points or grain centers into ‘‘physically 
small’’ cells 2 with positions indicated by sub- 
script capital letters. Each is further subdivided 
into ‘“‘geometrically small’”’ subcells Aw (positions 
designated by small letters) which may become 
infinitesimal. Let Ng denote the number of 
centers contained in Qg. The average 


N(Q)=Ne/2q 


will be referred to as the ‘“‘average density around 
Q” or, simply, the ‘‘density around Q.” The 
absolute size of Q is arbitrary within the limits 
set by the following conditions: (1) It shall be 
large enough to contain sufficiently many centers 
so that changing its dimensions continuously 
alters the density within only slowly. (2) It shall 
be small enough so that the average number of 
centers in the subcells Aw does not vary appre- 
ciably from one side of @ to the other. A further 
analysis of the concept of average density leads 
to some very interesting questions which we shall 
not discuss here. For the present purpose, the 
rough verbal restrictions given in (1) and (2) 
will suffice. We merely note that, for a sufficiently 
fine-grained distribution, the Q will be small 
enough to serve, approximately, as elements of 
integration. We shall assume this in the following, 
replacing the individual 2 by dQ whenever a 
summation is required. 

It is clear that, though the average density 
may vary in a complex manner from one region 
of space to another, the centers within a cell may 
be very simply distributed. A particularly im- 
portant case, which includes most physical situ- 
ations, is where the numbers in the subcells are 
distributed about the local average N according 
to the laws of chance; i.e., the probability of 
finding m centers in any Aw is given by the 
Bernoulli distribution function: 


( 


The distribution is then said to be “locally 
random.” This includes the many important 
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systems where the gradient of grain density is 
very sharp, even infinitely sharp, as when the 
nuclei are distributed upon certain surfaces 
(walls or grain boundaries of the old phase). We 
shall investigate the case of local randomness in 
detail below. 

We return now to consideration of the quan- 
tities discussed in the previous section. It is clear 
that practically useful descriptions of frequency 
and grain distributions will be in terms of average 
densities. Also, a theoretical treatment is likely 
to yield information only about such average 
densities. In related fashion the probable or 
average “extended volume density around a point 
P” V(t, P) may be defined by the average 
number of extended grains which have grown 
through dQp.t This may be expressed in terms 
of the integral of N over the entire region of 
influence of P 


Vit, P) = f Nt’, Q)dQo, (3) 


where ?’ is a retarded time something like that 
appearing in field theories (electromagnetism, 
etc.). The retardation t—?’ is the time required 
for the “propagation” or growth of a grain with 
linear velocity G from its origin at Q to its ter- 
minus P at time ¢. G is a function of temperature, 
concentration, etc., i.e., in general of position and 
time. Expressing Q in terms of spherical coor- 
dinates r, 6, g about P as origin, we obtain the 
more explicit form (letting dJ=sin 6d6d¢) 


Vit, P)= f f f (t’, r, 0, (3’) 
0% Yo 


where * denotes the radius of the region of 
influence for P (i.e., t’=0). If external conditions 
remain constant in time so that G is a function 
only of position p, 6, ¢ then 


T It is evident that, upon integration over any macro- 
scopic volume, the density so defined yields the correct 
total extended volume. Thus, for instance, the average 
extended volume in a macroscopic unit volume is the 
integral over the volume 


This is the integral equivalent of formula (7) in K.P.C. I. 
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If the temperature is changing so that G is also 
a function of the time w, ¢’ must be found from 
the conditions that, in the integral of the dif- 
ferential equation 


dp= —G(p, 8, 9, w)dw, 


w=t at p=0 and w=?’ at p=r. This then ex- 
presses ¢’ as a function of r, t’(r, 0, g; t). On the 
other hand, if we so wish, the differential equa- 
tion may be integrated to give r as a function of 
t’, r(0, ¢, t’; t), and r may be replaced in (3’) by 
t’ as an independent variable. 

For anisotropic grains G depends also on the 
grain orientation, so that N(¢’, Q) in the integral 
of (3) should be replaced by an average over the 
densities in the various orientations, each asso- 
ciated with a different value of the retarded 
time i.e., 


or 

V(t, [roe r, 0, g)dnr*drdy. 
(4’) 


This may be applied to various cases of distribu- 
tion of the grain centers, e.g., uniformly random 
distribution in space; random distribution on . 
foreign surfaces (ingot walls in solidification, elec- 
trode plates in electrolytic deposition, etc.); 
random distribution on internal grain surfaces of 
the old phase as in grain boundary nucleation. 
We shall not discuss the general problem of 
evaluating (4’) here, but three limiting cases 
previously considered! recommend themselves 
to our attention. In the first, where the growth 
is polyhedral, i.e., approximately equal in all 
directions, t,’, is the same for all directions, and 


f N(n, ty’, Q)dn= NU’, Q) 


with G in ¢’ representing the rate of radial growth 
of the pseudospherical grain. In this case, as is 
to be expected, (4) reduces to (3). 

In the second case, that of plate-like growth, 
the grain grows primarily parallel to a single 
plane, i.e., its (small) thickness h does not change 
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appreciably while it grows planarly from its 
center with radial velocity G. If the plate grows 
in the shape of a circular sector with constant 
angle y (radians), the cross section of the ad- 
vancing face of the plate at distance r from its 
origin in yrh. Supposing, also, that the distribu- 
tion of orientations is random, the probability of 
a grain with center at Q intercepting the element 
of volume dQ; at P, is the ratio of the area inter- 
cepted by the plate on the sphere of radius 7 to 
the area of the sphere itself, i.e., 


hyr/4ar?=hy/4ar, 


where 7 is the distance PQ. Thus, out of the total 
Nt’, Q) of grains which had appeared around 
Q at t’, only the fraction 


hy 


is effective. Similarly, in the third case, that of 
lineal or needle-like growth at the rate G with 
approximately constant cross section a, the 
probability of a grain, originating at Q, inter- 
cepting dQp is 

a 


4rr? 
and the effective orientation-averaged grain 


density is 


a 
—WNi', Q). 


For brevity, let us designate the factor in front 
of N in (4’) by or*/4x so that 


F(t, P) J J "PIN, 1, 0, 


where 


o=4r k=3 for polyhedral growth 
o=hy k=2 for platelike growth 
o=a k=1 for lineal growth. 


Integrating (4’’) by parts with respect to 7, and 
remembering that 


/k=i(t, t’) 
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for the three types of growth, we have 


V(t, P)= ff t’) 


ON aN at 
— }drdy 
or ar 


or replacing 7 by ?’ as an independent variable 


aN aN 
x(—o+— Java. 
or ot’ 


In the special case where N is independent of 
position in the region of influence of P (i.e., 
approximately uniform grain density), this re- 
duces to 


t aN 
V(t, P)= f a(t, yw. (5) 


which is equivalent to (6) of K.P.C. I with the 
mechanism of granulation assumed there. 

(3) or (4) expresses V in terms of N (and ?¢’, 
whose relation to ¢ depends on the law of growth 
a(t, t’) of the individual grain). It remains to 
express F also in terms of N and é. This may be 
done as follows: A grain of volume @ at time / 
comes from a nucleus which granulated at a 
previous time ¢’. This time can be determined 
from # and ¢ when the law of growth A(t, ¢’) is 
known. Now F(i, t)dv equals the number of 
grains which granulated at the earlier time ¢’ in 
an interval dt’ corresponding to di, i.e., 


— t)—dt’ 
ot’ ot’ 
From this we obtain the extended microstructure 


formula 
0 
—F (i, t)di= - (6) 
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We shall now see that, for very general cir- 
cumstances, it is possible to express V simply in 
terms of V. But first we must consider the pos- 
sible types of grain interaction. When two grow- 
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ing grains of a new phase impinge, several things 
may happen. They may coalesce to form one 
larger grain, as frequently happens with gaseous 
bubbles and liquid drops. They may adhere to 
form a compound grain, ceasing growth at the 
common interface while continuing to grow 
normally elsewhere; this is usually true of 
crystalline grains. Again, a third mechanism may 
operate, especially if the surrounding phase is 
fluid, i.e., the interacting grains may rebuff each 
other, pushing away and continuing to grow as 
though there were no impingement. In this last 
case V will be a good approximation to V up to 
the last stages of transformation. In some situ- 
ations, all three processes might operate to 
varying extents, together with motion of the 
grains (and nuclei) due to gravimetric or con- 
vective forces. This would complicate the trans- 
formation kinetics with considerations of me- 
chanics. Here, however, we shall confine our- 
selves to stationary particles with the adhering- 
grain interaction alone. This assumption is 
usually valid for the growth of crystalline phases. 

In the previous paper it was shown that, when 
the distribution of nuclei is uniformly random, 
V is given by 


V(t, P)=1—exp (— Vit, P)). (7) 


We shall here show that (7) is equally true when 
the distribution of nuclei (and initially present 
grains) is only locally random, and V(t, P) is 
given by (3). With this generalization most of the 
important nucleated reactions encountered in 
practice are included. 

Assuming the initial nuclei .to be locally 
random, the extended grain population at any 
later time is also distributed with local random- 
ness. Now, the ratio of non-overlapped to total 
volume of any random region, around any point 
P, is equal to the probability of untransformed 
volume there, i.e., 1— V(t, P). As in K.P.C. II, 
let us choose for the random region the increment 
of transformed volume in an element of time. 
The non-overlapped part is dV, the total dV, 
and we have 


dV/dV=1-V. 
Integrating, we obtain Eq. (7). 


An alternative derivation may be given as 
follows: Letting the average extended grain 
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density around Q, at time ?’, be designated as 
above by N(t’, Q), we may calculate the prob- 
ability that the element Aw,, at g, does not 
contain any grain centers. By the Bernoulli dis- 
tribution this probability (m=0) is 


@, 
If we let Aw, become infinitesimally small, we get 


lim =exp (—N(t’, Q)dw,). 
Aw 


This may be used to calculate the density of 
untransformed volume 1—V around a point P, 
which is simply the probability that the volume 
element around P is completely untransformed. 
This probability is the product of probabilities 
of zero grain center content in each of the 
elements dw, (at its associated retarded time 2’) 
in the region of influence of P, so that we have 


1-V(t,P)= TI 


all Q, q for ’>0 
=exp (— SN(t’, Q)d29) =exp (— Vit, P)) 


and we have justified Eq. (7). 

We have now expressed V in terms of V. The 
latter, as we have seen, can be calculated from a 
knowledge of N and @. Likewise, assuming local 
randomness, N, v, and F may all be expressed 
ultimately in terms of N and @. Thus, taking 
unit volume around a given point, the ratio of 
actual to extended granulations in an element of 
time will evidently equal the fraction of untrans- 
formed volume at that time, i.e., 


exp (—N(?’, 


or, integrating, 


.t ON 
Nit)= exp (— 


t aN 


Again, as in deriving (7), the ratio of the non- 
overlapped increment dv of an average grain (of 
size v) to its extended increment di, in an element 


|_| | 
— 
aN aN _ oN 
—=(1— V)—=exp (— V)— 
ot ot ot 
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of time dt, equals the untransformed fraction at 
that time, i.e., 


dv(t, t’) 
di(t, t’) 


=1-—V(t)=exp (— 


Integrating, we get 
00(t, t’) 
ot 


dt, (9) 


v(t, = f exp (— 


where the values of the integrand are taken along 
the (space-time) route of growth of the grain 
around its center. 

The expression for F may be obtained by 
reasoning similar to that used for (6). An average 
grain of volume v at time ¢ comes from a granu- 
lating center at a previous time ?¢’ which can be 
determined as a function of v and ¢t when the law 
of growth (9) of the actual (impinged) grains is 
known. F(v,¢t) equals the number of actual 
grains which made their appearance at time ?’ 


in the interval dt’ corresponding to dz, i.e.,* 
dv aN aN 
— F(v, t)—dt' = —dt' = exp (— V(t’))—d?’. 
at’ at’ 


Now, from (9) 


ov t t’) 
f exp (— 
ot’ ot’ dt 


so that the microstructure formula becomes 


aN 
Fle, )=— exp (~ 


f exp (— V(t)) (10) 
ot’ dt 


* Actually this is an approximation which is truer the 
more nearly ¢’ corresponds to a unique volume in the actual 
grain population, i.e., the narrower the size frequency dis- 
tribution of grains of the same ‘‘age,”’ a function we have 
not analyzed here. This approximation makes (10) less and 
less valid toward the end of transformation. For a more 
detailed analysis of this and related points see a forth- 
coming paper in Philosophy of Science (Technical Supple- 
ment), January (1941). 
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Remarks on the Fluorescence, Phosphorescence and Photochemistry of Dyestuffs 


James Franck, University of Chicago, Chicago, Illinois 
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Rosert LivinGston, University of Minnesota, Minneapolis, Minnesota 
(Received July 26, 1940) 


Many of the apparently conflicting facts of the photo- 
chemistry as well as of the phosphorescence and of the 
fluorescence quenching of dyes can be given a rational and 
unified interpretation if it be assumed that an elec- 
tronically excited dye molecule can go over, by a process of 
internal conversion, to the electronic ground state of a 
reactive energy-rich tautomer. The endothermic reversion 
of such a tautomer to the electronically excited state of the 
original molecule would explain the weak, temperature 
dependent phosphorescence observed in some systems. 
In the presence of a suitable reducing agent the reactive 


tautomer may be reduced to a semiquinone. This process 


and the likely subsequent reactions are sufficient to ex- 
plain the majority of dye-sensitized photo-oxidations as 
well as the photo-bleaching of dyes by reducing agents. 


T is well known that many dyestuffs which are 
able to fluoresce in solution are efficient 
sensitizers for photochemical reactions, while 


In the case of chlorophyll, and possibly a few other dyes, 
it appears to be necessary to assume that the tautomer and 
a normal dye molecule can undergo a process of dispro- 
portionation. The assumption of this step leads directly 
and simply to an explanation of the rapid reversible 
bleaching of chlorophyll which is consistent with the other 
known photochemical and optical properties of this sub- 
stance. In addition to the more familiar reaction steps 
which have been assumed as follow reactions in setting up 
the various mechanisms, it is necessary to include a step 
in which the HO; radical reduces a partially oxidized dye 
molecule (i.e., radical) back to the normal state. The 
several reaction steps assumed in these various schemes 
appear to be energetically feasible. 


most of those which show no sign of fluorescence 
in solution are extremely inefficient sensitizers." 
1 Hurd and Livingston, J. Phys. Chem. 44, 865 (1940). 
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The evident and commonly accepted expla- 
nation of these facts may be stated as follows. 
The ability of a substance to fluoresce is proof 
that the excitation energy of its electronic system 
remains intact long enough for this energy to be 
emitted, wholly or in part, as light. However, the 
oscillation energy which is associated with the 
electronic energy is quickly dissipated by impacts 
with the solvent molecules. If 7 is the natural 
lifetime of the excited state, as calculated (to a 
good approximation) from the integrated ab- 
sorption coefficient, and a is the fluorescent yield, 
then the product ar is the actual lifetime, which 
would be measured by a direct observation. A 
collision between an excited dye molecule and a 
molecule of a reactive substance may result in the 
utilization of the energy of excitation as energy of 
activation of a chemical reaction, with resultant 
quenching of the fluorescence. The probability of 
such collisions, and therefore the efficiency of 
utilization of the energy of excitation for chemical 
purposes, increases with an increase in the actual 
life time (a7) of the excited state. In the case of 
those dyes which show no trace of fluorescence, 
even under conditions where the transfer of 
energy of excitation into chemical energy is 
extremely improbable, one must assume that the 
light energy is dissipated into the degrees of 
freedom as thermal energy directly after the 
absorption act. The light energy is thus lost for 
chemical purposes. 

This interpretation, while generally correct, 
offers no explanation of the great difference in 
ability to fluoresce which the several dyes exhibit. 
Furthermore, it should be amplified in some 
respects and restricted in others. For example, 
the energy absorbed by some nonfluorescent dyes 
is not directly and completely dissipated as heat. 
It should also be noted that a collision between 
an excited dye molecule and one of the substrate 
is not a sufficient condition for the immediate 
occurrence of a chemical reaction. This point has 
been discussed more fully by Franck and 
Herzfeld,? who suggested the analogy between an 
excited dye in state of impact with a reactive 
molecule and an activated complicated molecule 
which is about to undergo a monomolecular 
dissociation. Just as there is in the second case a 


* Franck and Herzfeld, J. Phys. Chem. 41, 97 (1937). 


finite interval between the time when the com- 
plex molecule acquires the energy of activation 
and the time when it uses this energy to dissoci- 
ate, so there is a lapse of time after the dye 
molecule absorbs the light before it can transfer 
the absorbed energy to a second molecule which 
is able to undergo a chemical reaction. It is an 
interesting consequence of this view that some 
cases are to be expected where the fluorescence 
will not be completely quenched either by very 
frequent collisions with the chemically sensitive 
molecule or even by the formation of a complex 
with this molecule. This residual or ‘‘unquench- 
able” fluorescence will be strongest in those cases 
where the time required for the attainment of 
that configuration which leads to the occurrence 
of the chemical reaction is greatest. 

To arrive at any conclusion as to why the 
fluorescent yields of the various dissolved dves 
differ so greatly it is necessary to consider the 
several ways in which they can lose their exci- 
tation energy. If the dye is dissolved in a pure 
chemically stable solvent, whose molecules are 
unable either to react chemically with the excited 
dye molecules or to take over the excitation 
energy into their own electronic systems, a 
fluorescent yield less than unity indicates that 
excitation energy is either converted into heat 
or is stored in the dye molecule, very probably by 
a transition into a tautomer of higher energy 
content. It is known, from numerous experi- 
mental studies of the quenching of the fluores- 
cence of mono- and diatomic molecules, that 
energy cannot be converted,*‘ by collisions of the 
second kind, from electronic energy directly into 
rotational or translational energy in amounts 
much greater than the thermal (equilibrium) 
energy of these degrees of freedom. Furthermore, 
the direct conversion of electronic energy into 
oscillatory energy is possible only under the 
special condition‘ that the potential energy curves 
of the acceptor molecule are greatly changed by 
the forces (electric or chemical) exerted by the 
excited molecule. If such forces existed between a 
dye molecule and the molecules of the solvent 
(which are constantly in contact with it) the 


3 Franck, Naturwiss. 14, 211 (1926). 
4 Nordheim, Zeits. f. Physik 36, 496 (1926). 
a a and Eucken, Zeits. f. physik. Chemie B20, 460 
). 
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Fic. 1. Schematic 
representation of 
the energy levels 
and transitions cor- 
responding toa sys- 
tem where internal 
conversion which 
results in the de- 
gradation of the 
energy of excitation 


is possible. 


positions and intensities of its absorption bands 
would be greatly changed compared to those of 
the same dye in the gaseous state or dissolved in a 
solvent which does not interact with it. If the 
absorption bands of a dye are the same when it is 
dissolved in several different solvents, it may be 
safely assumed that no such interaction occurs, 
since it is very improbable that the interactions 
and resulting changes in the absorption spectra 
would be the same for the different solvents. 
Since many dye solutions are known which are 
negligibly fluorescent but which do not show such 
changes in their absorption bands, it is necessary 
to assume that electronic energy can be trans- 
ferred within the dye molecule itself. This could 
occur by a sort of “‘internal impact of the second 
kind.” 

While such acts cannot occur in simple mole- 
cules, a process having these characteristics and 
restricted to complicated molecules has been 
described by E. Teller. This process may be 
described as an internal conversion of the exci- 
tation energy of the electronic system into 
oscillatory energy. We shall refer to it briefly as 
“internal conversion.” Internal conversion is 
possible only for those complicated molecules 
whose quantum states are so numerous that their 
spacing is comparable to their natural width. It 
can occur only when the molecule is in a special 
configuration which is common to the excited and 
ground states. In terms of the potential energy 
diagram, this transition occurs when the system 
is represented by a particular point lying on the 
line of intersection of the surfaces corresponding 
to the normal and excited states. The time re- 
quired to reach this special constellation will 
differ greatly for different dyes. It seems very 


6 Teller, J. Phys. Chem. 41, 109 (1937). 


probable that the magnitude of this interval will 
be strongly effected by relatively small changes 
in the structure of the molecule. Since the act of 
fluorescence takes place during this interval, the 
intensity of fluorescence must also be very 
sensitive to small changes in molecular structure. 

In condensed systems when the activated 
molecule is continuously in a state of multiple 
impact, internal conversion readily leads to the 
degradation of electronic energy into heat (see 
Fig. 1). However, it is neither necessary nor 
probable that internal conversion always leads to 
the complete dissipation of the energy of exci- 
tation as heat. Indeed, it is quite probable that 
the violent oscillations, which are a consequence 
of internal conversion, lead to the formation of a 
tautomer of the original dye molecule; the keto- 
enol transformation is a possible example of such 
a tautomeric change. 

Both the physical and chemical properties of 
the molecule may be considerably altered by such 
a tautomeric shift in the position of a hydrogen 
atom. For example, if, as a consequence of this 
shift, the system of conjugated double bonds is 
interrupted, the absorption bands may be entirely 
changed. Of course, if the shift does not affect the 
resonating system (i.e., if it occurs in another 
part of the molecule) practically no change will 
occur in ‘the absorption spectrum. The behavior 
of chlorophyll dissolved in an organic liquid will 
be discussed later as an example of the latter 
type. The energy relations of a molecule of this 
general type are illustrated by Fig. 2. In this 
diagram, the line N represents the energy of the 
ground state of the dye; E£, its first excited state; 
and N’, the ground state of the tautomer. If N’, 
as in the figure, happens to be well below E but 
still much higher than N, then the molecule in 
the tautometric state will have optical properties 
which will be similar to those which it would have 
in a metastable electronic state, with the im- 
portant difference that the tautomer may have a 
much longer life. This difference in half-life is 
likely to be especially great when the dye 
molecule is surrounded by solvent molecules. It 
is evident that a tautometric level, like N’, could 
be responsible for the phenomenon of phosphor- 
escence. The intensity of the phosphorescence 
produced in this way should decay exponentially. 
Behavior of this type has been observed under 
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certain circumstances’~*® for several fluorescent 
dyes. This phosphorescence is of the ordinary 
type in that its wave-length distribution is the 
same as that of the fluorescent light, and the 
probability of the transitions increases rapidly 
with temperature. On the other hand, there is 
(according to the Franck-Condon principle) a 
small but finite probability of a transition be- 
tween the states N’ and N. At low temperature, 
where transitions from N’ to E (induced by 
collisions of the first kind) are rare events, the 
direct transition, N’—N, becomes relatively im- 
portant. The experimental counterpart to this 
change is the appearance at low tempera- 
tures?1%! of a long-lived fluorescence whose 
wave-length maximum is shifted toward the red, 
relative to the ordinary fluorescent band. 

It has been customary to explain these optical 
phenomena in terms of metastable states. While 
metastable states undoubtedly occur in dye 
molecules which contain many resonating groups, 
the observed long lives of the phosphorescence as 
well as some of the kinetic data make it appear 
very probable that tautometric rather than 
metastable states are responsible for these phe- 
nomena. As was mentioned before, the time, ar, 
which elapses between the act of absorption and 
internal conversion is strongly affected by small 
changes in the internal motions of the molecule. 
These motions are determined only by the 
thermal (equilibrium) oscillations, since the ex- 
cess of oscillatory energy which the molecule 
acquires at the moment of absorption is quickly 
dissipated by collisions with solvent molecules. 
Lowering the temperature decreases the thermal 
energy of oscillation, and so decreases the 
probability of internal conversion and therefore 
favors direct fluorescence. This conclusion is 
substantiated by the familiar experimental ob- 
servation that the fluorescent yield of dye solu- 
tions increases as the temperature is lowered. 
While the present hypothesis leads directly to the 
conclusion that ar, and therefore the fluorescent 
yield, should be strongly influenced by inter- 

7 Wavilow and Lewschin, Zeits. f. Physik 35, 920 (1926). 
10 (Is and Vinokurov, Physik. Zeits. (USSR) 10, 

® Schischlowski and Wavilow, Physik. Zeits. (USSR) 5, 
379 (1934). 

© Kautsky and Hirsch, Ber. 65, 401 (1932). 


" Compare also Wick and Throop, J. Opt. Soc. Am. 25, 
318 (1935). 
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molecular forces, it is unable to predict in general 
the direction of this influence. However, experi- 
ments show that any circumstances, such as a 
viscous solvent or adsorption on a solid surface," 
which hinders the oscillatory motion of the dye 
molecule acts like a low temperature in increasing 
the intensity of fluorescence.” 

Internal conversion which leads to the forma- 
tion of a tautomer may profoundly influence the 
chemical behavior. Since the tautometric shift 
very probably involves a shift of a hydrogen 
atom, it is to be expected that in some cases the 
tautomer will contain a very labile or loosely 
bound hydrogen atom. If, in this way the energy 
necessary to dissociate the molecule into a radical 
and a hydrogen atom is reduced below 50 kcal.," 
an appropriate collision with an oxygen molecule 
will result in the formation of an HO; radical and 
a partially oxidized dye radical. 

Although the most important applications of 
the ideas presented here are to the photochemical’ 
reactions of dyes a detailed discussion of these 
complicated reactions is beyond the scope of this 
paper. A brief outline of the established facts of 
dye photochemistry and of their possible inter- 
pretation is presented as an appendix. There are, 
however, two groups of experiments which de- 
serve special attention. 

It is particularly noteworthy that the hy- 
pothesis of tautomer formation is able to explain 


Fic. 2. Schematic 
representation of 
the energy levels N’ 
and transitions cor- 
responding to a sys- 
tem where internal 
conversion leading 
to tautomer forma- 
tion is possible. 


12 Kautsky and Hirsch, Ber. 64, 2677 (1931). 

For a discussion of this fact see Lewis and Calvin, 
Chem. Rev. 25, 273 (1939). It is also possible that in some 
cases, particularly in solutions where the concentration of 
the dye is high, that viscous solvents may increase the 
fluorescent yield by hindering ‘‘self-quenching.” 

4 There appears to be no quantitative evidence in re- 
gard to the heats of dissociation of such tautomers. How- 
ever, there is evidence that the energy necessary to disso- 
ciate stable molecules into radicals and hydrogen atoms 
depends upon the structure of the molecule; see for example 
Taylor and Smith, J. Chem. Phys. 8, 543 (1940); and Bent 
and Cuthbertson, J. Am. Chem. Soc. 58, 170 (1936). 
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the fact that the fluorescent yield of chlorophyll 
in an inert organic solvent is only slightly 
changed'® by the addition of oxygen and a 
reducing agent, at such concentrations that the 
quantum yield of the chlorophyll-sensitized 
photo-oxidation is approximately unity. Since 
the yield of fluorescence remains constant and is 
of the order of a few percent, the process which 
limits the lifetime of the excited dye molecule 
must be the same regardless of whether photo- 
oxidation is or is not taking place. The observed 
quantum yield of the photosensitized oxidation 
of one indicates that the greater part of the 
energy of excitation of the dye molecule remains 
available for chemical purposes (sensitization) 
after the molecule has (practically) lost its 
ability to radiate. Both of these facts are con- 
sistent with the assumption that the chlorophyll 
molecule undergoes internal conversion resulting 
in the formation of a reactive tautomer, which is 
capable of initiating those secondary steps which 
are responsible for the observed chemical changes. 
A brief outline of these secondary steps and their 
relation to the available experimental facts is 
given in the appendix. 

A second group of experiments, performed by 
Kautsky,!® deserve special treatment; not only 
because of their importance in the theoretical 
development of dye photochemistry, but also 
because of the controversies which have existed 
in regard to their interpretation. In these experi- 
ments Kautsky prepared two samples of pow- 
dered silica gel; he allowed one to adsorb the 
fluorescent dye trypaflavine, and the other, the 
easily oxidizable leuco base of a photochemically 
inert dye. A mixture of these two dry powders 
was placed in an evacuated glass vessel. When 
the mixture was illuminated the trypaflavine was 
excited to fluorescence and phosphorescence, but 
there was no chemical action. In the presence of 
0.002 mm of QO: the phosphorescence was 
quenched and a noticeable oxidation of the leuco 


% Compare Franck and Levi, Zeits. f. physik. Chemie 
B27, 409 (1934). This experimental result has been con- 
firmed by the present authors (unpublished work) who 
found that the fluorescent yield of chlorophyll in a 10-° M 
solution saturated with air and containing allyl thiourea at 


a concentration of 0.5 M was about 80 percent (or 85 per- © 


cent) of the yield observed in a 10-* M solution in pure 
air-free solvent (acetone). 
16 Kautsky, de Bruijn, Neuwirth and Baumeister, Ber. 


66, 1588 (1933). 


FRANCK AND R. LIVINGSTON 


base to the corresponding dye could be recognized 
by the change in color. Oxygen pressures greater 
than 0.003 mm greatly retarded this oxidation. 

In an attempt to explain these results Kautsky 
assumed that an adsorbed dye molecule can exist 
in a long-lived metastable state which is capable 
of importing its energy to an Oz molecule (by a 
collision of the second kind) with the formation 
of a metastable '= Oz molecule. He also assumed 
that this '2 Oz molecule could act as a reactive 
oxidizing agent, the ‘‘oxygen carrier’ of the 
reaction. 

This mechanism was improved by Weiss,!” who 
suggested that the excited dye molecule would 
split off a hydrogen atom (either spontaneously 
or aided by a collision of the second kind) and 
that this atom could react with an oxygen 
molecule to form an HO, radical capable of 
acting as an “oxygen carrier.” This mechanism 
has the advantage that it eliminates the necessity 
of postulating the existence of the very improb- 
able long-lived metastable dye molecule and of 
extraordinarily efficient deactivating collisions 
between normal and metastable oxygen mole- 
cules.!8 However, it has the disadvantage that it 
assumes that a dye molecule which has been 
excited by the absorption of a quantum of visible 
light is capable of spontaneously splitting off a 
hydrogen atom, a process which is energetically 
impossible. 

These several difficulties may be avoided by 
assuming that the excited dye molecule (RH*) 
transforms spontaneously into a reactive tau- 
tomer (HR). This tautomer may revert either to 
the original dye molecule (RH) or, more rarely, 
to the excited molecule (RH*) ; the latter process 
accounting for the faint phosphorescence which 
Kautsky observed. In the presence of oxygen the 
following reaction is probable. 


The HO; radical, so formed, may either diffuse to 
and be reduced by the leuco base or be lost by the 
back reaction 


HO:+R—-RH-+0:. 


17 { . Weiss, Naturwiss. 35, 610 (1935). 

18 Tt is also very doubtful if 12 O2 molecules are ever of 
importance in photo-oxidations; Compare Gaffron, Ber. 
68, 1409 (1935); Biochem. Zeits. 264, 251 (1933); Zeits. 
f. physik. Chemie B37, 437 (1937); Rollefson, Chem. Rev. 
17, 427 (1935). 
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PHOTOCHEMISTRY OF DYESTUFFS 


TABLE I. Pholochemical properties of typical dyes. 


PaRARO-| Fiuo- | METH- 
SANA- RES- YLENE | CHLORO- 
DYE LINE CEIN BLUE PHYLL 

Photosensitizes oxidations - + + + 
Fluorescence quenched by | no fluo- + + - 

reducing agents rescence 
Reversibly photo-bleached - + ad 

by reducing agents 
Reversibly photo-bleached (-) (-) + 

by oxidizing agents 


This mechanism is also consistent with the fact 
that the efficiency of the oxidation of the leuco 
base first increases to a maximum and then falls 
off rapidly as the pressure of oxygen is increased. 
The rate of oxidation of the leuco base is directly 
proportional to the rate of escape of HO: radicals 
from the capillary spaces in which the fluorescent 
dye molecules are adsorbed on the “inner 
surfaces” of the silica gel granules. To a first 
approximation this rate of escape is directly 
proportional to the concentration of HO: radicals 
in the capillaries and inversely proportional to 
the total gas pressure. Instead of escaping from 
the capillaries, the HO: radicals may be lost by 
reacting with the partially oxidized dye mole- 
cules. High total pressure would favor this 
reaction relative to the escape from the capillary, 
if the rate of the reaction were practically inde- 
pendent of the pressure, as would be the case if 
adsorbed HO; radicals were free to migrate upon 
the silica gel surface. It is also possible that the 
HO: radicals disappear within the capillaries by a 
bimolecular process of decomposition. However, 
if this latter reaction were relatively important a 
noticeable photobleaching of the fluorescent dye 
should occur in the presence of Oz, a process 
which Kautsky did not observe. 


APPENDIX: AN OUTLINE OF THE INTERPRETATION OF THE 
EXPERIMENTAL FACTS OF PHOTOCHEMISTRY OF DYES 
IN TERMS OF THE PRECEDING DISCUSSION 


The observations on the photochemistry of dyes and 
pigments can be interpreted in a manner which is con- 
sistent with the foregoing discussion of the optical prop- 
erties of complex molecules. Published studies of the 
photochemistry of dye solutions include the investigation of 
photosensitized oxidation-reduction reactions, of photo- 
bleaching of dyes by reducing agents, and of (at least in the 
case of chlorophyll) rapid-reversible photo-bleaching in the 
absence of a substrate. 

Since the various dyes differ in their behavior in respect 
to these reactions, it is convenient to divide them into 
several classes. The majority of dyes are photochemically 
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inactive and practically nonfluorescent. Pararosanaline 
may be taken as representative of this group. Fluorescein 
and the related dyes, eosine, rose bengale, and rhodamine, 
constitute a second group. Methylene blue is typical of a 
third group which includes thianine and, presumably, other 
thiazine dyes. The fourth group consists of chlorophyll and 
possibly the porphyrines. In Table I, which summarizes 
some of the photochemical characteristics of these types of 
dyes, a minus sign (—) signifies that the indicated process 
either does not occur or occurs only very inefficiently, while 
a plus sign (+-) signifies that the process occurs efficiently. 

The photo-bleaching referred to in the fifth row is the 
reversible bleaching produced by oxidizing agents such as 
Fe+++, This is, of course, quite different from the slow 
irreversible bleaching which solutions of all dyes undergo 
when they are illuminated in the presence of oxygen. 

To conserve space all of the reaction steps which appear 
to be required for any of the photochemical reactions of the 
several classes of dyes are grouped together in the following 
list. It should be remembered, however, that the explana- 
tion of no single set of experiments requires more than a 
small part of this rather formidable array of reaction steps. 


1 DH + hy — DH* 

2 DH* DH + hy, (fluorescence) 

3 DH* — HD (internal conversion ac- 
companied by tau- 
tomerization) 

4 DH* — DH (internal conversion fol- 
lowed by collisions 
with the solvent mole- 
cules) 

50: + DH*~D +HO0O, 

6 AH: DH* DH: + AH 

7 HD — DH* (reversion, responsible 
for phosphorescence) 

8 DH + HD — DH: + D 

9 AH: + HD — DH, + AH 

1000 +HD-—-D_ 

11 2DH: = D~ + DH; (reversible, semiquinone 
equilibrium) 

12D + DH: ~2DH 

13 O: + DH: ~ DH + HO, 

14 AH + DH: = DH; +A (reversible) 


15 AH + DH: ~ DH + AH; 

16 AH: + D DH +AH 

17 HO.+ D DH +0, 

18 AH; + HO: ~ AH + 

19 AH + H.,02— Oxidation products 
20 AH +O: ~A_  +HOs.. 


The special symbols used in this scheme have the 
following significance: DH—dye, DH*—electronically ex- 
cited dye, HD—tautomer of dye, DH2—semiquinone, 
D—partially oxidized form of dye, DH;—leuco base, 
AH,—reducing agent, A and AH—oxidized or partially 
oxidized forms of the reducing agent, O—an oxidizing agent 
and Re—its reduced form. 

Steps 5, 8, 10, 12, 16, and 17 appear to be required only 
for reactions involving chlorophyll. Since steps 8, 12, and 17 
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are somewhat unusual it may be worthwhile to point out 
that these steps are energetically feasible. Neither the heat 
of dissociation of the labile hydrogen of chlorophyll, nor the 
heat of association of normal chlorophyll with atomic 
hydrogen are known. We may consider, therefore, that the 
mechanism is a probable one if it is possible to assign 
plausible values to these energy terms which make steps 
8, 12, 13 and 17 exothermic. It is very probable that the 
heat of dissociation is less than 90 cal.!° and the heat of 
association less than 53 cal.,2° which is a probable value for 
the hydrogen atom affinity of O2. There is a fairly wide 
range of values which satisfy these conditions. For ex- 
ample, if the heat of association is taken as +45 cal. and 
the heat of dissociation as — 80 cal. the heats of reactions 
8, 12, 13, and 17 will be +3, +35, +8, and +27 cal., 
respectively. 


Quenching of fluorescence 


The quenching of the fluorescence of the dyes of the 
thiazine and fluorescein types is caused chiefly by step 6; 
although in some cases the added substance may, by 
complex formation or otherwise, increase the probabilities 
of steps 3 and 4 relative to 2. 

In the case of the nonfluorescent dyes, the probability of 
the occurrence of either step 3 or step 4 must be much 
greater than that of step 2. The fact that these dyes are 
usually very poor sensitizers indicates that if (by step 3) 
they form the tautomer, HD, it is either energetically very 
close to the normal molecule or else it is very short lived. 
There are, however, a few ‘“‘nonfluorescent” dyes which 
noticeably, although rather inefficiently, sensitize oxidation 
reduction reactions. It would not be at all surprising if 
further investigation were to reveal the existence of other 
nonfluorescent dyes which are efficient sensitizers. 


Photo-bleaching 


The bleaching of fluorescein and thiazine dyes which 
occurs when a solution containing a reducing agent is 
illuminated, is due to step 6 (or steps 3 and 9) followed by 
steps 11, 14, and 15. When little or no bleaching occurs, as 
in the case of the fluorescein dyes,”! step 15 is much faster 
than step 14. The reverse is true when bleaching occurs very 
readily, as in solutions containing methylene blue and allyl 
thiourea.2 The regeneration of the color in the dark, as 
with methylene blue and phenyl hydrazine,” probably is 
caused by the reverse step 14 followed by 15. Rapid 
reversible bleaching which is proportional to the square 
root of the intensity of the absorbed light, such as occurs in 
solutions containing thianine and KI,** probably involves 
steps 6 and 15, or 3, 9, and 15. 

The irreversible photo-bleaching of chlorophyll by re- 
ducing agents is a very inefficient process. If it occurs at all 


19 Rollefson and Burton, Photochemistry (Prentice-Hall, 
1939), p. 422. 

20 Bray, J. Am. Chem. Soc. 60, 82 (1937). 

21 Livingston and Hurd, J. Phys. Chem. 45 (1941), in 
press. 

22 Mudiovcic, Zeits. f. wiss. Phot. 26, 171 (1928). 

3 Holst, Zeits. f. physik. Chemie A179, 172 (1937). 

*% Weiss, Trans. Faraday Soc. 35, 48 (1939). 


it may be caused by an inefficient process similar to step 11. 
The reversible photo-bleaching by oxidizing agents such as 
Fet**% jis due to steps 3, 10, and 16. The oxidized form, D 
is yellow and photolabile, and probably differs fom the 
oxidized form produced by step 8. When the oxidizing 
agent is Fe+**, steps 10 and 16 may be rewritten as follows: 


10’ 
16’ Fe(H,O)*++D —DH-+Fe(OH)**. 


The relatively slow bleaching of chlorophyll by Fe*** in 
the dark probably involves a slow thermal reverse of step 7 
followed by 10 and 16. The rapid reversible bleaching which 
occurs in the absence of a substrate?* is caused by steps 3, 8 
and 12. Since the bleached chlorophyll is regenerated by a 
second-order process, we may conclude that the color of the 
reactive tautomer (HD) is practically the same as that of 
the normal molecule and that the actual bleaching is a 
consequence of a bimolecular process, probably step 8. The 
inhibition of this process by O2 is due to 13 and 17. Since 
this action of O2 probably does not reduce the steady state 
concentration of D, this substance must have about the 
same color as the normal molecule DH. The bleached form 
must be DHz. The inhibition by Fe*+t is due to steps 16 
and 15. 


Sensitized photo-oxidations 


Photo-oxidations sensitized by fluorescein and thiazine 
dyes!21.27.28 are the results of steps 6 (or 3 and 9), 13, 
18, and 19. The sensitized oxidation of benzidine by 
chlorophyll***° probably follows a similar course, but the 
chlorophyll sensitized oxidation of substrates like allyl 
thiourea and isoamylamine* is due to steps 3 and 8 (or 9), 
13 and follow reactions such as 18, 19, and 20. 


Destructive oxidation of a dye 


Photo-bleaching in the presence of oxygen (i.e., de- 
structive oxidation) of these dyes (as well as chlorophyll) is 
greatly accelerated by the presence of traces of oxidizable 
compounds.” This action is probably due to steps 6 (or 3 
and 9) and 13, followed by a relatively inefficient reaction 
between the dye molecule and the HOz radical. The 
inefficiency of this latter process is indicated by the strong 
protective action of moderate concentrations of the same 


reducing agents which induce the bleaching when present as 


traces.*8 


— and Weiss, Proc. Roy. Soc. A162, 251 

2% Porret and Rabinowitch, Nature 140, 321 (1937). 

27 Schneider, Zeits. f. physik. Chemie B28, 311 (1935). 

28 Spealman and Blum, University of California, Pub. in 
Physiol. 8, 147 (1937). 

29 Noak, Biochem. Zeits. 183, 153 (1927). 
as and Levi, Zeits. f. physik. Chemie B27, 409 
' 81 Gaffron, Ber. 60B, 755, 2229 (1927). 

32 Weber, Ber. 69, 1026 (1936). 

33 The possibility of a direct reaction between DH*, or 
HD, and O2 cannot be excluded. Indeed, some such 
mechanism is probably necessary to explain the high yields 
of the photo-bleaching of the cyanine dyes by oxygen. 
Bhattacharyya and Dahr, J. Ind. Chem. Soc. 6, 145 (1929). 


F 
| 
tk 
| be 
te 
tic 
wl 
he 
| en 
dy 
the 
dei 
ma 
pu: 
suc 
der 
are 
| tha 
stat 
stal 
pur 
sucl 
thei 
heai 
| | resp 
defi 
the 
| cons 
| the | 
line 
| | 
| This 
integ 
(exce 
a pat 
Now 
are us 
| 
if 
| 


FEBRUARY, 1941 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 9 


Notation for the Derivatives of the Two Types of Line Integral in Thermodynamics 


G. TUNELL 
Geophysical Laboratory, Carnegie Institution of Washington, Washington, D. C. 


(Received October 19, 1940) 


N thermodynamics, quantities of two types 

are constantly dealt with, namely: quantities 
that depend only on the momentary state of the 
body or system (such as pressure, volume, 
temperature, energy, and entropy), and quanti- 
ties that depend on the series of states through 
which the body or system passes (viz., work and 
heat). On the one hand, energy and entropy 
enter the mathematical formulations of thermo- 
dynamics as line integrals that depend only on 
the initial and final values of the variables 
defining the state of the system; standard mathe- 
matical symbols altogether adequate for the 
purposes of thermodynamics are available for 
such line integrals,! and also for their partial 
derivatives. On the other hand, work and heat 
are defined in thermodynamics as line integrals 
that depend on the path (i.e., the entire series of 
states through which the body or system passes) ; 
standard mathematical symbols adequate for the 
purposes of thermodynamics are available for 
such line integrals, but the same is not true of 
their derivatives. The derivatives of the work and 
heat line integrals are total derivatives with 
respect to the variable chosen as the parameter 
defining the path of the integral. For example, if 
the path (of a homogeneous fluid system of 
constant mass and composition) be defined by 
the equation p=f(t), the derivative of the work 
line integral is given by the following equation 


This equation for the derivative of the work line 
integral holds for any path in the 7, p plane 
(except a straight line parallel to the p axis). For 
a path parallel to the T axis in this plane the 
equation reduces to 


Now in thermodynamics many siieiods planes 
are used (such as the #, V plane, the T, V plane, 


‘Cf. G. Tunell, J. Phys. Chem. 36, 1744-1771 (1932). 


(1) 


(2) 


the E, S plane, the T, S plane, etc.) and a concise 
but unambiguous notation for the derivatives of 
the work and heat line integrals along paths 
parallel to the coordinate axes in each of these 
planes has been very much needed. In order to 
understand Eq. (2) it is necessary to have the 
supplementary information, not indicated in the 
equation itself, that the work line integral is 
being differentiated along a path in the 7, p plane 
parallel to the T axis; the derivative of the work 
line integral along a line parallel to the T axis in 
the T, S plane is not equal to p(dV/dT), but it 
would still be represented by the same symbol 
dW/dT in the notation of Eq. (2). A symbol that 
would convey the necessary supplementary in- 
formation is (dW/dT), resp. (dW/dT)s. The 
symbol (dz/dx), was formerly used for the partial 
derivative of a function of two independent 
variables, and for this reason the writer previ- 
ously has hesitated to suggest that it be used to 
fill the need in thermodynamics under discussion. 
However, since Legendre and Jacobi introduced 
the symbol dz/dx for the partial derivative of a 
function of two or more independent variables, 
the symbol (dz/dx), has fallen into disuse for this 
purpose, and it has seemed to the writer for some 
time that it might well be used in thermodynamics 
with the new meaning mentioned above. The 
reason for writing this article is that in Eucken’s 
Grundriss der physikalischen Chemie? and in 
Sherwood and Reed’s Applied Mathematics in 
Chemical Engineering? the symbol (dQ/dT)y is 

*Eucken, Grundriss der physikalischen Chemie (Aka- 


demische Verlagsgesellschaft M.B.H., Leipzig, 1934). In 
the fourth edition (1934) Eucken writes, on page 88, the 


aT. 
The same usage also occurs in earlier editions. 

3 Sherwood and Reed, Applied Mathematics in Chemical 
Engineering (McGraw-Hill, 1939), p. 163. Sherwood and 
Reed write the equations 


C= 77), 


dQ= (59) at+(2), dv=CadT +hdv. 


and 
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dQ aU 
dw OV 
|_| 
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used precisely to represent the derivative along a 
path parallel to the T axis in the T, V plane of 
the heat line integral, and it seems desirable to 
urge that the symbol with straight d’s, paren- 
theses, and subscript be accepted generally in 
this sense. Eucken and Sherwood and Reed have 
introduced this symbol in this sense so unob- 
trusively that its significance may escape general 
notice. 

The need in thermodynamics for such a 
symbol as (dQ/dT), is emphasized by a compari- 
son of certain second derivatives. In the case of a 
function of the state of the body there exists a 
well-known relation between the cross deriva- 
tives, which in the case of the energy E is as 
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but the partial derivative with respect to T of the 
total derivative of Q with respect to p along a 
line parallel to the p axis is not equal to the 
partial derivative with respect to p of the total 
derivative of Q with respect to T along a line 
parallel to the T axis. The symbol dQ/dT is also 
inadequate because the inequality would then be 


written 
( ~) ( dQ 
aT dp), \apdT/r 
which is quite ambiguous without supplementary 
information in the text. Actually the derivatives 


along paths parallel to the coordinate axes of the 
line integrals that depend on the path (work and 


follows: heat) are themselves functions of all the inde- 
pendent variables that determine the state of the 
body or system‘ and admit partial derivatives 


with respect to each of these variables. This fact 


VE 
aTap apaT 


(3) 


whereas, with the functions that depend on the _ js well expressed by the notation ma 
series of states that the body passes through, one dey 
has in place of the preceding equation an ha dQ | He 
inequality, which in the case of the heat Q is as aT\ dp] rly be 
‘ follows: poi 
; It may be suggested finally that confusion be- the 
| h following Eucken and re 
|—(-) —(—) | (A) tween the new usage (following Eucken an 
dT\dpJ/ Lop\dT/ Sherwood and Reed) of the symbol with straight ‘ 
, d’s, parentheses, and subscript, and the old 
id A consideration of the preceding equation and usage, which has been discarded since the time of bet 
: inequality shows at once that neither the symbol Jacobi, is very unlikely. In papers or books in “n 
r (9Q/8T)» nor the symbol dQ/ aT is suitable for hich this symbol is used in the new sense there at 0 
4 the purposes of thermodynamics. The symbol will also occur derivatives of functions of two or st 
\ (8Q/8T), cannot be used because according to more variables, which will of course be written ron 
J accepted mathematical usage at the present time with rounded @’s. Thus where both types of pa 
* would stand for the partial derivative of @ symbol are found in one publication it will be ness 
function Q of the two independent variables T Clear that the symbol with straight d’s, paren- tion, 
and p with respect to 7, whereas Q is not a theses, and subscript is used in the sense of a ‘ “s 
function of T and 9; ar, Cae would im- derivative along a line parallel to a coordinate hae 
mediately expect on seeing the symbol (8Q/8T)» axis of a line integral that depends on the path. In 
that an equation could be written between aueidiaanis pron 
: 4 A point well stated and 1 hasized by Clau- : 
ig —— Section 2, Part 9 (Friedrich Vieweg und Sohn, Braunsch- nd 
weig, 1867), p. 7. 
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these are terse and contain few figures, and especially 
few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
Contributions to this section should not exceed 600 words in 
length and must reach the office of the Managing Editor not 
later than the 15th of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
authors. The usual publication charge ($3.00 per page) will 
not be made and no reprints will be furnished free. 


Influence of the Temperature of Crystallization 
on the Melting of Crystalline Rubber 


NORMAN BEKKEDAHL AND LAWRENCE A. Woop 


National Bureau of Standards, Washington, D. C. 
January 8, 1941 


Data presented in this communication show an instance 
in which melting of a crystalline material is very much 
dependent on the temperature at which the crystals have 
been formed. It is well known that many substances in 
which crystallization is relatively slow can be crystallized 
at different temperatures in a range below the melting 
point, but no effect of the crystallization temperature on 
the temperature of melting seems to have been previously 
reported. The quantitative results for crystalline rubber, 
the material under investigation, are shown in Fig. 1. 

The crystallization’ of unvulcanized rubber in the un- 
stretched state has been found to occur at temperatures 
between about —40°C and 13°C. The time required for 
crystallization is about one year at 13°C, about ten days 
at 0°C, and a few hours at — 20°C. Below —40°C the mo- 
bility is presumably insufficient for the formation of crys- 
tals. Crystallization and fusion are accompanied by changes 
in volume, heat capacity, light absorption, birefringence, 
x-ray diffraction, and mechanical properties such as hard- 
ness. The volume decreases about 2.5 percent on crystalliza- 
tion, and the magnitude of the change is little influenced, 
if at all, by the temperature. Fusion, as measured by the 
volume increase, is found to be independent of the rate of 
heating, and to occur over a range of five or ten degrees. 

In the present work, volume change as measured in a 
mercury-filled dilatometer was taken as the criterion of 
crystallization. Smoked sheet rubber was held at various 
fixed temperatures until no further volume change could 
be detected. After each crystallization the temperature was 
raised, and observations were taken for a graph of the 
volume as a function of temperature. The temperature 
corresponding to the beginning of melting could be located 
as that at which the observations began to deviate from 
the curve representing the normal thermal expansion of 
crystalline rubber. Likewise, at the end of the melting 
range the observations began to follow the curve repre- 
senting the normal thermal expansion of amorphous rubber. 
These two limiting temperatures could be located within 
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about one-half a degree and are plotted in Fig. 1. The 
shaded area between the curves represents the melting 
range. 

It is to be observed from the figure that the higher 
the temperature of crystallization, the higher are the tem- 
peratures at which both the beginning and ending of 
melting occur. The range of melting itself appears to be- 
come somewhat narrower at the higher crystallization tem- 
peratures. From the fact that the broken curve does not 
intersect or touch the shaded area, it appears that there 
is no temperature at which crystalline and amorphous 
rubber are in equilibrium at atmospheric pressure. If the 
crystallization occurs between about 6° and 10°C, the 
melting of the crystals begins at a temperature about 4° 
above that at which the crystals were formed. This in- 
terval is greater than 4° for crystallization at temperatures 
above or below this range. In previous work in this labora- 
tory a melting range from 6° to 16°C has been reported! for 
crystallization temperatures near 0°C. Van Rossem and 
Lotichius? found that crystals formed at —10°C melted 
between 2° and 14°C. Both these sets of observations are 
included in Fig. 1. 


TEMPERATURE, °C 


-20 -/0 5 10 4S 


TEMPERATURE OF CRYSTALLIZATION, °C 


Fic. 1. Melting range of crystalline rubber as a function of the tem- 
perature of crystallization. The lower points represent the beginning of 
melting, the upper points the ending of melting. The broken curve 

the crystallization temperature as ordinate for easier comparison 
= = range of melting temperature in each case. Van 
ichius. 


m and 


The results here presented undoubtedly explain many 
previously reported discrepancies in the temperatures of 
melting of crystalline rubber. They may also be expected 
to furnish an explanation of certain effects with stretched 
rubber recently observed by Treloar,? who, lacking the 
present data, felt that he was not justified in making the 
assumption that ‘the melting point of the crystals could 
depend on the temperature at which they are formed.” 

There seems to be no reason why the effect discussed 
here should be limited to rubber, and, if investigations are 
made, this effect may be found in many other organic 
compounds of high molecular weight. 

1N. Bekkedahl, J. Research Nat. Bur. Suns. 13, 411 (1934) RP717. 
Reprinted in Rubber Chem. Tech. 8, 5 (19. 

2A. Van and J. Lotichius,  eceachoede 5, 2 (1929). Reprinted 
in Rubber Chem. "Tech. 2, 378 (1929). 


R. G. Treloar, Trans. Fa raday Soc. 36, 538 (1940). Reprinted in 
Rubber Chem. Tech .13, 795 (1940). 
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194 LETTERS TO THE EDITOR 


The Reaction Between Hydrogen and Oxygen: 
The Upper Explosion Limit and the Reac- 
tion in Its Vicinity! 

GUENTHER VON ELBE, Coal Research Laboratory, Carnegie Institute of 
Technology, Pitisburgh, Pennsylvania 


AND 


BERNARD LEwis, Physical Chemistry Section, Central Experiment Stalion, 
Bureau of Mines, Pittsburgh, Pennsylvania 


November 23, 1940 


We have found that water vapor effects a marked lower- 
ing of the upper explosion limit of hydrogen and oxygen. 
This is illustrated by the following typical results for a 
2 : 1 mixture at 530°C, H.O being admitted to the reaction 
vessel from a reservoir. 


Percent H,O 0 3.35 6.37 11.03 
Upper limit, mm Hg 85.6 54.9 41.8 29.4 


Nine determinations with water percentages ranging from 
3 to 11 gave an average value of kgn,o0/ken, equal to 14.3, 
with a maximum error of 0.7. These results demonstrate 
that the water effect is no different in kind than the influ- 
ence of other gases as described by the explosion condition 


2ke = ken, (Hz) +h60,(O2) +kex(X), 


where X is water or any other inert gas. This equation 
corresponds to the reaction scheme 


1. OH+H.=H,0+H 

2. H+0.=OH+0 

3. O+H:=OH+H 

6. O., X) =HO2+(H2, Oz, X). 


The equation numbers follow those given elsewhere.? 

The ratio kgo,/ken, was determined at 530°C as 0.35 
+0.01 by observing the upper limit for mixtures of hydro- 
gen and oxygen containing 20, 33.3, 50, and 66.7 percent 
hydrogen. The value showed no dependence on mixture 
composition. Similarly, k¢x/keu, for other gases were found 
to be independent of the percentage present; for He, A, 
N. and CO; they are 0.36, 0.20, 0.43 and 1.47, respectively. 
The values for Os, He, A and Ne agree well with former 
determinations.’ 

On slowly approaching the upper limit from the high 
pressure side one may observe, very close to the limit, 
short bursts of reaction that die out rapidly, as is illus- 
trated in Table I for a mixture containing 20 percent He 
at 530°C. 


TABLE I. Initial pressure =137.7 mm (upper explosion limit slightly 
below this). 


37 100 


Time, sec. 0 15 23 30 
—dO2/di, mm/min. 0.93 0.80 0.40 0.40 0.05 


The water thus formed pushes the limit toward lower 
pressures. That the water alone is responsible for this 
shift was shown by the precise agreement between the 
observed shift after allowing a known amount of water to 
be formed by reaction, and that calculated from the 
efficiency of water of 14.3. 

A typical curve of the change of reaction rate with 
pressure above the upper explosion limit is shown in Fig. 1. 


The Pyrex spherical reaction vessel, 7.2 cm in diameter, 
was heavily coated with KCl. The rates proved to be very 
well reproducible over a period of many weeks. They pass 
through a minimum, as expected. On the left side of the 
minimum the reaction inhibits itself (see above), and only 
initial rates are given. 

We have also carried out upper-limit determinations in 
the same vessel from 400°C to 570°C. The values are con- 
siderably higher than those found by other investigators.‘ 
This and the occasional poor reproducibility found by 
these investigators may be explained by the water-vapor 
effect. 

The foregoing observations furnish an explanation of 
most of the questions raised by Oldenberg and Sommers.® 
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RATE OF PRESSURE DECREASE, MILLIMETERS PER MINUTE 


40 120 200 360 
PRESSURE OF 2H2~Oz MIXTURE, MILLIMETERS OF MERCURY 


Fic. 1. Reaction rate of mixtures of 2H2+Oz2 at 530°C in spherical 
KCl-coated Pyrex vessel, 7.2 cm in diameter. @ indicates duplicated 
observations. 


The latter, working at 568°C, were unable to observe a 
minimum rate on approaching the upper limit. However, 
at this temperature the rates near the minimum are some 
100 times greater than at 530°C, thus tending to obscure 
the branch of the curve to the left of the minimum by 
self-inhibition. Self-heating of the mixture, even of the 
order of 1°, should also contribute to this because the 
coefficient k2 increases sharply with temperature; therefore, 
the transition between slow and explosive reaction becomes 
very abrupt. Oldenberg and Sommers’ observation that 
the limit shifted from about 180 mm to about 120 mm after 
4 percent of the mixture reacted is in striking agreement 
with the value of 14.3 for the water-vapor efficiency, the 
calculated shift being identical. It is also evident that the 
increase of the upper limit on coating the vessel with KCl, 
which these authors report, is caused by the decreased rate 
of water formation. The same limit can be found in both 
vessels by operating sufficiently fast. 

We are making a detailed study of the kinetics of the 
hydrogen-oxygen reaction and shall publish further details 
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later, including a discussion of the remaining points raised 
by Oldenberg and Sommers. 


1 Published by permission of the Director, Bureau of Mines, U. S. 
Department of the Interior, and the Director, Coal Research Labora- 
tory, Carnegie Institute of Technology. 

2 Lewis and von Elbe, a Flames and Explosions of Gases 
Ca Univ. Press, 1938), p. 3 

. H. Grant and C. N. Hidshelwood, Proc. Roy. Soc. (London) 
AlaL “29 A aad A. A. Frost and H. N. Alyea. J. ie. Chem. Soc. 55, 


3227 (193. 

4H. , — and C.N. Hinshelwood, Proc. Roy. Soc. (London) 
A122, 610 (1929); A. A. Frost and H. N. Alyea, reference 3; G. H. Grant 
and C. N. Hinshelwood, ibid. 

5Q. Oldenberg and H. S. Sommers, J. Chem. Phys. 8, 468 (1940). 


Coriolis Perturbations and Molecular 
Dimensions in Germane GeH, 


Ta-You Wu 


Department of Physics, National — Ass. University, Kun-ming, 
ina 


November 25, 1940 


In a recent paper Murphy! has calculated the structure 
of the v4 fundamental in SiH, and in GeH, by evaluating 
the matrix elements of the coriolis perturbations between 
v, and the doubly degenerate v2. It was found that these 
perturbations satisfactorily explain the observed com- 
plexities in these bands and in particular the unsymmetrical 
shapes of the Q branches. While the agreement between 
the observed and the calculated structure is good in the 
case of SiH,, it is less so in the case of GeHy. The purpose 
of this note is to suggest that the poorer agreement is 
probably due to the use of a vlaue for the moment of 
inertia of GeH, which is highly questionable. 

Steward and Nielsen estimated J at 7.0X10~ g cm? 
for GeH, from the spacings Av3= 5.56, Ava~6 cm™ accord- 
ing to the theory of Teller and of Johnston and Dennison. 
This value of J yields the value 1.27A for the Ge-H distance 
which has been pointed out by the writer to be unreason- 
able on comparison with the Si-H distance in SiH,.2 On 
account of the complexities in the structure of v4 due to 
the perturbation by v2, the Avs has no strict meaning and 
must hence not be used in determining J. On the other 
hand, the assumption of pure valence and deformation 
forces in such molecules as CH,, SiH4, GeH, is probably a 
good approximation. On this assumption and from the 
observed value of Av; in the v3 band which is not perturbed 
by v2, we obtain the following dimensions (Table I). 
These values of I are also only approximate since the in- 
ternal angular momenta £’s obtained on the assumption of 
pure valence force system are approximate, as can be seen 
in the case of CH, in which the C—H distance is known 
to be 1.09A. Their relative values are, however, more 
reasonable than those obtained by employing Av; and Av, 
since one would expect the Ge-H distance to be greater 
than Si-H. It is hoped that the use of I9.8X 10 g cm 


TABLE I. 
From Av3 AND Av4 
(1 —£s)h 
Avs 4x*Avs | I X-H 
CH, 9.77 | 0.100 | 0.400 |5.09 X1074° | 1.07A |/5.47 X1074| 1.11A 
SiH, | 5.65 | 0.046 | 0.454 1.45 8.9 1.41 
GeH, | 5.56 | 0.018 | 0.482 |9.77 1.49 7.0 
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for GeH, in the perturbation calculations will improve the 
agreement between the calculated and the observed struc- 
ture of v4. 


1G. M. Murphy, J. Chem. Phys. 8, 71 (1940). 
2T.-Y. Wu, Vibrational Spectra and Structure of Polyatomic Molecules 
(National University of Peking), pp. 230-231. 


Conductivity of Calcium Salts 


ARTHUR S. JENSEN, Physics Department, University of Pennsylvania, 
Philadelphia, Pennsylvania 


AND 


M. RutH UNaAnGst, Pharmacology Department, Hahnemann Medical 
College, Philadelphia, Pennsylvania 


December 16, 1940 


It is of interest to the medical world to prevent convul- 
sive effects that usually follow the administration of certain 
anesthetics. Drs. R. Beutner and G. P. Miley,’ ? Depart- 
ment of Pharmacology, Hahnemann Medical College, 
observed that the simultaneous injection of certain calcium 
salts with these anesthetics inhibited convulsions. Dr. 
Beutner*® suggested that the action of the salt could be 
correlated with its calcium ion content, that a measure- 
ment of the salt’s electrical conductivity would be of 
interest to the formulation of a theory of the salt’s physio- 
logical action. 

Since there were found in the literature almost no data 
on the electrical conductivities of the salts used in the 
animal experiments, we measured them, using the stand- 
ard Kohlrausch bridge method with 60-cycle current and 
a conductivity cell with Pt electrodes, at 18°C, and the 
following results (Table I) were obtained: 


TABLE I. 


@ IN MHO/CM 
0.01 N 


SALT 0.1 N 0.001 N 

Ca Gluconate 2.52 X1073 4.97 X1074 6.92 X1075 
Ca Benzoate 2.68 7.25 8.46 

Ca Lactate 3.79 6.40 9.50 

Ca Levulinate 4.75 6.55 8.13 

Ca Salicylate 5.35 7.50 9.49 

Ca Chloride 8.62 10.1 12.5 

Ca Bromide 9.10 10.7 10.1 


Calcium gluconate held special interest to the pharma- 
cologist, so further detailed measurements were made on 
it (Table IT). 


TABLE II. 
CONCENTRATION AVERAGE CONDUCTIVITY 

18°C 0.1 N 2.52 X1073 mho/cm 

0.05 N 1.62 X10-3 

0.01 N 4.97 X1074 = 

0.005 N 2.75 X1074 

0.001 N 6.92 X10-5 
25°C 0.1 N 3.24 ot mho/cm 

0.01 N 6.50 X10 


Deviation for these readings was within +3 percent. 
These measurements were made in the Pharmacology 
Department of the Hahnemann Medical College. 


(1938). Beutner and G. P. Miley, Proc. Soc. Exp. Biol. Med. 38, 279 
(1935). Beutner and G. P. Miley, Proc. Soc. Exp. Biol. Med. 42, 547 
3 R. Beutner, Am. J. Physiol. 129, 489 (1940). 
4 Curr. Researches in Anes. and Analgesia 19, 132 (1940). 


